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ABSTRACT 


This  dissertation  is  a  study  of  the  description  and  effects 
of  particle  interactions  in  ionized  gases.  The  principal  results  are: 

(i)  An  expansion  theorem  for  the  linearized  Fokker-Planck  collision 
operator  for  each  component  of  a  two- component  fully  ionized  gas,  and 

(ii)  A  description  of  photon  scattering  from  a  partially  ionized  gas. 

It  is  shown  that  the  Fokker-Hlanck  collision  operator  generates 
a  complete,  continuous  set  of  velocity-space  eigenfunctions,  for  which 
high-speed  asymptotic  forms  are  found.  Since  the  set  is  continuous, 
the  expansion  formula  has  the  form  of  a  generalized  Fourier  integral. 

The  effect  of  neutral  atoms  on  the  spectrum  of  photons 
scattered  by  electrons  in  a  partially  ionized  gas  is  shown  to  be 
primarily  a  reduction  in  height  and  increase  in  width  of  the  two 
electron-plasma  "wings."  The  scattered  photon  spectrum  is  described 
for  several  characteristic  cases. 


BLANK  PAGE 


I.  INTRODUCTION 


The  purpose  of  this  dissertation  is  to  investigate  certain 
aspects  of  the  kiretic  behavior  of  ionized  gases.  The  emphasis  here 
is  primarily  on  the  description  and  effects  of  particle  interactions. 
Two  principal  results  are  obtained.  The  first,  developed  in  Chapter 
III,  is  an  expansion  theorem  for  the  linearized  Fokkei  Planck  collision 
operator  for  a  two  component  fully  ionized  gas.  Tie  second,  developed 
in  Chapter  IV,  is  a  description  of  photon  scattering  from  a  partially 
ionized  gas. 

The  function  of  Chapter  I  is  twofold:  it  presents  a  brief 
discussion  of  the  present  state  of  the  art  in  the  treatment  of  classi¬ 
cal  many-particle  systems,  thereby  setting  the  stage  for  the  rest  of 
the  work;  in  addition  it  contains  an  outline  of  this  work  designed  to 
guide  the  reader.  Chapter  II  contains  a  survey  o!-  some  of  the  more 
pertinent  results  of  plasma  kinetic  theory. 

The  reader  interested  primarily  in  the  results  of  this  disser¬ 
tation  may  choose  to  skip  the  first  two  chapters  and  proceed  directly 
to  Chapters  III  and  IV,  which  are  essentially  self-contained. 

1.  The  Ilany-Particle  Problem 

The  theory  of  classical  many-particle  systems  may  be  studied 
from  three  points  of  view.  One  may  begin  with  the  macroscopic,  or 
fluid,  equations  with  parameters  such  as  density,  mass  velocity,  and 
temperature  as  independent  variables,  and  involving  various  transport 
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coefficients,  e.g.,  viscosity,  heat  conductivity,  etc.  Examples  are 

the  Euler  and  Navier-Stokes  equations^  ‘«i  altogether  different  approach 

involves  the  use  of  more  fundamental  and  general  microscopic  formalisms. 

On  the  one  hand,  one  may  work  with  equations  describing  the  evolution 

of  one  particle  distribution  functions,  the  well-known  Boltzmann  equa- 
2 

tion  being  a  prime  example.  On  the  other  hand,  one  may  employ  equa¬ 
tions  relating  one-,  two-,  etc.  :*article  distribution  functions  such  as 

3 

the  hierarchy  of  equations  derived  from  the  Liouvi.le  equatio-,  which 
latter  treats  the  evolution  of  the  distribution  function  for  all  N 
particle:  in  the  system. 

The  fluid  equations  are  generally  considered  adequate  for 
treating  wide  classes  of  problems  in  gas  dynamics.  In  this  range  the 
microscopic  theory  would  not  yield  significantly  different  results.  In 
fact,  subject  to  certain  conditions  which  define  their  range  of  applica¬ 
bility,  the  fluid  equations,  together  with  explicit  formulas  for  the 

various  transport  coefficients  entering  into  them,  are  derivable  from 

2  4 

the  microscopic  theo-y.  * 

There  are,  however,  many  important  situations  in  which  the 
macroscopic  theory  does  not  give  a  correct  description.  In  general, 
this  occurs  when  the  length  or  time  scales  characterizing  phenomena  of 
interest  are  not  long  as  compared  with  the  scales  on  which  the  micro¬ 
scopic  quantities  fluctuate.  Examples  are  he  propagation  of  high- 
frequency  or  short  wavelength  waves,  and  behavior  near  boundaries.  In 
these  cases  one  must  properly  begin  with  the  microscopic  equations. 

The  solution  of  the  Boltzmann  equation  (or,  more  generally,  any 
equation  involving  only  the  single  particle  distribution  function)  is, 
in  general,  a  matter  of  considerable  difficulty  even  in  caseu 


corresponding  to  the  physically  simplest  situations.  Significant 
progress  has  been  confined  almost  entirely  to  the  study  of  two  limiting 
cases  in  which  two  different  approximation  procedures  can  be  applied. 

A  criterion  for  the  range  of  validity  of  the  approximate  methods  is 
provided  by  the  comparison  of  a  characteristic  time  T'  or  length  L  for 
the  relevant  process  with  the  average  time  f  or  mean  free  path  Lc 
between  particle  encounters. 

For  high  densities  ( T»  *c  or  L»  Lc)  the  Chapman-Enskog 
4  5 

theory  *  may  be  used.  The  first  approximation  of  the  theory  consists 
in  assuming  collisions  to  be  sufficiently  frequent  to  maintain  a  local 
thermodynamic  equilibrium.  The  next  approximation  corrects  the  distri¬ 
bution  function  by  terms  proportional  to  gradients  in  temperature  T  , 
flow  velocity  ^  ,  and  density  n  ;  this  corresponds  to  the  fluid  equa- 
tions  with  transport  coefficients  for  heat  conduction,  viscosity,  and 
diffusion.  This  high  density  region  is  in  fact  the  range  in  which  the 
fluid  equations  provide  an  adequate  description.  Higher  approximations 
of  the  Chapman-Enskog  theory  lead  to  correction  terms  proportional  to 
higher  derivatives  of  T}  £  )  n  .  The  successive  approximations  of  the 

4 

Chapman-Enskog  theory  correspond  to  an  expansion  of  the  distribution 
function  in  powers  of  the  mean  free  path  Lt  .  For  example  if  we  con¬ 
sider  sound  waves  with  wavelength  L»Lct  the  first  and  second  approx 
mations  are  already  sufficient  to  give  all  significant  features  of  the 
process.  When  L  becomes  comparable  to  Lc  however,  it  is  necessary  to 
go  to  the  tMrd  and  even  higher  approximations  to  obtain  adequate 
results;  the  third  approximation  already  involves  formidable  labor  onH 

3 

has  been  used  to  solve  only  the  simplest  problems.  Consideration  of 
higher  approximations  is,  in  any  case,  of  doubtful  value  since  the 
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entire  procedure  breaks  down  in  just  the  range  where  the  contributions 
from  these  higher-order  terms  becomes  important,  A  different  approach, 
using  expansions  in  terms  of  Hermite  polynomials  in  velocity  space,  has 
been  given  by  Grad.^  He  uses  some  moments  of  low  order  in  addition  to 
the  usual  ones  representing  n,  ^  and  T .  The  procedure  involves  a  gain 
in  simplicity  over  the  Chapman-Enskog  theory  but  is  still  quite  compli¬ 
cated.  In  any  event,  it  is  basically  a  high-density  theory. 

The  opposite  limiting  case  of  low  densities  (  r<<  \  or  L«  i^) 

has  been  studied  using  iterative  schemes  beginning  with  the  solution  of 

7  8 

the  "collisionless"  equation.  *  In  the  case  of  ionized  gases  immersed 
in  strong  electromagnetic  fields,  such  theories  have  been  used 
extensively.^  As  with  the  approximation  schemes  employed  in  the  high- 
density  case,  these  iterative  procedures  become  unwieldy  if  more  than 
one  iteration  is  necessary. 

It  would  naturally  be  very  desirable  to  have  a  method  capable 
of  treating  the  microscopic  equations  over  the  whole  range  from  low  to 
high  densities.  Unfortunately  the  describing  equations  are  generally 
non-linear,  and  even  when  linearized  are  extremely  intractable,  a  prime 
source  of  difficulty  being  the  term  representing  interparticle  colli¬ 
sions.  Relatively  little  work  has  been  done  in  this  intermediate  density 
region,  often  referred  to  as  the  kinetic  regime  Attention  has  generally 
focused  on  mathematical  properties  of  the  collision  operator,10  on  test 

5 

particle  treatments,  or  on  a  numerical  solution  of  the  kinetic  equa¬ 
tion  in  a  few  simplified  situations. 11  A  notable  exception  is  in  the 

12  13 

work  of  Chang  and  Uhlenbeck.  '  These  authors  treated  the  propagation 
of  small  amplitude  sound  waves  in  a  monatomic  gas  composed  of  i  ms 
interacting  via  an  inverse-fifth  power  force  law,  i.e.,  Maxwell 
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molecules.  They  were  able  to  show  that  the  linearized  Boltzmann  colli¬ 
sion  operator  for  this  case  generated  a  complete  set  of  velocity-space 
eigenfunctions.  Upon  expanding  the  perturbed  distribution  function  in 
terms  of  these  eigenfunctions,  they  obtained  from  the  linearized 
Boltzmann  transport  equation  an  infinite  set  of  coupled  algebraic 
equations  which  they  solved  by  successive  approximation.  Their  results 
were  in  quite  good  agreement  with  experimental  observations  on  collec¬ 
tions  of  neutral  atoms. 

In  view  of  the  difficulties  involved  in  solving  the  microscopic 
equations  in  the  kinetic  region,  considerable  interest  has  recently 
been  focused  on  the  mathematical  properties  of  the  terms  representing 
collision  effects.  Motivation  in  this  direction  has  been  based  in  part 
on  the  feeling  that  a  knowledge  of  the  spectral  properties  of  the 
(linearized)  collision  operator  would  lend  insight  into  the  kinetic 
behavior  of  the  system. Grad1^  has  considered  the  linearized 
Boltzmann  collision  operator  for  particles  interacting  via  the  genera, 
inverse  power  force  law, 

F  =■  */r  * 

where  K  is  a  constant  and  f  is  the  interparticle  separation.  Grad 
found  that  in  order  to  obtain  mathematical  results  it  was  necessary  to 
assume  the  interparticle  force  extended  over  a  finite  range;  i.e., 
angular  integrations  in  the  collision  integral  were  truncated  at  small 
deflections.  On  the  basis  of  this  assumption  was  able  to 
show  that  the  spectrum  consists  of  two  parts  ..  uiscrete  spectrum  and 
a  continuous  spectrum.  The  latter  is  bounded  away  from  zero  for  "hard" 
potentials  {  S>5  ),  and  approaches  zero  for  "soft"  potentials  (  5  <  S'  ). 
For  the  special  case  of  the  Maxwell  molecule  (S  »5")  there  is  only  a 
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discrete  spectrum.  For  5<  3  Grad  was  unable  to  find  the  spectrum. 
Recently  Fer2iger  has  used  Grad's  results  for  inverse  power- law 
molecules  to  show  that  the  linear  Boltzmann  collision  operator  generates 
a  complete  set  of  eigenfunctions.  The  form  of  the  eigenfunctions  was 
not,  however,  obtained.  In  any  event,  the  spectral  and  completeness 
properties  v/ere  not  obtained  for  the  Coulomb  potential  (5  =2). 

A  major  problem  in  plasma  physics  is  that  of  determining  the 
properties  of  an  isolated  hot  plasma;  any  material  probe  introduces 
impurities,  while  for  a  fully  ionized  plasma  the  emitted  radiation  is 
only  moderately  informative,  any  line  structure  arising  from  undesirable 
impurities. 

The  interaction  of  an  incident  beam  of  radiation  with  a  plasma 
has  proved  to  be  a  useful  method  for  determining  the  electron  density; 
since  the  plasma  acts  roughly  as  a  dielectric  with  coefficient 

<=  *  I"  )  60  ?e  -  Livv\z'/**) 

transmission  is  cut  off  below  the  plasma  frequency.  Radiation  above 
the  plasma  frequency  may  also  be  used  as  a  plasma  probe,  and  several 
experiments  have  used  the  modificiation  in  phase  velocity  produced  by 
the  dielectric  coefficient  £  as  a  measure  of  electron  density.  ^ 

It  has  been  known  for  some  time  that  the  scattering  of  photons 
or  material  particles  from  a  system  of  interacting  particles  yields 

1 7 

detailed  information  on  the  structure  of  the  scattering  system.' 

With  the  advem  f  intense  light  sources  such  as  Dulsed  ruby  lasers, 
considerable  attention  has  been  given  to  the  scattering  of  photons  by 
free  electrons  in  ionized  gases.  Several  authors  have  presented 
analyses  of  this  phenomenon,  usually  basing  their  descriptions  on 
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semi-intuitive  derivations,  and  employing  collisionless  kinetic 

1 8 

theories.  Notable  among  these  are  Salpeter,  and  Rosenbluth  and 

19  20 

Rostoker.  Shortly  thereafter  considerations  of  relativistic  effects 

21 

and  nonlinear  scattering  appeared. 

The  first  experimental  observation  of  photon  scattering  from 

22 

an  ionized  gas  was  reported  by  Bowles,  who  observed  the  scattering 
of  a  radar  beam  from  the  ionosnhere.  More  recently  many  workers  have 
reported  the  measurement  of  optical  photon  soectra  produced  by  scat¬ 
tering  from  ionized  gases  in  the  laboratory.  The  observations  are 
generally  in  remarkable  qualitative  and  quantitative  agreement  with 
theoretical  predictions,  which  is  a  rarity  in  plasma  physics.  An 

illustration  of  this  close  agreement  is  given  in  the  recently  reported 
23 

work  of  Anderson. 

Provided  the  photon  wavelength  A  is  of  the  order  of  the  Debye 
length  XD  or  larger,  the  scattered  spectrum  is  characterized  by  a 
narrow  central  peak  located  at  the  incident  frequency,  and  by  two 
symmetrically  placed  satellite  peaks  separated  from  the  central  peak 
by  A a> -scope.  ,  the  electron  plasma  frequency.  The  central  peak 
reflects  the  strong  coupling  of  the  electrons  to  the  ions  character¬ 
istic  of  long  wavelength  plasma  phenomena,  while  the  satellites  are 

attributed  to  the  resonant  scattering  of  photons  from  longitudinal 

18  19 

electron  plasma  oscillations.  ’  As  the  photon  wavelength  becomes 
large  in  comparison  with  the  Debye  length,  the  satellites  become 

19 

narrower  and  rapidly  increase  in  height.  This  has  been  attributed 
to  a  decrease  in  the  effect  of  Landau  damping  on  long  wavelength  plasma 
oscillations. 

24  25 

Recently  Ron,  Dawson,  and  Oberman  and  Fante  have  computed 
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the  scattered  photon  spectrum  for  a  fully  ionized  gas  including  the 

effect  of  collisions.  They  find  the  difference  between  their  results 

*  »  | 

and  the  collisionless  treatments  to  be  very  small,  of  order  A  ,  where 

3 

A.*r\^ois  generally  very  large;  in  fact  A  is  generally  assumed  large 

for  the  various  theoretical  models  to  be  valid. 

In  recent  years  considerable  effort  has  been  expended  on  the 

production  and  diagnosis  of  gases  that  are  only  partially  ionized. 

Examples  range  from  low  temperature  gas  discharges  to  relatively  high 

temperature  (  0e*)  and  high  density  plasmas  generated  by  the  laser 

bombardment  of  solids.  Since  photon  scattering  has  been  proven  to  be  a 

most  useful  tool  in  the  diagnostics  of  fully  ionized  gases,  it  is 

natural  to  expect  this  usefulness  may  be  extended  to  include  systems 

containing  significant  numbers  of  neutrals.  In  addition,  since  the 

Thomson  scattering  cross-section  for  electrons  is  several  orders  of 

magnitude  larger  than  the  ion  cross-section  or  the  Rayleigh  scattering 

26 

cross-section  for  neutrals,  we  would  expect  that  photon  scattering 
from  electrons  should  be  observable  even  when  neutral  densities  exceed 
electron  densities.  The  p  .ry  difference  between  fully-  and  partially- 
ionized  gases  in  this  respect  would  then  be  in  the  effects  of  neutral 
atoms  on  the  scattering  process,  since  charge-neutral  collision  frequen¬ 
cies  may  often  be  considerably  larger  than  their  Coulomb  counterparts- 

2.  Outline  of  this  Work 

The  purpose  of  this  dissertation  is  twofold.  In  the 
first  part  of  this  work  we  obtain  the  spectrum  and  prove  an  expansion 
theorem  for  the  linearized  ^okker-Planck  collision  operator  for 
particles  interacting  via  an  inverse-square  force  law.  In  the  second 
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part  of  this  work  we  employ  a  simple  collision  model  to  study  the 
effects  of  collisions  on  the  spectrum  of  photons  scattered  from  a 
partially  ionized  gas.  The  plan  of  this  paper  is  as  follows. 

In  Chapter  II  we  present  a  brief  survey  of  plasma  kinetic 
theory,  with  particular  emphasis  on  the  description  and  effects  of 
charged  particle  interactions.  Due  to  the  long  range  nature  of  the 
Coulomb  interaction,  when  an  ionized  gas  is  in  the  kinetic  regime,  a 
great  number  of  particles  are  "colliding"  simultaneously.  In  this  case 
the  simple  binary  collision  models  employed  to  treat  collections  of 
neutral  particles  are  often  deemed  to  be  inadequate,  As  a  consequence 
many  attempts  have  been  made  to  develop  suitable  kinetic  descriptions 
for  plasmas,  and  a  comparison  of  the  various  treatments  has  often  led 
to  confusion.  The  purpose  of  this  chapter  is  to  compare  a  few  of  the 
better  known  kinetic  models  in  an  attempt  to  cast  some  light  on  their 
similarities  and  differences. 

In  Part  1  of  Chapter  II  we  give  a  very  brief  phenomenological 
derivation  of  the  Boltzmann  collision  integral,  and  of  the  Fokker-PIanck 
collision  operator  for  inverse-square  law  forces.  Emphasis  is  on  the 
difference  between  short-  and  long-range  interactions. 

In  Part  2  of  Chapter  II  we  briefly  discuss  the  hierarchy  of 
equations  generated  by  the  Liouville  equation.  The  effect  of  corre¬ 
lations  between  particles  is  studied  by,  on  the  one  hand,  neglecting 
correlations  altogether  and,  on  the  other*  retaining  two  particle 

correlations.  In  the  first  case  we  obtain  the  Boltzmann  equation  for 

28 

short-range  interactions,  and  the  Vlasov  equation  for  long-range 

interactions.  In  the  second  case  we  obtain  a  kinetic  equation  developed 
29  30  31 

by  Lenard,  Guernsey,  and  Balescu.  With  proper  assumptions  this 
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equation  reduces  to  the  Fokker-Planck  equation. 

In  view  of  the  mathematical  difficulties  involved  i.i  the 

solution  of  the  various  kinetic  descriptions,  it  is  often  advantageous 

to  replace  the  more  accurate  and  less  manageable  collision  descriptions 

by  a  model  that  simplifies  the  solution  of  the  kinetic  equations.  In 

Part  3  of  Chapter  II  we  consider  such  a  simplified  collision  model, 

32 

generally  referred  to  as  the  Krook  model. 

The  ultimate  test,  of  any  theory  lies  in  a  comparison  of  the 
predictions  thereof  witn  experimental  observation.  Due  to  the  scarcity 
of  relevant  experimental  and  theoretical  information  in  plasma  physics, 

,t  is  often  instructive  to  compare  the  results  of  the  various  theories. 
One  hopes,  in  so  doing,  to  acquire  physical  insight  into  both  the 
structure  of  the  theories  and  the  as  yet  unobserved  properties  of  nature. 
In  the  last  three  parts  of  Chapter  II  we  make  such  comparisons  for  a 
few  illustrative  cases. 

In  Part  4  we  review  Some  recent  numerical  treatments  of  relaxa¬ 
tion  to  equilibrium  in  velocity  space.  In  some  cases  it  is  possible 
to  compare  the  Fokker-Planck,  Lenard-Guernsey-Balescu,  and  Krook  descrip¬ 
tions.  From  the  information  currently  available  we  observe  a  negligible 
difference  between  the  predictions  of  the  first  two  treatments.  The 
limitations  of  a  single  parameter  Krook-type  model  are  discussed  and 
compared  with  other  results. 

In  Part  5  of  Chapter  II  we  discuss  the  phenomenon  of  electron 
runaway  and  the  contribution  of  collisions  to  plasma  transport  para¬ 
meters.  The  failure  of  the  Krook  model  in  the  description  of  runaway 
is  noted,  and  the  results  of  a  computation  using  the  Fokker-Planck 
description  are  discussed.  The  majority  of  this  section  is  devoted  to 


a  comparison  of  the  different  collision  descriptions  in  the  computation 
of  the  plas.iia  electrical  conductivity.  For  field  frequencies  small 
compared  to  the  collision  frequency,  the  Krook,  Boltzmann,  and  Fokker- 
Planck  results  are  identical.  For  frequencies  above  the  collision 
frequency  but  below  the  electron  plasma  frequency  co^  ,  the  Boltzmann 
and  Fokker-Planck  results  differ  only  slightly.  For  frequencies  above 
tOpe  these  collision  descriptions  breax  down;  the  reasons  for  this 
failure  are  discussed.  For  high  frequencies  conductivity  computations 
baf ed  on  the  Vlasov  equation  and  the  first  two  members  of  the  BBGKY 
hierarchy  give  similar  results,  and  match  the  results  of  the  collision 
description  just  below  . 

In  the  final  section  of  Chapter  II  we  present  a  brief  review 

of  some  recent  work  on  collisional  effects  in  plasma  collective 

behavior.  Since  relatively  little  work  has  been  done  in  this  area, 

only  a  few  comparisions  of  the  different  collision  descriptions  are 

possible.  In  general  it  is  found  that  for  wavelengths  long  compared 

to  the  Debye  length  damping  is  primarily  collisional;  for  wavelengths 

33 

cf  the  order  of  the  Debye  length  or  less  Landau  damping  predominates. 
Moreover,  if  a  plasma  is  inherently  stable,  collisions  increase  the 
damping  of  small  amolitude  oscillations.  In  contrast  if  a  plasma  is 
unstable,  collisions  may  increase  the  growth  rate  of  the  i ns tabi lit”. 

In  Chapter  III  we  obtain  the  spectral  properties,  and  deve.^p  a 
expansion  theorem  for  the  linear  Fokker-Planck  collision  operator  for  a 
two-component  fully  ionized  gas.  In  Part  1  of  this  chapter  we  take 
advantage  of  the  small  electron-ion  mass  ratio  to  decouple  the  equa¬ 
tions  for  each  species.  For  convenience  we  then  concentrate  our 
attention  on  the  ion  collision  operator  and  later  discuss  the  extension 
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of  our  results  to  the  electron  case. 

To  develop  the  expansion  formula  we  follow  the  standard  method 
of  assuming  solutions  to  the  kinetic  equation  of  the  form 

rfct)  -  $xfc)  e*f  (- M)  . 

This  reduces  the  equation  to  the  form 

LU)  g/^  -  -  A 

where  L  is  a  three-dimensional  integrodifferential  operator.  In 
Part  1  we  show  that  ReX>o  and  as  we  would  expect  physically. 

In  Part  2  we  introduce  a  spherical  harmonic  expansion  which  replaces 
the  three  dimensional  equation  by  an  infinite  set  of  uncoupled  equations, 

where  is  a  singular  integrodifferential  operator.  These  are 

then  cast  into  a  self-adjoint  form  in  Part  3  by  introducing  a  suitable 
algebraic  transformation  on  the  functions  ^  •*/*,)  .  With 

ocundary  conditions  obtained  by  combining  the  original  kinetic  equation 
with  the  conservation  laws,  we  proceed  to  find  the  spectrum 
which  is  continuous  and  for  /-°,  I  consists  of  all  A > o £  *•  c/ 
and  for  X-Z  consists  of  all 

34 

Although  the  spectral  resolution  theorem  implies  the  existence 
of  an  expansion  theorem  for  self-adjoint  operators,  there  always  remains 
the  task  of  constructing  the  expansion  explicitly.  We  turn  to  this  task 
in  Part  4.  Since  L ^  is  singular  at  ^so  and  vr-* o©  ,  we  temporarily 
replace  the  interval  o  5  vr-o©  by  the  interval  o<ir1^irsvl<oeanrj  use 
a  result  of  Tamarkin  to  show  that  generates  a  complete  ortho- 

normal  set  on  this  interval.  To  return  to  the  original  interval  and 
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thus  obtain  the  desired  expansion  theorem,  we  use  the  above  completeness 
property  together  with  an  extension  of  the  theory  of  singular  differen- 

07 

tial  equations.  *  This  finally  yields  a  continuous,  normalizable 

set  of  functions  2^#  (V,  A*./*)  that  is  complete  with  respect  to 
functions  square  integrable  in  velocity  space.  Since  the  set  is  con¬ 
tinuous,  the  expansion  has  the  form  of  a  generalized  Fourier  integral. 

For  it  » (lG/y*)lf\e  have  found  asymptotic  forms  of  the  expansion 
functions.  Finally,  in  Parts  5  and  6  the  extension  to  the  electron 
kinetic  equation  is  discussed,  as  well  as  certain  implications  of  our 
results. 

In  Chapter  IV  of  this  dissertation  we  develop  a  theory  of 
photon  scattering  from  partially  ionized  gases.  Tho  starting  point  for 

38 

the  present  discussion  is  the  description  developed  recently  by  Osborn, 

which  treats  photon  scattering  from  a  fully  ionized  gas.  The  primary 

38 

concern  in  this  work  was  with  the  establishment  of  a  relationship 
between  the  observed  distribution  of  scattered  photons  and  the  dynamical 
and  statistical  characteristics  of  the  scattering  plasma.  It  was 
assumed  that  (i)  the  dynamical  variables  of  the  plasma  obeyed  the 
classical  equations  of  motion,  and  that  (ii)  the  plasmas  in  question 
were  sufficiently  highly  ionized  that  the  presence  of  neutrals  could  be 
neglected. 

In  Chapter  IV  we  examine  the  second  assumption  described  above 

Our  motivation  in  this  direction  is  twofold.  First,  it  is  not  clear 

a  priori  when  the  presence  of  neutrals  will  be  truly  negligible. 

Further,  when  their  effect  is  significant,  it  must  then  be  incorporated 

39  18 

in  a  description  of  the  scattering  process.  Vineyard,  Salpeter, 

40 

and  Feyer  have  discussed  the  contribution  of  neutrals  to  the  plasma 
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scattering  function,  but  each  based  his  discussion  on  semi -intuitive 
arguments  and  in  no  case  was  a  quantitative  description  developed. 

It  was  shown  in  reference  38  that,  neglecting  relativistic 
effects,  the  cross-section  presented  by  an  electron  to  a  photon  having 
frequency^  and  direction  for  the  scattering  of  that  photon  into 
( to 1  co  '  +  d  f  )  and  (  Q  \  -D- '  +■  d  A- '  )  is  given  by 


where  (Tj-  is  the  Thomson  cross-section.  The  so-called  scattering 
function  aJ)  is  given  by 


s- 


2  irh){ 


where 


s" 

_  o o 


tz 


-  itA^ 


6“  =  <  ■ 

The  normalization  is  such  that  /Ve  is  the  total  number  of  electrons  in 
the  scattering  volume,  and  the  function  is  the  Fourier-transformed 
configuration  space  electron  density  operator,  with  transform  variable 


H  . 

In  Part  1  of  Chapter  1/  we  briefly  review  some  of  the  more 

38 

pertinent  aspects  of  reference  38.  Following  Osborn,  we  employ  a 
classical  representation  of  the  density  operators  and  develop  and  dis¬ 
cuss  various  aspects  of  the  classical  scattering  function  Certain 
properties  of  a  thermal  equilibrium  plasma  relevant  to  the  computation 
of  the  scattering  function  are  also  discussed. 

In  general  neutrals  may  influence  the  scattered  photon  spectrum 
in  two  ways.  The  first,  and  perhaps  more  obvious,  contribution  is 
significant  when  the  number  of  photons  scattered  by  neutrals  into  the 
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frequency  range  of  interest  is  not  small  compared  with  the  number 
scattered  by  electrons.  Secondly,  neutrals  may  bias  the  scattered 
spectrum  through  their  influence  on  the  electron  density  operator.  We 
investigate  this  latter  effect  in  Part  2  of  Chapter  IV,  We  approach 
the  problem  along  two  somewhat  different  paths.  The  first  is  a  simple 
extension  of  the  treatment  in  ref.  38,  to  include  a  neutral  species  in 
the  description.  We  show  that  this  approach  does  not  lead  to  a  noti cable 
contribution  from  neutrals  unless  exceptionally  high  neutral  densities 
are  present.  At  this  point  we  take  a  somewhat  different  tack,  employing 
the  simple  Krook  model  to  represent  interactions  with  neutrals,  and 
leading  eventually  to  a  modified  scattering  function  displaying  signifi¬ 
cant  contributions  from  neutrals  at  neutral  densities  at  least  five  or 
six  oruers  of  magnitude  lower  than  in  our  earlier  treatment. 

Part  3  of  Chapter  IV  is  devoted  lo  a  presentation  and  discus¬ 
sion  of  the  classical  scattering  function,  with  a-..d  without  the  effects 
of  neutrals,  for  several  experimental  configurations.  Our  results  are 

oq  I o 

in  agreement  with  earlier  predictions  ’  ’  in  that  the  effect  of 

neutrals  is  to  heighten  and  narrow  the  central  spectral  peak,  while 
lowering  and  broadening  the  so-called  electron  plasma  w\.gs. 

In  Part  4  we  discuss  the  implications  and  limitations  of  the 
approximations  and  assumptions  employed  in  the  development  of  our 
results:  a  few  recent  experiments  are  cited  as  illustrations.  Exten¬ 
sions  of  the  present  theory  are  suggested  where  such  modifications  may 
lead  to  significant  differences. 


II.  A  BRIEF  SURVEY  OF  PLASMA  KINETIC  THEORY 
1.  Phenomenological  Kinetic  Equations 


One  of  the  earliest  and  still  most  successful  descriptions  of 

1  2 

a  collection  of  free  particles  is  the  famous  Boltzmann  equation  * 

which  expresses  how  FA  changes  in  time  due  to  streaming  and  to 
encounters  with  other  particles.  Here  F  *  A,  is  the  expected 
numbe.  of  type  A  particles  in  c/V|dV,  about  f,  jfj  at  time  t,  and 
'$♦  is  the  force  on  a  particle  of  type  A  and  repre- 

sents  the  time  rate  of  change  of  f  produced  by  inter-particle 
encounters.  For  our  purposes  we  assume  that  only  elastic  encounters 
are  important. 

1  2 

In  the  classical  theory  of  non-u. ti form  neutral  gases  * 
(ifVst),  is  taken  to  be  the  net  number  of  particles  of  type  A 
entering  the  phase  space  volume  element  c/3F,  c i*\rf  per  unit  time  due 
to  instantaneous  binary  encounters.  Thus  if  jj  is  the  relative 

velocity  of  a  colliding  pair,  then  the  flux  of  particles  of  type  8 
having  velocity  yv  incident  on  any  particle  of  type  A  having  velocity 


is 

—  i 


d\  P  -t ,t)  I  Vf,  — 


If  <J>6  (<^)  &)  is  the  differential  scattering  cross-section,  then  the 

number  of  particles  of  type  A  scattered  out  of  the  phase  space  volume 
^  3  • 

element  d  f,  d  vr,  into  da  during  dt  by  collisions  with  type  B 
particles  is 
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A,JV.  ^ 8^, , r, , t)  <j, <rAfe  (],e)Jajt  , 

The  number  scattered  into  this  same  volume  in  dt  is 


As  a  consequence  of  momentum  and  enrrgy  conservation  in  the  collision 
2 

we  have 


=  JVJVI 

and  hence 

(dV^  Jo.irw  (^0)  3  [P  Yft'lFVO  -  FVr.)  FV^Vt 

the  well-knwon  Boltzmann  collision  integral, 

A  binary  encounter  is  often  described  in  terms  of  the  impact 
parameter  b  .  the  distance  of  closest  approach  if  no  interaction  is 
present.  This  latter  quantity  is  related  to  the  differential  scattering 

3 

cross-section  tnrough  the  relation 


where  J-Q-  = 

The  derivation  of  the  Boltzmann  collision  integral  rests  on 
three  basic  assumptions,  and  as  yet  there  is  considerable  debate  as  to 
their  significance: 

(i)  The  collection  of  particles  is  sufficiently  rarefied  so 
that  only  binary  collisions  need  be  considered; 

(ii)  The  probability  of  finding  two  particles  in  J  V 
about  r  simultaneously  is  proportional  to  the  prod  ct  of 


their  individual  distribution  functions; 

(iii)  the  force  is  distinct  from  the  interparticle 
forces,  and  it  does  not  affect  the  collision  process. 

Assumption  (i)  implies  that  the  interparticle  potential  is  so  localized 
that,  for  a  sufficiently  dilute  gas,  the  chances  of  finding  more  than 
two  particles  within  "range"  of  each  other  simultaneously  is  negligible. 

While  this  assumption  seems  reasonable  for  dilute  gases  composed 
of  neutral  atoms  or  molecules,  it  causes  considerable  difficulty  in  the 
treatment  of  charged  particles.  The  difficulty  arises  because  electro¬ 
static  forces,  being  proportional  only  to  the  inverse  square  of  the 
distance,  permit  many  particles  to  be  within  range  of  each  other  at  a 

given  time.  To  see  this,  note  that  the  effective  interparticle  potential 

4 

for  a  fully  ionized  gas  not  far  from  thermal  equilibrium  is 

V  e 

\  4 

where  the  Debye  screening  length  is  given  by 


where  0  is  the  temperature  in  energy  units.  It  follows  that  the  number 
of  particles  interacting  simultaneously  is  roughly 

Y\  A  p  =  A 

or  (for  a  singly  ionized  gas) 

/i.  ~  Jje/8nel  )k~,/\] 

4 

which  is  usually  quite  large.  In  contrast,  charged  particles  will 

3 

suffer  large  angle  deflections  only  when  the  impact  parameter  is  of  the 
order  of  eye  ,  the  distance  of  closest  approach  in  a  head-on  collision. 
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Since  interparticle  spacings  are  ^  the  fraction  of  particles 

making  such  encounters  at  any  instant  is 


Thus  the  binary  collision  assumption  seems  inadequate  since  it  does 

not  appear  to  account  for  the  many  overlapping  long-range  encounters. 

5  6 

This  was  recognized  by  Chapman  ’  in  the  calculation  of  transport 
coefficients  for  an  ionized  gas.  Chapman  foetid  that  integrals  over  the 
impact  parameter  diverged  for  b-^eo  due  to  the  long  range  nature  of  the 
Coulomb  force.  He  overcame  this  difficulty  by  cutting  off  the  integra¬ 
tion  at  b  ~  ,  and  assuming  that  the  ♦'esultant  force  on  a 

particle  due  to  more  distant  encounters  could  be  represented  by  an 
internal  electrostatic  force  in  the  streaming  term  of  the  Boltzmann 
equation. 

Cohen,  Spitzer,  and  Routly^  sought  to  overcome  Chapman's  diffi- 

8  9 

culties  by  adopting  a  treatment  developed  by  Chandrasekhar  ’  in  a 

study  of  stellar  dynamics.  Chandrasekhar's  work  was  based  on  Jeans'^ 

demonstration  that  when  particles  interact  through  inverse-square 

forces,  the  cumulative  effect  of  the  weak  deflections  resulting  from 

the  relatively  distant  encounters  is  more  important  than  the  effect  of 

occasional  large  deflections.  Chandrasekhar  noted  the  strong  similarity 

between  the  Brownian  motion  of  a  colloidal  particle  and  the  motions  of 

particles  interacting  via  inverse-square  forces;  his  treatment  is  based 

11  12 

on  a  description  of  Browman  motion  due  to  A.  D.  Fokker  and  M.  Planck.. 

While  Cohen,  et  al/  restricted  themselves  to  slightly  anisotropic  velocity 

distributions,  the  general  case  of  arbitrary  distribution  functions  was 

1 0 

considered  by  Rosenbluth,  MacDonald,  and  Judd.  J  We  present  a  brief 
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9  1 3 

derivation,  following  Chandrasekhar  and  Rosenbluth,  et  al , 

Assume  that  there  exist  time  intervals  At  long  enough  for  a 

particle  to  suffer  a  large  number  of  weak  deflections  but  short  enough 

for  the  net  mean  square  change  in  velocity,  </Ajr|1‘>  ,  to  be  small 

compared  with  the  mean  square  velocity.  Let  PA{k}A  if)  denote  the 

probability  that  in  the  time  interval  At  a  particle  of  type  A  having 

velocity  jr  undergo"'  a  displacement  As  .  Assuming  that  P*  dC3S  not 

depend  explicitly  on  time,  the  distribution  function  for  the  species 
g 

is  then  given  by 


FVl.St-)  '  PVr-4iVr).  (2.D 


Since  At  and  As  art  ~oth  assumed  small,  the  integrand  is  expanded  in 
a  Taylor's  series: 


Using  the  fact  that 


(jV&z)  P  ViT;*/}  -  I  , 

equation  (2.2)  gives  for  the  time  rate  of  change  o*  f-  "resuHing 
from  the  cumulative  effect  of  small  deflections 
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where 

\a£]as  PV/; 'wVr, 

=  (JS^1  pV*,1 *£  (2,4) 

and  so  forth. 

Since  jj/i «/■)  represents  the  mean  change  in^j/-  resulting 
from  encounters  duriny  dtr,  we  write 

^1A=  2  jjV'  FVrOfd/i-^g^,©)^  tt 

6  ' 

s  <aOA  4t  (2.5) 

and 

■  I  Ov,FVr)0-tt«;,^lA')3^4JrAi,f 
4  A  8 

=  <ar*r>A2it)  (2’6) 

etc,  for  higher  order  terms,  where  ^  s  l£~£'l  is  the  magnitude  of 

1 3 

the  relative  velocity.  Rosenbluth,  et  al_  evaluated  the  Fokker-Planck 
coefficients  (2.5)  and  (2.6)  assuming  a  two-particle  Coulomb  cross- 
section.  Since  the  integrals  over  JX  diverge  logarithmically  it  large 
impact  parameters,  they  followed  Cohen  et  aj7  and  cut  off  the  integra¬ 
tion  at  b^x  -v  Ap  ,  the  Debye  screening  length.  After  some  alge¬ 
braic  manipulation  Rosenbluth  et  al  obtain 


1^1  jl  lT 

^  eft /A  0 


where 
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and  is  the  reduced  mass.  Assuming  the  logarithm  to  be  a  suf¬ 

ficiently  slowly  varying  function  of  the  relative  velocity,  the  approxi¬ 
mation  is  made 

%  et 


r  ^ 

1  A&  - 


•aJ£  J  \e_ 


(2.8) 


v* 
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(It  can  be  shown  that  this  approximation  introduces  an  error  of  less 
fhan  one  percent  in  the  determination  of  transport  coefficients  based 
on  (2.7).) 

In  the  derivation  of  (2.7)  only  those  terms  porportional  to 
Ia\  AAfthave  been  retained  from  the  expansion  (2.3)  and  in  the  evalua¬ 
tion  of  the  Fokker-Planck  coefficients  (2.5)  and  (2.'').  All  other  terms 
can  be  shown13  to  be  down  by  a  factor  (A*  Thus  the  Fokker- 

Planck  equation  may  be  viewed  as  an  expansion  in  powers  of  the  ratio  of 
mean  kinetic  energy  to  potential  energy  at  a  separation  of  X  0  »  and  is 
sound  for  sufficiently  large,  or  roughly 

O^/jur)  >  |0“H  * 

Equation  (2.7)  may  be  transformed  into  a  more  symmetric  form. 

Since 

l  JdV  FVr  V  =  -  t  (dV' 

_  i  (  jj  i  ;re/V')  A 

-  i)d  *  Tr7-  •  >r>r  J 

(2.7)  can  be  rewritten  in  the  form 


(&)  = 
1  yt'c'ii 


I  i  r..  ^jJv{FVr)f  -  i  g, B 

8  ~ 
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This  fo.'m  or  the  Fokker-Planck  equation  was  published  by  Landau  in 
1936. 


While  the  above  analysis  treats  the  effects  of  a  large  number 

of  overlapping  small-angle  deflections  occurrinq  in  the  time  interval 

tt  ,  the  form  of  the  probability  employed  in  (2.5)  and  (2.6)  still 

assumes  that  these  small  angle  deflections  are  themselves  due  to  binary 

encounters.  Moreover,  we  have  retained  the  assumption  (ii)  of  the 

joint  probability  being  proportional  to  the  product  of  the  two  singlet 

probabilities  and  have  thereby  ignored  possible  correlation  effects. 

In  addition,  the  present  description  explicitly  ignores  the  effects  of 

the  relatively  infrequent  large  angle  deflections.  Both  Cohen'7  and 
13 

Rosenbluth  and  their  co-workers  suggested  the  inclusion  of  a  Boltzmann 
collision  operator  for  impact  parameters  below  an  unspecified  critical 
value.  At  the  same  time,  they  suggested  the  effect  would  usually  be 
negligible.  Finally,  we  note  in  passing  that  the  Fokker-Planck  equation 

4 

as  displayed  above  may  be  obtained  by  Taylor  expandina  the  integrand 
of  the  Boltzmann  collision  integral  v.ith  Coulomb  cross-section,  cutting 
off  the  integrals  at  -  A  p  ,  and  retaining  only  dominant  terms. 

Several  authors  have  presented  descriptions  which  take  into 
account  the  electrostatic  properties  of  the  plasma. 16,17,18  yhe 
is  to  consider  a  test  particle  as  being  subject  to  local  fluctuating 
electric  fields,  and  then  calculating  the  Fokker-r lanck  coefficients  on 


this  basis.  In  all  cases  the  results  are  quite  similar,  and  are  tanta¬ 
mount  to  incl  .ding  a  dielectric  constant  in  the  functions  in  (2.7). 

An  advantage  of  this  work  is  that  the  Debye  length  enters  the  descrip¬ 


tion  as  a  natural  cutoff  distance,  ,‘ithout  having  to  be  introduced  in 


an  ad  hoc  fashion  as  before.  In  contrast,  the  treatment  yields  a 
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divergence  at  small  impact  parameters  due  to  an  improper  handling  of 

close  encounters.  To  overcome  this  a  cutoff  is  postulated  at  an  impact 

parameter  of  the  order  of  the  distance  of  closest  approach*^  e  y<9  . 

19 

Hubbard  surmounted  this  latter  difficulty  by  retaining  the  entire 
infinite  series  in  the  expansion  (2.2). 


2.  Kintetic  Equations  Derived  From  the  Li ou vi lie  Equation 


The  construction  of  kinetic  theories  on  the  basis  of  phenomeno¬ 
logical  considerations  naturally  raises  questions  concerning  the 
validity  and  range  of  applicability  of  the  various  descriptions.  Per¬ 
haps  the  most  satisfactory  scheme  for  surmounting  these  difficulties  is 
to  begin  with  the  most  complete  (and  intractable)  description  available, 
the  Liouyille  equation. 

For  simplicity  we  consider  a  collection  of  /V  indistinguishable 
particles,  occupying  a  volume  V  ,  with  no  external  forces.  The 
generalization  to  a  multicomponent  system  is  straightforward  but  tedious, 
and  will  be  indicated  by  reference  where  appropriate. 

We  define  ( t)  )  C  r  /,  N  as  the  N -particle 

distribution  function  such  that 

F„  F  J  W 

i  =  i 

is  the  probability  of  finding  the  i ^  particle  in  V}-  for  each 

of  the  N  particles,  all  at  time  t  .  Clearly 


(.If  JV.d V,- F„  =  I. 

/  L  *■  / 


on 

According  to  Liouville's  theorem,  is  governed  by  the  equation 


I  f 


1 5, 

tr. 


I  ■= 


y  ■  \in 

fe >*• 


n 


t 


~  o 


(2,10) 
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where 


f.  =  r. 

' 1  ) 


*  =  -  i  h  Y&), 


Ykt)  -  Z  v"'Yu.'-Xji) 

l  i 

and  is  the  two-body  potential.  Thus  we  have 

\  +  -i  ^r.  ** Z_  L- 

dl  L~\  f* i 


*r(  ‘  v.  =c 


(2.11) 


The  information  contained  in  the  Liouville  equation  (2.10)  or 
(2.11)  is  all  inclusive  but  in  general  is  inaccessible  because  of  its 
complexity.  A  method  for  extracting  useful  information  has  been  devel¬ 
oped  independently  by  Bogolyubov,  Born,  Green,  Kirkwood,  and  Yvon,  and 

is  summarized  ir  Montgomery  and  Tidman  or  in  any  good  text  on  statis- 

21 

tical  mechanics  such  as  de  Boer  and  Uhlenbeck, 

To  develop  the  so-called  BB(/KY  hierarchy,  we  define  reduced 
probabiliiy  distributions 

5»VS(/  ,S<N.  (2.12) 

/  i*  S-H 


Multiplying  (2.11)  by  and  integrating  as  in  (2.12)  we  obtain  an 
expression  for  Fs  which  involves  ■  In  particular,  for  s  -i  and 
S  =  2  ue  have 
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While  the  Liou  ille  equation  uniquely  determines  the  evolution  of  PN 
from  the  initial  condition  ^(t-o)  t  the  derived  relations  for  F^ 
do  not  uniquely  determine  the  evolution  of  F^  from  Fj  fc'0').  Instead, 
the  interactions  of  our  S  particles  with  the  remaining  M-s  particles 
are  summarized  in  terms  of  a  typical  S-H  particle  and  we  require  a 
knowledge  of  to  solve  the  problem. 

We  shall  now  see  how,  for  a  given  set  of  circumstances,  the 
relations  for  and  can  be  simplified.  We  will  make  assumptions 
consistent  with  the  ohysical  conditions  which  we  desire  to  treat  and  on 
the  basis  of  these  assumptions  introduce  approximations  which  are 
designed  to  retain  the  pertinent  elements  of  information. 

Our  first  approximation  is  of  quite  general  validity.  We  assume 
N  to  be  a  large  number  and  expand  in  terms  of  N'1  ,  At  the  same 
time,  we  allow  the  volume  of  the  system  to  become  arbitrarily  large, 
but  such  that  hi/V-tn  remains  finite.  This  removes  interactions  with 
the  boundary,  and  reduces  (2.13)  and  (2.14)  to 


A.  The  Absence  of  Correlations.  The  Boltzmann  and  Vlasov  Equations. 
Writing  r,)XV)t)  in  the  form 


Fz 


V, 

I  v'  1 ) 


Oj  -o  ,  the  s,  tem 

(2.15),  (2.16)  reduces  to  an  equation  for  F,  , 


where  ^  is  a  correlation  function,  and  taking 


u 


t  vr 

^  i 


>£. 


w*-  vFjo 


F,A).  (2.17) 


It  is  interesting  to  consider  (2,17)  in  two  contrasting  cases.  In  the 
first,  we  consider  a  dilute  gas  with  short-range  interparticle  poten¬ 
tials  V,z  ,  such  that  the  effective  range  of  the  potential  is  much 
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smaller  than  the  mean  interparticle  separation.  Following  Schonberg 
23 

and  Osborn  ,  we  rewrite  (2.17)  in  the  form 


which  is  equivalent  to  (2.17)  since  the  additional  terms  are  zero  as  may 

23 

be  verified  by  partial  integration.  Following  Osborn  we  now  introduce 
a  set  of  assumptions  designed  to  lend  plausibility  to  the  succeeding 
argument. 


(i)  The  potential  V^is  appreciable  only  over  a  suf¬ 
ficiently  well  defined  region  of  radius  such  that  the 

integration  in  (2  18)  is  effectively  confined  to  the 
volume  rx  I  <  cl. 

(ii)  The  length,  (k  ,  is  small  compared  to  the  mean  inter¬ 
particle  spacing,  i.e,,  <*.  r,/3«  {  ,  This  implies  that 

<K.  can  be  chosen  such  that  the  probability  of  finding  more 
than  two  particles  within  i  distance  a.  of  each  other  at 
any  time  is  negligible,  i.e.,  we  choose  so  that  the 
binary  collision  assumption  is  in  some  sense  justified. 

(iii)  The  above  choice  of  a,  leads  to  a  sufficiently 
restricted  region  for  the  spatial  integration  in  (2.18)  so 
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that  P,  F,  F,  {snjrtl  t) 

for  all  tx  defined  by 

(iv)  The  product  Fffi)F,fe)  does  not,  vary  appreciably 
over  time  intervals  of  order 

r  =  a,  //  if,  -  I  . 

(v)  Two  particles  within  a  distance  a.  of  each  other 
will  be  presumed  to  be  interacting  so  strongly  that  they 
may  be  regarded  as  effectively  decoupled  from  their 
envi ronment. 


Assumption  (v)  permits  us  to  interpret 

i  3  V'1 


II  ***  I 


in 


-L  }¥-  -  -  a  *  -  a 

^  c.  -xi  - r 


(2.19) 


where  ftj  and  av  are  the  accelerations  experienced  by  particles  1  and 
2  respectively  throughout  the  duratiun  of  their  close  encounter.  The 
approximation  (2.19)  plus  assumptions  (iii)  and  (iv)  then  enables  us 
to  approximate  che  integrand  in  (2.18)  as  follows: 


*  i  tt  •  f;  U, ,  sr(|  t)  f;  (x, ,  t)  * 

-  -  S/&,  »*)(;£.>  ft,*)  -  F,fo1y1,OF/x„!(\,t) 

Evidently  we  may  interpret  the  velocities  (  jr*  )  as  the 

precollision  velocities  of  a  pair  of  particles  entering  into  a  strong 
binary  interaction,  whereas  (  i£l;^fv)  are  the  post-collision  velocities 
of  the  same  pair.  Introducing  the  notation 

(V  -  *  y,  \r  -  a  r)  *  At,  v' ) 

\  —  f  ~  I  }  v  ^  x  L  /  l  }  —  \*  / 
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we  see  that  we  nay  now  approximate  (2.18)  as 


(c/5 


hr.  -  ^ 


1  OFlCr,)ri)t)'fffcnr,t)FJS»ty)l(2  20) 

* 


Note  now  that  nothing  in  the  integrand  of  (2.20)  depends  on  £v  ,  Thus 
the  space  integration  is  readily  performed  if  we  introduce  the  variable 
change  R  =  r,  -rv  whence 


pH  =  j  .  (2,21) 

I  l£l<A 

If  we  take  /?  to  be  the  distance  between  the  two  particles  during  their 
close  encounter  and  introduce  a  cylindrical  coordinate  system  ' see 
figure  (2.1)  )  with  1  -axis  parallel  to  the  pre-collision  relative 

velocity,  then  (2.21)  may  be  written  as 

It  *./t. 

j  (  Ul if  j  •  (2  22) 

<fi'o 

The  quantity  is  clearly  the  center  of  mass  differential  scat¬ 

tering  cross-section  introduced  earlier,  which  we  may  write  as  o ~  Jfl 
Thus  (2.21)  may  be  written  as 


A 


j(T  Jfl 
jn 


(2.23) 


,  =  CL 

ir.-rA* 

It  is  understood  that  the  restriction  on  the  range  of  impact  parameters 
as  displayed  in  (2.22)  implies  a  corresponding  limitation  on  the  range 
of  the  angular  integration  in  (2.23).  In  these  terms  then,  equation 
(2.18)  becomes,  for  the  dilute  short-range  approximation, 
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Figure  2.1 

Details  of  Binary  Interaction 
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Equation  (2.24)  is,  apart  from  normalization  differences,  the  Boltzmann 
equation  in  the  absence  of  external  forces. 

Finally,  let  us  consider  again  equation  (2.17)  when  V11  is  the 
Coulomb  potential.  In  this  case,  as  we  have  seen  earlier,  the  notion 
of  a  binary  interaction  length  4,  is  indistinct  when  the  number  of 
particles  in  a  Debye  sphere  is  large.  !*ie  could  pick  a  length  a.  say  of 
order  a  few  times  the  distance  of  closest  approach,  but  the  choice  is 
vague.  In  this  case  i:,«?  potential  term  would  be  divided  into  a  "binary" 
term  and  a  "collective"  term,  this  latter  corresponding  to  the  force 
term  in  the  Boltzmann  equation.  For  example  if  we  consider  only  collec¬ 
tive  effects,  then  we  may  regard  potential  term  in  (2.17)  as 
producing  an  internal  electric  field  which  is  a  function  of  r  alone 
lie  define  this  field  by 


1V11 

>1, 


and  write  (2.17)  in  the  form 


15 

u 


IL  -e_g 


(2.25) 


which  is  known  as  the  Vlasov  equation. 


B.  The  Inclusion  of  Correlations.  The  Lenard-Guemsey-Balescu  Equation. 

It  is  clear  from  the  foregoing  discussions  that  the  neqlect  of 
correlations  in  the  description  of  an  ionized  qas  is  a  simplification 
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which  is  difficult  to  justify,  and  which  leads  to  problems  such  as  the 

divergence  at  large  impact  parameters.  Since  the  8 6 &KY  hierarchy 

(2.15),  (2.16)...,  includes  correlations,  we  might  hope  that  the 

kinetic  equations  we  derive  from  it  will  be  intrinsically  free  from 

this  divergence,  In  the  following  we  present  the  salient  points  of  an 

important  kinetic  equation  derived  from  the  hierarchy  which  includes 

24 

two-particle  correlations .  We  outline  the  work  of  Lenard  and 

25  26 

Guernsey,  which  is  based  on  earlier  results  of  Bogolyubov. 

27 

We  begin  with  Mayer  cluster  expansions  of  the  distribution 
functions,  similar  to  that  used  in  the  previous  section.  We  have,  with 
(!■)  =  (tC,X{)  , 

=  FiAl  (2.26) 

and 

F3(V,3,t)=  * 

+  F-  h +•  F,  72,+')  j  0,s,i)  +■ 
f  t  h(i ,2,5,0  <2-27> 

where  h  is  a  three- particle  correlation  function.  We  work  in  ;he 
so-called  plasma  limit,  r\  Ao>>!  ,  and  we  further  assume  the  three- 
particle  correlations  to  be  negligible  (K*o).  We  consider  only  homo¬ 
geneous  plasmas,  still  with  no  external  forces,  so  that  the  various 
are  independent  of  position  and  the  correlation  functions 

depend  on  JN  and  £■  only  in  the  combination  r  =. 

26 

Bogolyubov  assumed  that  as  a  consequence  of  the  assumption 

a 

y\Xq  >^>  I  ,  the  correlation  functions  would  vary  much  more  rapidly 
in  time  than  any  of  the  distribution  functions  p ,/«')"£)  .  He  thus 
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suggested  that 

(1)  The  distribution  functions  Fj(R£)  may  be  considered 
time  independent  in  solving  for 

(ii)  The  asymptotic  value  of  the  correlation  function, 

,  represented  as  a  functional  of  the  F,  {i\t) >  may 
be  used  in  solving  for  these  latter  quantities. 

(iii)  As  a  consequence  of  (ii)  the  information  contained 
in  the  initial  condition  would  not  appear  in  the  Fj 

and  hence  we  may  take  ^  /V,  o  )  =•  O  . 

Employing  these  assumptions,  Bogolyubov  was  able  to  reduce  (2.15) 
and  (2.16)  to 

y  •  i  Jjv  ^ .  jjk  «.»> 

where  &  is  a  time  independent  function  of  P(:  R) ,  and  is  the  solution 
of  the  integral  equation 

<*(£,£,  iO  *  K(z,sr,  -rJ  •  F&.)  ~  F^()J  4‘ 

r>F^)  w,  .  \  -)F fa)  rrK  u-  ~)l 

'  TTe(^a>  (2  29, 


whose  kernel  is 


,  (  ~tJ.  WhlrZ tl) 
i  vr  ' 

Jo  **• 

The  normalization  is  jdVFsn  and  we  take  \/,l»  eVr  . 
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Lenard  and  Guernsey  reduced  the  system  of  equations  (2.28), 


(2.29)  to  a  single  equation  for  F  by  following  Bogolyubov's  suggestion 


,»_26 


-  34  - 


of  introducing  spatial  Fourier  transforms.  Their  result  was  also  ob- 

28 

tained  by  Balescu  using  a  diagram  technique; 


(2.30) 


where  Q  is  a  symmetric  second  rank  tensor  whose  components  are  given 
by 


)l 


(2.31) 


Here  fc  »  i/it  i.  and 

“  '^)-i(J3re-^V'VO=  »£ 


xi  -  \ 

Z'(\ 


The  integrand  in  (2.31),  aside  from  the  -function,  becomes 
independent  of  “fc  for  k  +  O  (large  impact  parameters),  so  that  the 
integral  converges  in  th's  limit.  Thus  the  incorporation  of  correla¬ 
tions  in  the  kinetic  description  removes,  in  a  natural  way,  the  trouble¬ 
some  divergence  at  large  impact  parameters.  For  ‘k  large,  the  integrand 

behaves  as  k~Z  and  hence  the  integral  diverges  for  k-*00  .  This  is 

17  18 

the  same  divergence  encountered  in  the  electrostatic  treatments  ’  of 
the  Fokker-Planck  equation  discussed  above,  and  is  due  to  the  neglect 
of  the  three-particle  correlation  h  when  two  of  the  three  are  of  +he 
order  of  €*’/ 9  apart.  In  this  case  h  is  the  same  order  of  magnitude 
as  the  binary  correlation  cj  ,  and  cannot  properly  be  ignored.  To 
achieve  convergence  in  (2.31)  it  is  therefore  necessary  to  make  a  short 
range  cutoff,  i.e.,  integrate  only  within  the  sphere 


s  K  *  e/e 1 . 


-  35  - 


24 

Lenard  simplified  (2.31)  somewhat  by  showing  that,  if  one 

neglects  contributions  to  Q ,<<  coming  from  speeds  greater  than  a  few 

times  (ZQ/\ Q  reduces  to 
)  ** 


(2.321 


The  approximation  introduced  in  going  from  (2.31)  to  (2.32)  requires 
that  f/y)  must  be  small,  for  velocities  greater  than  a  few  times  the 
thermal  speed,  as  compared  with  its  values  at  lower  speeds.  Equation 
(2.30)  with  (2.32)  is,  aside  from  a  missing  factor  of  2,  just  the 
Landau  form  (2  9)  of  the  Fokker-Planck  equation. 

We  might  consider  the  foregoinq  as  a  rather  convincing  argument 
supporting  the  use  of  the  Fokker-Planck  equation  for  ionized  gases  with 
distribution  functions  satisfying  the  above  condition,  and  which  satisfy 
the  Bogolyubov  assumptions  (i)  through  (iii).  These  assumptions  are 
violated,  for  example,  by  systems  in  which  F  and  ^  vary  on  the  same 
time  scale.  Examples  are  the  interaction  of  high-frequency  waves  with  a 
plasma,  or  rapidly  growing  instabilities. 

The  generalization  of  equation  (2.30)  for  a  multi compt  ant 

4 

system  may  be  found  in  Montgomery  and  Tidman.  The  generalization  to  a 

30 

constant,  uniform  magnetic  field  was  worked  out  by  Rostoker,  and  is 

considerably  more  complex  than  (2.30),  (2.31). 

Working  with  a  multiple  time  scale  theory  developed  by  Bogolyubov 

31 

for  certain  problems  in  nonlinear  mechanics,  Frieman  and  Book  have 
developed  a  kinetic  equation  for  homogeneous  field-free  systems  that  is 
free  of  divergences  for  all  impact  parameters.  For  small  impact  para¬ 
meters  their  result  resembles  the  Boltzmann  collision  integral;  elsewhere 
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it  is  similar  to  equation  (2.30) 

The  development  of  plasma  kinetic  equations  is  an  active  field, 
and  it  is  beyond  our  purpose  here  to  survey  the  topic  in  its  entirety. 

In  concluding  this  section  we  mention  only  briefly  some  of  the  problems 
of  current  "'nterest  in  this  area,  and  related  work. 

In  our  derivation  of  equation  (2.30)  we  assumed  a  homogeneous 
system  with  no  external  forces.  We  noted  earlier  that  the  phenomeno¬ 
logical  kinetic  equations  generally  assume  that  inhomogeneities  and 
local  force  fields  do  not  significantly  interfere  with  interparticle 
encounters.  While  these  may  oftrjn  be  valid  assumptions,  we  would  expect 
them  to  fail  when  length  and  time  scales  characteristic  of  inhomogene¬ 
ities  or  of  local  forces  ai\  small  compared  with  A  0  or  ^0/trr  (the  time 
required  for  the  establ ishmenc  of  a  Debye  screening  cloud  about  each 
particle'  where  \rr  is  the  thermal  speed  of  the  particles.  At  this 
writing  a  kinetic  equation  that  includes  these  generally  neglected 
effects  is  not  available.  Many  attempts  have  been  made  to  overcome 
these  difficulties.  Perhans  the  most  noteworthy  are  in  the  work  of  Bchm 
and  Pines,^  Rostoker  and  Rosenbluth,33  Dupree,^4  and  Frieman.^5  While 
none  of  these  authors  have  presented  a  treatment  satisfactory  to  ail, 
their  work  lends  considerable  insioht  into  the  problem  at  hand,  and 
often  presents  novel  and  significant  mathematical  tools.  This  work  is 
certainly  a  good  starting  point  for  the  interested  leader. 

3.  A  Simple  Collision  Model 

Me  have  so  far  been  concerned  with  the  development  of  a  kinetic 
equation  for  the  description  o *  ionized  gases.  While  a  satisfactoiy 
theory  is  yet  unavailable,  it  is  apparent  from  our  earlier  discussion 
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that  the  more  rig.  'ous  the  treatment,  the  less  tractable  it  if,  in  terms 
of  analytic  solution.  In  point  cf  fact,  this  problem  is  not  unique  to 
the  study  of  ionized  gases.  The  Boltzmann  transport  equation  has  so 
far  been  solved  only  for  a  few  special  cases,  namely  the  inverse-fifth- 
power  force  law  suggested  by  Maxwell,  and  for  certain  scattering  ker- 

Of: 

nels  relevant  to  neutron  transport  theory.  In  view  of  these  difficul¬ 
ties,  it  is  often  advantageous  to  replace  the  mere  accurate  and  less 
manageable  collision  descriptions  by  a  model  that  simplifies  the  solu¬ 
tion  of  the  kinetic  equations.  With  the  current  paucity  of  experimental 
information  on  fully  ionized  gases,  such- a  model  may  often  be  a  good 
starting  point  for  the  interpretation  of  the  little  information  that  is 
available.  We  thus  consider  briefly  a  collision  model  designed  to 
satisfy  the  conservation  laws  and  an  H-theorem,  and  which  considerably 
simplifies  the  mathematical  ana’ysis. 

One  of  the  earliest  of  these  models,  generally  referred  to  as 
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the  Krook  model,  was  developed  by  Bhatnagar,  Gross,  and  Krook.  A 

similar  but  somewhat  simpler  model  was  suggested  independently  by 
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Welander.  The  first  Krook  model,  for  a  single  component  system,  is 


($£)  =  0-  [-  Ffr,  *;*)  t  N(s,  t)  $  (s,  ; f,  t) )]  (2  r) 
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and 


3  efr.t) 
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where  ^  and  0  define  the  flow  velocity  and  kinetic  temperature  (in 

enerqy  units)  at  (t}t)  ,  and  o-  is  the  parameter  of  the  model. 

The  dimensions  of  N<r  are  inverse  time,  and  <r  is  generally 

chosen  to  yield  an  appropriate  cc.lision  frequency  on  the  basis  of 

39 

phenomenological  considerations.  For  an  ionized  gas,  <r  is  gener- 
39 

ally  chosen  such  that 


(T  - 


A  *hA 


While  the  Krook  model  is  highly  nonlinear,  it  is  considerably 
simpler  than,  e.g.,  the  Boltzmann  collision  operator,  since  the 
distribution  function  enters  in  (2.33)  in  a  simple  way:  as 

F)  ^  c/V  F  )  ^c/VjrF^  ^c/Vir1F, 


For  the  particular  case  of  small  amplitude  perturbations  near  equili¬ 
brium,  the  linearized  form  of  the  Krook  model  permits  solution  of  the 
kinetic  equation  in  closed  form  for  several  interesting  cases. ^7’40 
To  linearize  (2.33)  we  write 
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where  €  is  a  small  parameter,  and  where 
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We  thus  obtain,  neglecting  terms  in  €  , 

-  <r  N„  jj 


(2.34) 


The  extension  of  the  Krook  model  to  a  two-component  system  was 

40 

first  given  by  Gross  and  Krook,  together  with  an  application  of  the 

linearized  version  to  oscillations  in  a  fully  ionized  gas.  More 

41  42  43 

recently  Sirovich,  Liboff,  and  Oppenheim  have  presented  similar 
collision  models  for  a  general  multicomponent  system.  While  all  of  these 
models  are  similar  in  form,  they  are  not  identical. 
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The  relationship  between  the  Krook  model  and  the  more  sophisti¬ 
cated  collision  models  rests  primarily  on  intuitive  grounds.  In  one 

special  case,  however,  a  more  direct  relationship  has  been  demonstrated. 

44 

For  the  case  of  a  Maxwell  molecule,  Grose  and  Jackson  showed  that  the 
linearized  Boltzmann  collision  operator  yields  the  linearized  version 
of  (2.33)  if  all  the  non-zero  eigenvalues  of  the  collision  operator  are 
approximated  by  a  single  constant  which  appears  as  tr  in  (2.33). 

In  Figure  (2.2)  we  have  indicated  for  convenience  the  relations 
between  the  kinetic  theories  discussed  above. 


Fig.  2.2  Relations  Between  Kinetic  iheories  Discussed  in  Text 


4.  Relaxation  in  Velocity  Space 

The  ultimate  test  of  a  theory  lies  in  a  comparison  of  the  pre¬ 
dictions  thereof  with  experimental  observation.  At  the  present  writing, 
the  scarcity  of  relevant  experimental  and  theoretical  information  is  a 
major  source  of  difficulty  in  kinetic  physics.  This  ->s  particularly 

4 

so  in  the  physics  of  fully  ionized  gases  where  on  the  one  hand,  the 
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maintenance  of  a  plasma  and  reliable  observation  techniques  still  present 
many  unsolved  problems  and  on  the  other  hand,  tractable  theories  allow¬ 
ing  for  inhomogeneities,  rapid  temporal  variation,  and  boundary  effects, 
are  generally  unavailable. 

In  the  presence  of  these  difficulties  it  may  be  instructive  to 
compare  the  results  of  the  various  theories,  with  the  hope  that  this 
would  lend  physical  insight  into  both  the  structure  of  the  theories  and 
the  as  yet  unobserved  properties  of  nature.  In  this  and  the  following 
two  sections  we  intend  to  make  such  comparisons  for  a  few  illustrative 
cases.  Our  division  of  these  topics  into  three  separate  categories  is 
not  intended  to  imply  that  they  are  mutually  exclusive;  it  has  been 
effected  for  convenience  alone. 

The  problem  of  determining  how  a  homogeneous  expanse  of  gas 
behaves  as  it  approaches  equi liu  'urn,  i.e.  "relaxes,"  is  perhaps  the 
simplest  problem  in  the  kinetic  theory  of  gases.  The  problem  is  of 
interest  here  because  it  focuses  attention  on  the  collision  operator. 
°erhaps  the  simplest  description  of  the  relaxation  process  is  that  ob¬ 
tained  from  the  homogeneous  isotropic  Krook  model  (2.?5).  Thus  if  fD 
is  the  value  of  the  distribution  function  at  t-0  ,  then 

ut)  -  N0F,+  [f.  -  H,  "•<rt  (2.36) 

where  F0  is  given  in  (2.34).  The  characteristic  time  (^''appearing  in 
(2.36)  is  generally  referred  tc  as  a  relaxation  time.  Due  to  the 
complexity  of  the  more  sophisticated  models,  a  study  of  relaxation  via 
analytic  solution  has  generally  not  beer;  ccnieved.  While  the  simple 
Krook  relaxation  time  may  often  be  a  sufficiently  accurate  estimate  of 
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relevant  time  scales,  it  is  easily  recognized  that  the  description  of 
a  relaxation  process  by  a  single  parameter  could  often  be  misleading. 

For  example,  the  rate  at  which  a  given  distribution  becomes  isotropic 
in  velocity  space  if  initially  anisotropic,  could  be  significantly 
different  than  the  rate  at  which  it  relaxes  to  a  Maxwellian. 

One  method  of  estimating  such  rates  without  actually  solving 
the  kinetic  equation  is  to  single  out  and  consider  a  single  "test" 
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particle  in  the  gas.  This  procedure  has  been  employed  by  Spitzer, 

who  analyzed  various  aspects  of  the  relaxation  of  the  electron  and  ion 

components  of  an  ionized  gas,  such  as  (i)  removal  of  angular  anisotropy, 

(ii)  energy  exchange,  (iii)  loss  of  energy  of  a  particle  by  "dynamical 
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friction."  Bohm  and  Aller  have  similarly  presented  a  detailed 
analysis  on  the  relative  importance  of  electron-electron  collisions  in 
establishing  the  velocity  distribution  of  electrons  in  gaseous  nublae 

4 

and  stellar  atmospheres.  Montgomery  and  Tidman  perform  a  test  particle 
analysis  by  assuming  all  particles  except  the  test  Darticle  have  a  known 
(equilibrium)  distribution.  The  kinetic  equation  (in  this  case  Fokker- 
Planck)  is  then  "linearized"  about  the  test  particle  "distribution," 
and  velocity  moments  of  the  linear  equation  are  obtained.  Relaxation 
times  are  then  obtained  by  defining  these  as  the  ratio  of  the 
velocity  moment  M  in  question  to  its  time  derivative,  i.e., 

M 

i't  I  =  - 7— T  - 

The  relaxation  times  obtained  from  the  test  particle  approach 
generally  depend  on  the  initial  speed  of  the  test  particle,  the  relative 
temperatures  of  the  species  present,  and  the  relative  masses.  For  a 
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two  component  ionized  gas  with  equal  electron  and  ion  temperatures,  it 
is  found  that  electrons  become  isotropic  primarily  through  collisions 
with  ions;  collisions  with  electrons  play  a  small  role  in  ion  relaxa¬ 
tion,  and  the  relaxation  to  equilibrium  of  an  isotropic  electron  dis¬ 
tribution  is  primarily  due  to  encounters  with  other  electrons. 

Although  tie  qualitative  conclusions  reached  in  a  test  particle 

treatment  should  generally  be  correct,  they  do  not  display  all  of  the 

information  available  in  the  kinetic  equation.  Hence,  in  lieu  of  an 

analytic  solution  of  the  kinetic  equations,  several  authors  have 

presented  numerical  treatments  for  various  situations,  MacDonald, 
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Rosenblutn,  and  Chuck  have  presented  a  numerical  solution  of  the 
Fokker-Planck  equation  (2.7)  for  an  isotropic  electron  gas  imbedded  in 
a  positive  neutral izina  background.  They  assumed  an  initial  Gaussian- 
shaped  distribution  peaked  in  the  vicinity  of  the  speed  (2®/^ They 
found  the  time  required  for  the  distribution  to  come  within  a  few  per¬ 
cent  of  the  final  Maxwellian,  throughout  the  range  from  zero  energy  to 
several  times  the  average  energy,  is  about  ten  times  the  self-collision 
time  defined  by  Spitzer  (the  mean  time  required  for  a  thermal  particle 
to  eventually  suffer  a  9(T  change  in  direction  due  to  the  cumulative 
effects  of  many  small  angle  encounters  with  like  particles).  MacDonald, 
et  al_  also  found,  as  could  be  expected,  that  it  takes  considerably  longer 

to  fill  out  the  high  velocity  "tail"  of  the  Maxwell  distribution. 
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Recently  Wu,  Levans,  and  Primack  have  studied  numerically  the 
relaxation  of  a  two-component  plasma  with  initially  anisotropic  electron 
and  ion  temperatures,  and  with  initially  isotropic  (but  unequal)  compon¬ 
ent  temperatures ,  using  the  Lenard-Guernsey-Balescu  equation  (2.30), 
(2.31).  They  assumed  that  the  distribution  functions  maintain  a 
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Maxwellian  character  throughout  the  relaxation  process,  having  the 
form 
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The  results  are  in  qualitative  agreement  with  the  predictions  of  the 

test  particle  theory,  except  when  Qe/&{  >•/£>,  in  this  case,  the 

anisotropic  electron  temperature  relaxation  is  governed  by  collective 

47 

phenomena.  Since  Wu,  ext  al^  constrained  their  distribution  functions 

to  the  form  (2.37),  a  comparison  of  their  isotropic  relaxation  results 
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with  those  of  MacDonald,  et  a\_  would  not  be  fruitful.  We  can  only 
note  that  the  collective  effects  manifest  in  (2.30),  (2.31)  are  impor¬ 
tant  in  the  relaxation  process,  under  certain  conditions  of  anisotropy 
A  direct  comparison  of  the  Fokker-Planck  and  Lenard-Guernsey- 

Balescu  equations  has  been  achieved  in  another  numerical  relaxation 
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study,  performed  by  Dolmsky.  Dolinsky  solved  both  equations  for 
several  different  initial  conditions,  for  an  isotropic  electron  gas  in 
a  neutralizing  background,  A  comparison  of  the  solutions  showed  a 
difference  of  less  than  two  percent,  for  all  speeds  and  for  all  time. 
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5.  Transport  Phenomena;  Electrical  Conductivity 


Many  of  the  interesting  phenomena  in  kinetic  theory  involve 
systems  that  are  inhomogeneous  or  are  subject  to  external  fields.  As 
we  discussed  in  Chapter  I,  the  role  of  particle  interactions  in  such 
cases  may  or  may  not  be  important  relative  to  other  phenomena,  depending 
on  the  nature  of  the  system  under  consideration.  A  simple  and  yet 
interesting  illustration  is  the  phenomenon  of  electron  runaway,  which 
occurs  when  an  ionized  gas  is  subject  to  a  sufficiently  strong  electric 
field. 
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Kruskal  and  Bernstein  have  studied  electron  runaway  using  a 
transport  equation  with  Fokker-Planck  collision  operator.  For  simplicity 
they  neglected  electron-electron  collisions,  and  assumed  the  electron- 
ion  mass  ratio  to  be  zero.  Their  analysis  leads  to  a  decomposition  of 
velocity  space  into  three  regions,  for  electric  fields  greater  than  a 
critical  value.  In  the  first  of  these,  the  low  velocity  domain,  the  form 
of  the  electron  distribution  function  is  dominated  by  collisions  and 
hence  almost  isotropic.  The  second  regior,  one  of  intermediate  velocity, 
is  characterized  by  "quasi -steady"  flow  in  velocity  space,  for  which  the 
low  velocity  region  provides  the  source.  Lastly  there  is  a  high  velocity 
region,  fed  by  the  intermediate  region,  in  which  the  electrons  acceler¬ 
ate  or  "run  awav"  almost  freely  under  the  action  of  the  electric  field, 
iy  a  very  weak  diffusion  due  to  collisions. 

The  phenomenon  of  runaway,  like  the  relaxation  of  high  speed 


electrons  discussed  earlier,  reflects  the  rapid  decrease  of  the  Coulomb 
cross-section  with  increasing  relative  velocity.  It  is  apparent  that 
a  simplified  collision  model  that  does  not  take  into  proper  account  the 
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nature  of  the  interacting  oarticles  would  here  lead  to  erroneous 
results. 

As  we  noted  in  Chapter  I,  there  are  many  interesting  circum¬ 
stances  wherein  the  macroscopic  properties  of  a  system  exhibit  only 
small  variations  in  times  of  the  order  of  the  inverse  collision 
frequency  Vc  or  in  space  over  distances  of  the  order  of  the  mean  free 
path  Lc  './vc  .  It  is  then  possible  to  approximate  the  kinetic 
descript  /  a  fluid  description  treating  macroscopic  quantities  such 

as  density,  mean  velocity,  pressure,  etc,,  where  flows  are  linearly 
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related  to  the  generalized  forces  driving  them.  ’  ’  For  example  the 
electric  current  is  given  by  the  product  of  the  electric  field  and  the 
conductivity.  For  an  isotropic  system  the  conductivity  is  a  scalar; 
more  generally  it  is  a  tensor. 

Since  we  are  concerned  in  this  work  with  the  description  and 
effects  of  particle  interactions,  it  is  instructive  to  consider  the 
calculation  of  transport  parameters  briefly.  While  all  transport 
coefficients  are  sensitive  to  particle  interactions,  the  phenomena  of 
interest  here  can  be  illustrated  by  a  consideration  of  the  electrical 
properties  of  a  plasma.  In  the  following  discussion  we  will  assume  for 
simplicity  that  the  system  being  considered  is  free  from  magnetic  fields, 
temperature  gradients,  and  i nhomogenei ties .  We  will  further  assume  that 
the  applied  field  is  spatially  uniform;  i.e.,  that  A » A  p  where  X  is 
a  length  characterizing  the  field. 

Considerable  attention  has  been  turned  in  recent  years  to 
determining  the  conductivity  of  a  fully  ionized  gas.  The  subject  is  not 
only  of  interest  as  a  problem  in  kinetic  theory,  but  is  also  of  practical 
impot  cance  in  that  from  a  knowledge  of  the  a.c.  conductivity  one  can 
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compute  immediately  the  absorption  coefficient  for  radiation  in  a  plasma 

4  59  51 

and  hence,  by  Kirchoff's  law,  the  emission  properties.  *  * 

The  earliest  calculations  of  electrical  conductivity  were  based 

2  39 

on  phenomenological  considerations.  *  Thus  one  simply  assumes  that 
the  current  carrying  electrons  suffer,  on  the  average,  equal  accelera¬ 
tions  by  the  electric  field  and  decelerations  due  to  collisions.  Using 

39 

such  considerations,  Spitzer  calculates  a  conductivity  assuming  all 
current  to  be  carried  by  the  electrons,  and  neglecting  interactions 
between  electrons.  He  finds  the  conductivity  cr  to  be 


where  N0  and  »  are  the  electron  number  density  and  mass,  respectively, 

and  ve  is  the  electron-ion  collision  frequency.  We  can  easily  obtain 

similar  results  with  a  simplified  Krook  model.  The  d.c.  conductivity 

of  a  fully  ionized  gas  has  also  been  comouted  using  a  Fokker-Planck 
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collision  operator.  Spitzer,  et ,  aj_  ’  obtained  a  numerical  result  very 
close  to  that  given  by  (2.38)  with  vc  the  "self-collision"  frequency 
defined  earlier. 

As  we  noted  above,  a  significant  contribution  to  the  collisional 
processes  in  a  p’asma  arises  from  long-range  Coulomb  encounters,  and 
the  duration  of  these  encounters  is  quite  sensitive  to  the  relative 
speed  of  the  pc. tides.  While  a  Krook-type  collision  model  seems  satis¬ 
factory  for  a  fully  ionized  gas  subject  to  low  frequency  or  d.c.  fields, 
we  would  not  expect  such  a  simple  representation  of  particle  interactions 
to  suffice  for  frequencies  eo  of  order  or  higher. 

The  earliest  treatments  of  the  a.c.  response  of  a  plasma 
attempted  to  overcome  these  difficulties  by  employing  velocity  dependent 
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collision  frequencies  or  collision  operators  that  described  a  diffusion 

2 

in  velocity  space.  These  approaches  were  not,  however,  founded  on  a 

consideration  of  the  nature  of  the  interacting  particles. 

Several  authors  have  computed  the  impedance  2  -  CT"1  of  a 

plasma  employing  the  various  kinetic  theories  considered  earlier  in  this 

chapter.  Their  results  provide  an  effective  means  for  comparing  Lhe 

various  theories.  In  Fig.  2.3  we  have  displayed  the  quantity  RAe/Roc 

where  £  is  the  real  part  of  the  plasma  impedance  2  ,  as  a  function  of 

cj/ojf  .  The  results  are  given  for  a  fully  ionized  hydrogen  plasma 

with  A  *  loy  ,  where  A  is  given  by  N0  ,  and  are  based 
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on  a  similar  display  due  to  DeWolf. 

For  frequencies  co  well  below  the  collision  frequency  the 

purely  resistive  impedance  is  constant  and  the  results  of  the  simple 
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Krook  theory  agree  with  the  Fokker-Pl anck  calculations.  The  low 
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frequency  resistivity  has  been  computed  *  both  including  and  excludinq 
(Lorentz  gas)  encounters  between  electrons.  The  effect  of  including 
these  is  to  increase  the  low  frequency  resistivity  bv  a  factor  ~1.7 
as  is  evident  in  the  figure. 
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Bernstein  and  Trehan,  Robinson  and  Bernstein,  Kauffmann, 
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and  Shkarofsky  have  obtained  the  a,c.  plasma  impedance  using  a  Fokker- 
Planck  collision  operator.  Their  results  are  summarized  in  Shkarofsky, 
Bernstein,  and  Robinson. ^  Marshall^  performed  a  similar  analysis 
using  the  linearized  Boltzmann  collision  operator.  None  of  these 
authors  included  the  effects  of  internal  "self-consistent"  fields;  i.e„, 
they  did  not  include  the  Maxwell  equations  in  their  analysis.  The 
results  for  the  Boltzmann  and  Fokker-Planck  collision  operators  agree 
within  a  few  percent,  and  this  difference  is  likely  due  to  different 
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computational  procedures.  For  frequencies  below  vc  these  results  match 
the  d.c.  results,  as  is  evident  in  the  resistivity  diagram.  (For 
frequencies  not  small  compared  to  the  collision  frequency  the  impedance 
has  a  reactive  part^*^*^*^  reflecting  inertial  effects  of  the  con¬ 
ducting  charges.)  As  the  frequency  increases  past  vc  the  resistivity 
increases  and  eventually  approaches  a  constant,  independent  of  ^  .  For 
frequencies  well  above  ve  electron-electron  collisions  are  seen  to  be 
insignificant  as  compared  with  electron-ion  collisions, 

62 

When  the  field  frequency  exceeds  the  plasma  frequency  Scheuer 
has  argued  that  the  resistivity  should  decrease,  in  contrast  with  the 
Fokker-Planck  or  Boltzmann  results. When  eo>cu?  Scheuer  suggested 
..hat  the  maximum  effective  impact  parameter  should  decrease  from  the 
Debye  length  \rJ/’oo?  to  the  length  v~T/oo  .  At  distances  larger  than 
\rT/co  ,  encounters  do  not  contribute  to  the  resistivity  since  they 
are  much  longer  in  duration  than  the  oscillations  themselves.  Dawson 
and  Oberman63’^  computed  the  high  frequency  impedance  of  a  plasma  using 
the  simple  Vlasov  equation  including  the  internal  electrostatic  field. 

go 

Their  results  (see  Fig.  2.3))  agree  with  Scheuer's  reasoning  for 
co>oo p  and  join  the  Fokker-Planck  and  Boltzmann  results^*^  for 
co  <  ^0  f  .  They  observed  a  slight  bump  in  the  resistivity  near  co~coy 
which  they  attribute  to  the  generation  of  longitudinal  plasma  oscilla¬ 
tions.  For  very  low  frequencies  their  resistivity  does  not  decrease, 
in  contrast  with  the  collisional  treatment. 

This  latter  difference  in  the  predictions  of  the  two  treatments 
miqht  be  interpreted  by  reasoning  as  follows.  For  frequencies  below  the 
plasma  frequency  the  dielectric  response  of  the  plasma  is  fast  compared 
with  the  period  of  the  imposed  oscillation.  For  frequencies  in  the 
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range  Vc<<  to  <  ajp  ,  collisions  are  too  slow  to  affect  the  plasma 
response  and  a  description  incorporating  only  dielectric  effects,  i.e,, 
the  Vlasov  equation,  yields  results  that  are  insensitive  to  frequency. 
For  lower  frequencies  the  collisions  become  important  and  the  response 
is  frequency  sensitive.  Evidently  there  is  a  range  of  frequencies 
below  cop  wh .  a  the  collisional  and  collective,  or  dielectric,  des¬ 
criptions  produce  similar  results. 

65 

Oberman,  Ron,  and  Dawson  have  computed  the  high  frequency 

conductivity  of  a  fully  ionized  plasma  by  solving  the  first  two  members 

(2.13),  (2.14)  of  the  BBGKY  hierarchy  using  a  method  due  to  Guernsey.*’*’ 

We  note  that  the  Bogoliubov  hypothesis  was  not  employed;  i.e.,  the  two 

particle  correlation  function  was  allowed  to  very  on  the  same  time 

scale  as  the  one-particle  distributions.  The  results  of  Oberman,  et 
65 

al  are  in  complete  agreement  with  the  predictions  of  the  much  simpler 
Vlasov  treatment.63*6^ 

In  concluding  this  section  we  note  some  other  computations  of 

plasma  transport  parameters  for  the  interested  reader.  In  references 

57  through  60  the  low  frequency  (  u><  oj*)  thermal  diffusion  and 

conductivity,  and  the  viscosity  have  been  computed  as  well  as  the 

electrical  conductivity,  for  a  plasma  having  small  temperature  and 

density  gradients  and  immersed  in  a  constant  uniform  magnetic  field. 

Kivelson  and  Dubois*’7  have  found  the  electrical  conductivity  for  finite 

fift 

wavelengths  using  the  kinetic  equation  (2.30)  of  Lenard,  et  al_.  Berk 
nas  obtained  the  conductivity  for  finite  wavelengths.  His  approach  was 
similar  to  that  of  Dawson  and  Oberman.*’3,6^  Oberman  and  Shure^9  used 
the  first  two  BBGKY  equations  as  in  ref.  65  to  compute  the  high  fre¬ 
quency  conductivity  with  a  magnetic  field.  The  high  frequency 
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electrical  conductivity  has  been  computed  quantum  mechanically  by  DuBois, 

Gil  insky,  and  Kivelson^0  and  by  Ron  and  Tzoar.^  Oberman  and  Ron^ 
extended  this  work  to  include  a  magnetic  field.  The  results  are  in 
agreement  with  the  classical  descriptions. 

We  summarize  briefly  the  effects  of  finite  wavelengths,  and 
magnetic  fields.  Fo>'  finite  wavelengths  it  is  convenient  to  refer  to 
the  phase  speed  of  tha  wave,  =  oj/'k  where  h  *  2tt/\  ,  For  or  i/y. 
or  less,  and  <*>+ *  ,  the  principal  contribution  to  the  conductivity 

is  electrostatic.®®^0  For  greater  speeds  collisional  effects  predom¬ 
inate.  At  high  frequencies,  eo>cLjf)  ,  collisions  are  unimportant.  In 
the  presence  of  a  magnetic  field*’0  the  conductivity  in  the  direction  of 
the  field  is  unaffected.  In  contrast  the  transverse  comoonents  decrease 
with  increasing  field  strength  until  in  the  limit  of  infinite  field 

strength  no  current  flows  across  the  field. 

73 

Finally,  Klevens,  Primack,  and  Wu  have  computed  the  a.c. 
conductivity  for  using  the  Lenard-Guernsey-Balescu  equation 

(2.30).  Two  specific  cases  are  considered  in  detail:  in  the  first,  the 
unperturbed  plasma  has  different  electron  and  ion  temperatures;  in  the 
second,  the  unperturbed  plasma  is  characterized  by  a  relative  drift 
between  electrons  and  ions.  For  the  first  case  they  find  that  for 
/ &’  or  \oH  ,  the  real  part  of  the  conductivity  becomes 

negative.  For  the  second  case  they  find  that  if  the  electron  drift 
speed  exceeds  1.37  times  the  electron  thermal  speed,  and  if  Bt'  >  , 

the  conductivity  is  again  negative. 

I 
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6,  Cillisional  Effects  on  Small  Amplitude  Plasma  Oscillations 

Until  recently,  most  studies  of  plasma  oscillations  have  been 
concerned  with  relatively  high  temperatures  and  low  densities,  as  in 
thermonuclear  devices,  or  with  very  weakly  ionized  systems  such  as  the 
ionosphere.  For  systems  in  the  first  category  the  collision  frequen¬ 
cies  are  generally  very  small  compared  with  the  oscillatory  frequencies 
of  interest,  and  this  is  used  as  a  basis  for  disregarding  collisional 
effects.  For  systems  in  the  second  category  collisions  with  neutrals 
often  predominate,  and  a  simple  Krook-type  model  is  employed  to  account 
for  these. 

In  recent  years  considerable  experimental  attention  lias  been 
given  to  fully  ionized,  relatively  low  temperature  plasmas  for  which 
the  foregoing  collisionless  or  simple  collision  model  assumptions  are 
thought  to  be  unsound.  In  addition,  in  the  study  of  high  temperature 
unstable  plasmas  it  has  been  recognized  that  an  inclusion  of  even  very 
weak  collisions  can  have  a  significant  effect  on  the  growth  rate  of  the 
instability.  For  these  reasons  there  has  appeared  an  incentive  to 
treat  collective  phenomena  including  the  effects  of  particle  inter¬ 
actions. 

In  the  following  discussion  we  present  a  brief  review  of  some 
recent  work  on  collisional  p  ects  in  plasma  collective  behavior.  Since 
relatively  few  theoretical  results  are  available,  the  work  summarized 
here  should  be  considered  as  a  first  step  in  the  direction  of  under¬ 
standing  these  phenomena.  In  keeping  with  the  objective  of  the  present 
treatise,  our  emphasis  is  on  the  nature  and  description  of  collisional 
effects.  The  multitude  of  possible  collective  modes  that  a  plasma  may 
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support  often  mokes  a  generalization  of  specific  results  very  difficult, 
and  few  attempts  to  do  so  are  made  here.  Since  a  considerable  effort 
has  gone  into  the  analysis  and  classification  of  collisionless  oscilla¬ 
tions,  the  reader  having  more  than  a  passing  interest  in  collective 

phenomena  would  probably  benefit  by  consulting  the  collisionless  litera- 

x-  x  56.74-83 
ture  first.  * 

In  the  analysis  of  small  amplitude  plasma  collective  phenomena 

it  is  frequently  convenient  to  Fourier-Laplace  transform  the  governing 

equations,  together  with  the  Maxwell  equations  for  tho  electromagnetic 

field,  and  then  to  solve  the  transformed  equations  for  the  internal 

electric  field,  from  which  ail  other  field  quantities  may  be  deU;  ■  i  J, 

83 

The  result  is  then  displayed  in  the  form 


E 


f  *• « 


(2.39) 


where  f  is  the  plasma  dielectric  tensor  (or  constant  for  isotropic 
systems),  the  elements  of  §*  are  the  cofactors  of  their  counterparts 
in  6  ,  and  tl*)  1*  a  vector  incorporating  the  initial  conditions.  The 
dependence  of  the  electric  field  is  given  by  the  inverse  Laplace  trans¬ 
form  of  (2.39),  and  since  £*  and  are  entire  functions  cf  S  and  & 

83 

for  many  interesting  cases,  one  is  usually  interested  in  the  zeroes 
of  the  determinant  }6  I  .  Thus  setting  this  quantity  equal 

to  zero  yields  a  relation  between  the  wave  vector  k  and  the  Laplace 
variable  $  =  ieo  f-Y  ,  and  hence  one  estimates  the  growth  or  decay  rates, 
etc.  of  various  collective  modes  (it  should  be  noted,  however,  that  the 
dispersion  relation  l£|  =  ©  does  not  necessarily  imply  a  one-to-one 
correspondence  between  frequency  and  wavelength  ). 
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In  the  following  review  the  reader  will  note  in  some  cases  a 
remarkable  similarity  between  the  colli sional  effects  on  plasma  collec¬ 
tive  behavior,  and  the  collisional  effects  on  transport  phenomena 
discussed  above.  This  is,  of  course,  mere  than  fortuitous.  It  can  be 
shown7^*®^  that  the  longitudinal  and  transverse  dielectric  and  conduc¬ 
tivity  tensors  are  related;  for  S-  iw  we  have7^ 

-  I  +  1*i  (i,  )/*? , 

6jt~)  =  I  ¥ 
r  *  >*•  i 

with  r  the  unit  dyadic. 

In  an  early  attempt  to  treat  longitudinal  plasma  oscillations 
including  collisions,  Bhatnagar,  Gross,  and  Krook'  employed  the  simple 
collision  model  discussed  earlier.  They  treated  a  one  component  plasma 
consisting  of  electrons  with  fixed  ions  with  no  external  fields,  and 
neglected  collisions  between  electrons.  (The  assumption  of  fixed  ions 
implies  the  wave  frequency  is  large  compared  to  the  ion  plasma  frequency.) 
Their  results  may  be  summarized  as  follows:  (i)  For  wavelengths  long 
compared  to  the  Cebye  length  or  the  mean  free  path  a  small  change  in  the 
oscillation  frequency  was  observed  as  the  collision  frequency  varied 
from  zero  to  infinity;  the  damping  was  slow  (i.e,,  S/co«l)  and  reached 
its  maximum  when  the  collision  frequency  equalled  the  plasma  frequency, 
(ii)  For  wavelengths  shorter  than  both  the  Debye  length  and  the  mean 
free  path  the  damping  was  heavy  and  was  primarily  electrostatic,  or 
Landau  damping. 

or 

Lenard  and  Bernstein  treated  the  problem  studied  by  Bhatnagar, 
37  38 

et  al  *  using  a  pseudo  Fokker-Planck  collision  operator  designed  to 


represent  a  diffusion  in  velocity  space,  and  which  conserved  electron 
number  density  and  yielded  the  Maxwell  distribution  for  the  equilibrium 
state.  Their  velocity  dependent  "diffusion  coefficients"  however, 
increased  with  velocity  in  contrast  with  the  true  Fokker-Planck  coef- 

37 

ficients.  Their  results  are  in  general  agreement  with  Bhatnagar,  et  al_. 

87  88  89 

Comisar,  Gorman  and  Montgomery,  Burgers,  and  Wu  and 

90 

Kievans  have  treated  collisional  damping  of  longitudinal  electron 

oscillations  in  a  one  component  plasma,  including  both  electron-electron 

87 

and  electron-ion  collisions.  Comisar  used  the  linearized  Fokker-Planck 

88  66 

collision  operator,  Gorman  and  Montgomery  used  Guernseys  reduction 

89 

of  the  first  BBGKY  equations,  Burgers  solved  a  Boltzmann- like  equation 

with  the  Debye  potential  replacing  the  Coulomb  potential,  and  Wu  and 
90 

Kievans  approximated  the  first  two  BBGKY  equations  and  then  employed 

a  uuernsey-like  reduction.  All  of  these  authors  obtained  similar 

results,  which  were  restricted  to  weak  collisions  and  long  wavelengths. 

The  results  may  be  summarized  as  follows:  (i)  A  wavelength-independent 

damping  constant  was  found  for  electron-ion  collisions,  (ii)  a  damping 

constant  porportional  to  k1  was  found  for  both  electron-electron  and 

electron- ion  collisions,  (iii)  collision  damping  dominated  Landau  damping, 

(iv)  electron-ion  collisions  dominate  the  damping,  and  (v)  a  small, 

wavelength  independent  correction  to  the  oscillation  frequency  was  found. 

In  each  of  the  first  two  cases  the  damping  constant  V  was  found  propor- 

39 

tional  to  the  respective  collision  frequencies  riven  by  Spitzer.  The 
work  of  Comisar^  has  been  extended  by  Buti  and  Jain^  to  treat  high 
frequency  transverse  plasma  oscillations.  Their  results  are  essentially 
the  same  as  Comisar' s. 

The  collisional  damping  of  electron  plasma  oscillations  is 


-  57  - 


easily  described  in  terms  of  momentum  transfer  out  of  the  collective 

modes,  due  primarily  to  electron-ion  collisions,  and  the  damping 

increases  with  increasing  electron-ion  collision  frequency.  In  contrast, 

we  might  expect  the  effects  of  collisions  on  low  frequency  ion  waves  to 

be  somewhat  different,  since  momentum  transfer  to  electrons  is  small, 

92 

Bhadra  and  Varma  have  investigated  collisional  damping  of  longitudinal 
ion  waves  using  a  simple  Kro^k  model  and  neglecting  ion-electron  colli¬ 
sions.  For  equal  electron  and  ion  temperatures,  the  damping  decreased 
moiiotonical ly  with  increasing  collision  frequency.  Their  interpreta¬ 
tion  of  this  result  is  that  since  collisions  do  not  transfer  momentum 

out  of  the  wave,  their  only  affect  is  to  enhance  the  propagation. 

93 

Kulsrud  and  Shen  have  investigated  the  propagation  c.  ion 
waves  using  a  Fokker-Planck  collision  operator  in  the  limit  of  weak  ion- 
ion  collisions.  They  found  the  spatial  damping  to  decrease  with  increas¬ 
ing  collision  frequency  as  with  the  time  damping  treated  by  Bhadra  and 
92 

Varma,  and  calculated  the  relation  between  wave  speed  and  collision 

frequency  for  comparison  with  experiments  on  ion  waves  performed  by 
94 

Motley  end  Wong.  Their  results  are  in  fair  quantitative  and  qualita¬ 
tive  agreement  with  the  experimental  results,  but  they  suggest  this  may 
be  only  fortuitous  since  they  attempted  to  extrapolate  a  time-damping 
theory  to  explain  spatial  damping  lengths. 

We  have  so  far  been  concerned  with  waves  in  isotropic  plasmas- 
42  95 

Liboff  and  Oppenheim  have  treated  longitudinal  electron  plasma 
oscillations  in  the  presence  of  a  constant  uniform  magnetic  field. 

Liboff  used  a  Krook  model  to  represent  collisions,  while  Oppenheim 
employed  a  pseudo  Fokker-Planck  collision  operator  similar  to  that  used 

fifi 

by  Lenard  and  Bernstein.  The  two  treatments  gi»e  similar  results 
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for  long  wavelengths  and  low  temperatures,  and  for  magnetohydrodynamic 
modes,  in  the  absence  of  the  magnetic  field.  The  results  differ,  how¬ 
ever,  in  the  case  of  "microscopic  Lanrior  resonance"  modes,  in  the  para¬ 
meter  range  where  wavelength  is  much  longer  than  both  the  Larmcr 

42  . 

resonance  and  the  collision  length.  Liboff's  Krook  model  gives  an 

infinite  number  of  damped  Lafmor  modes  only  at  propagation  precisely 

perpendicular  to  the  a-plied  magnetic  field.  Cppenheim's  model,  in 

contrast,  gives  an  infinity  of  Larmor  modes  at  arbitrary  directions, 

except  parallel  to  the  field.  Oppenheim  suggests  this  difference 

reflects  the  velocity-space  diffusion  property  of  his  collision  operator. 

The  damping  constants  found  by  Oppenheim  and  Liboff  were  quite  similar, 

being  proportional  to  the  collision  frequency  in  each  case. 

83 

It  is  wel1  known  that  small  amplitude  disturbances  of  a 
homogeneous  plasma  near  thermal  equilibrium  are  stable;  i.e.,  any  such 
disturbances  tend  to  decay  in  time.  In  addition,  this  inherent  stability 
is  not  affected  by  the  inclusion  or  exclusion  of  colli  si  onal  effects  in 
the  describing  equations,  or  by  the  imposition  of  a  uniform  magnetic 
fi°ld.  In  contrast  the  presence  of  currents  c r  spatial  gradients  is 
known  to  be  sufficient  to  induce  unstable  plasma  behavior. 

The  study  of  plasma  instabilities  is  a  relatively  new  field  but 
nevertheless  has  received  prominant  attention  in  regard  to  both  labora¬ 
tory  and  extra-terrestrial  phenomena,  prime  examples  being  the  contain- 

96  97  98 

ment  of  hot  plasmas  and  the  growth  mechanism  of  stellar  flares.  * 

Due  to  the  considerable  complexity  of  the  equations  employed,  the 

analysis  of  plasma  instabilities  has  genevally  been  restricted  to 

collis:onles$  treatment*. .  Only  within  th-:  past  two  years  have  attempts 

been  made  to  include  colli  si onal  effect  .i .  Wf  .ile  these  efforts  have 


been  few  in  number,  the  results  indicate  that  these  effects  may  have  a 
profound  influence  on  plasma  behavior.  Certainly  more  work  is  needed 
in  this  area. 

99 

In  an  early  treatment  including  collision  effects,  Kuckes 
analyzed  the  propagation  of  low  frequency  ion  waves  in  a  current  carry¬ 
ing  plasma  without  a  magnetic  field.  Using  a  simple  one-parameter 
collision  model,  he  showed  that  "collisional  effects  of  the  electrons 
can  lead  to  growth  mechanisms  for  these  oscillations,"  while  "the 

thermal  motions  of  the  ions  leads  to  a  damping."  More  recently 
100  93 

Bhadra,  and  Kulsrud  and  Shen  have  reported  studies  of  ion  acoustic 
waves,  employing  Fokker-Planck  collision  operators  in  an  iterative  weak 
collision  analysis.  Bhadra10^  treated  waves  propagating  parallel  to  a 
strong  magnetic  field  with  a  perpendicular  density  gradient,  and  Kulsrud 
and  Shen  assumed  a  homogeneous  plasma  with  a  small  external  electric 
field.  Bhadra  found  electron-electron  collisions  to  have  a  destabili¬ 
zing  effect,  while  electron-ion  'ollisions  tended  to  stabilize.  Kulsrud 
and  Shen,  in  contrast,  observed  electron-ion  collisions  to  decrease  the 
critical  current;  electron-electron  collisions  had  negligible  effect. 
Bhadra^  also  used  a  simple  Krook  model  for  purposes  of  comparison; 
he  found  only  a  slight  difference  in  growth  rates  under  some  conditions. 


III.  AN  EXPANSION  THEOREM  FOR  THE  LINEARIZED  FOKKER-PLANCk  EQUATION 


1.  Properties  of  the  Equation 


In  the  first  three  sections  of  this  chapter  certain  spectral 
properties  of  the  collision  operator  are  established.  While  these 
properties  (apart  from  trie  reality  of  the  spectrum)  are  not  necessary 
for  the  later  development  of  the  expansion  theorem,  they  are  both  useful 
by  themselves,  and  enable  certain  conclusions  to  be  drawn  regarding  the 
final  form  of  the  expansion. 

For  our  purposes  it  will  prove  convenient  to  write  the  Fokker- 
Planck  equation  in  the  Landau  form  (2.9); 


where 


is  a  positive  constant,  g  =  jr_vr'  ,  and  X  is  the  unit  dyadic. 

It  is  not  difficult  to  show  that  the  Maxwell  distributions  }  F* 
satisfy  (3.1)  for  ,  j  =  A)  6  .  In  the  vicinity  of  equili¬ 
brium  we  may  write  +  •  Neglecting  terms  quadratic 

in  f*  we  then  obtain  from  (3.1)  the  linear  equation 
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We  will  refer  to  the  quantity  as  the  perturbation  from  equili¬ 

brium. 

Equation  (3.2)  as  it  stands  is  ‘n  fact  a  pair  of  coupled  equa¬ 
tions  for  -f.  and  fB  .  Due  to  the  uite  small  value  of  the  electron- 
ion  mass  ratio  the  equations  are  however  only  very  weakly  coupled. 

Thus  for  example  me  effect  of  the  ion  perturbation  on  the  electron 
pirturbation  is  small  when  compand  with  the  effect  of  the  ions  and 

electrons  in  the  unperturbed  equilibrium  distributions.  Tn  addition, 

2 

it  can  be  shown  that  in  the  approximation  ,  he  ions  act 

like  a  single  component  gas.  In  the  following  we  will  consider  the 
equation  for  he  ions.  The  treatment  of  the  electron  equation  is  quite 
similar,  and  the  modifications  necessary  for  this  case  will  be  indica¬ 
ted  later.  We  have  then, 


f: 


Fj^L  =  2.  .  (js^'F^F1  r>±  ifi'7  n  <3-3> 

* 

We  will  for  convenience  drop  the  subscript  "i"  from  and 
.  If -f  satisfies  the  conditions 


(/-f  —0 

v^o  1 


^  F„  H  =  ° 


(3.4) 


it  is  possible  to  show  that  (3.4)  conserves  number,  momentum,  and 
kinetic  energy  densities. 

Introducing  in  (3.3),  we  find 

Multiplying  (3.5)  by  ^  and  integrating  over  tr  we  have,  after  a  parts 
integration, 
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The  second  term  on  the  right  in  (3,6)  vanishes  provided  satisfies 
the  second  of  conditions  (3.4)  and 


xr~*>o 


(3.7) 


Assuming  these  conditions  hold  we  exchange  y  and  vr'  in  (3.6),  noting 
that  Q(£t  £')  -  Q  (}£\  £ )  .  We  add  the  result  to  (3.6),  obtaining 


(3.8) 


ffjWVF„F4$ $7. 

Since  Q  is  a  real  positive  quadratic  form  it  follows  that  the  right 
side  of  (3.8)  is  real  and  positive  or  zero.  Hence  &v\^=0and  A  ^  O 

3 

Employing  standard  methods  we  can  find  from  (3.8)  the  most 
general  form  of  ^  when  A  -0  ; 


Q  / lT \  -  Q  \rZ-t  CL  'V"  -j-  Cl 

j  o  \  ++•  J  ^  j  ‘v'u  ^  3  y 


(3.9) 


with  a  and  arbitrary,  but  necessarily  independent,  constants 


2.  Expansion  in  Spherical  Harmonics 


From  (3.5)  we  have 
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To  perform  a  parts  integration  on  the  second  term,  consider  the  quantity 


(3.10) 


with  Fm  -  ,  we  have 

C F*//^ ' *  -Z*;jr'Ff/  *  Also  JT'*  9  *  Jf  •  ?  .  so  (3-10)  gives 


•  (3J1 

The  first  term  on  the  right  in  (3.11)  vanishes  if  ^  satisfies  (3.4). 
Using  the  relations 


Q  =  Hv  ^  ^ •  Q  =  -  2fl  (3  )> 

and  introducing  a  dimensionless  time  f  and  dimensionless  velocity  c  in 
(3.3), 


y(-  ) 


c 
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the  kinetic  equation  takes  the  form 


~  Hrr A  e  c  fe)  * 


(3.12) 


The  time  T  is  measured  in  units  of  the  "Spitzer  self-collision  time,"^ 
and  c  -  Is/  is  in  units  of  the  rms  thermal  speed. 

Equation  (3.12),  in  three  dimensions,  may  be  replaced  by  a  set 
of  uncoupled  equations  in  one  dimension  by  introducing  the  spherical 
harmonic  expansion 


to  A 

^  Z  {)*'**'  (C)  Xl 

K  A-o  +«--J 


We  find  (see  appendix  A) 


c/V  e 


*4 


jc Pc'e-‘'\;k-£'l  - 
(jVe-£'l|£-f'l  " 


•f,"1 

T  1 fe_  VWC 


4" 

HwT 


where 


(3.13) 

(3.14) 

(3.15) 


R 


* 


J\4*i 

t 


-c 
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with  ef-f  £0  *  ln~,,x  (  e/xtxpC-x1)  .  The  expressions  for 

and  5^  were  found  by  Rosenbluth  et  al_.  ,5  using  a  different 
method,  for  the  axially  symmetric  case  iv\  ■<*.  o . 

Combining  (3.12)  through  (3.15)  we  find 


X 'i* 


-  c 


-p1 


Z/+I  ' 


-‘ifr*.,!- 


—  — — —  -2-  • 
Z-  *Ul-l  t*£ 
4*. 


.c  e 


(3-16) 

^  v#  » 


Performing  the  indicated  angular  differentiations  in  (3.16)  and  then 
employing  the  orthogonality  property  of  the  spherical  harmonics,  we 
obtain  the  uncoupled  equation 

-'U**.  -  T"  J-fer  +  [T"'  +  - 


_  AU+O 


1+0  t'  2  ?„W  -/Ct+Ool 

7^  T  ^  ~  \  Kt  +  (  * 


c  sj  +  s„j 


(3.17) 
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The  index  w\  is  clearly  superfluous  and  will  be  deleted  in  the  following. 
Performing  the  primed  differentiations  in  (3.17)  we  find 

+  {  4se-*%  0-  FrO  f 

WMO+V  p«V‘W' V<0  - 

(2/+l)(2A+l)  I  d 


_  Jk. 

ZUf 


tcUa-i) 

Hi 

Jk— 

2i+l 


*  't: 

-  Cic'e-^l T'fLt'). 

Utl  J  u  4 


(3.18) 


For  boundary  conditions  we  will  use  conditions  (3.4),  which  were 
obtained  from  the  conservation  laws.  Although  (3.7)  is  stronger  than 
the  first  of  conditions  (3.4)  we  will  see  below  that  the  solutions  of 
(3.18)  which  satisfy  (3.4)  also  satisfy  (3.7). 

3.  Spectrum  of  the  Radial  Equation 

If  we  introduce  the  transformation 

~  c  ^  (c 


(3.19) 
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we  obtain  from  (3.18)  the  formally  self-adjoint  equation 


sL 

dc 


Qtl  +  3.20) 

4 


with 

f“^)=  771  ecP£)-  2^e‘el; 

=  (h  '?)  ^frfVy  +  (i-gk  -  ■p)e";l- 


Transforming  (3.4)  via  (3.19)  we  have 

-4>U  e~c^z  =0  (3.21) 

C  00  d  C 

"ft0  %  -  °  •  (3.22) 

We  will  later  show  (Appendix  B)  that  the  problem  (3. 20)- (3. 22)  is  self- 
adjoint. 

The  spectrum  of  (3.20)  is  that  set  of  numbers  such  that 

(3.20)  has  non-trivial  solutions  which  satisfy  (3.21)  and  (3.22).  We 
have  already  seen  that  the  must  be  real,  and  must  be  positive  for 
A't-Z.  and  positive  or  zero  for/<2. 
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We  can  find  the  spectrum  of  (3.20)  by  first  considering  the 
related  problem 


l]  f  [fy  f  (3.23) 

with  conditions  on  the  functions  identical  to  (3.21)  and  (3.22). 

Clearly  ?  ^  <JP/dc,  and  are  bounded  and  continuous  for  all  finite 
c  except  possibly  near  c-o  .  For  c « I  we  have 

Qjt  (c) ~  ^  -A (ah)  [fp.  '  f  . 

In  general  v;e  can  write 

P£)  ~  eT?  (  J*  e~c  ) 

'  o 

and  thus  ?(c)>o  for  all  c<<*>  .  It  follows  that  for  J  J  o  , 
(3.23)  has  a  regular  singular  point  at  . 

For  c  small,  (3.23)  has  the  asymptotic  solutions 

^(c.)  ~  e**\  C-*  (e«  o.  (3.2/1) 

The  first  of  these  satisfies  (3.22)  for  all  /  .  It  also  satisfies  the 
stronger  conditio.;  obtained  from  (3.7),  (c)  -  o  . 

The  second  solution  satisfies  neither  condition.  This  is  clear  for 
XJ  0  .  For  /(-  0  ,  the  second  solution  is  a  constant  which  cannot 
be  zero  since  the  solutions  (3.24)  are  linearly  independent. 

For  c  sufficiently  large  and  for  Jo  ,  (3.23)  takes  the 
asymptotic  form 


liillil li'tlt! iil  1*1' " . . . . . . .  "  . . . . . .t.H ililljl fllifl I H 
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I 


=a*i 


I 

I 


Cc 


We  find  (^^>0  ) 


)  -  A*  C  3A  <W  (V C*,X-  )  (3-25) 

with  ■=  . 

Given  /  and  ,  equation  (3.23)  has  only  one  solution  which 
satisfies  the  condition  at  c  -O  ,  This  solution  thus  contains  only  one 
arbitrary  constant  and  it  follows  that  A ^  and  in  (3.25)  are  not 
independent.  Whatever  the  relation  between  A*/  and  is,  (3.25) 
satisfies  (3.21)  for  all  positive  A^  .  Hence  the  spectrum  of  (3.20) 
contains  all  positive  A^  ,  for  each  X  . 

To  determine  the  spectrum  of  the  integrodifferential  equation 
(3.20)  with  (3.21)  and  (3.22)  we  note  that  K j  (*>£')  is  a  Hilbert- 
Schmidt  kernel  and  the  symmetric  integral  operator  in  (3.20)  is  conse- 

C 

quently  completely  continuous.  According  to  Weyl's  perturbation 
theorem, ^  the  addition  of  a  completely  continuous  symmetric  operator 
cannot  alter  the  continuous  spectrum  of  any  symmetric  operator  to  which 
it  is  added.  Since  (3.23)  is  self-adjoint  (see  Appendix  B)  it  is 
symmetric,  and  it  follows  that  the  spectrum  of  each  X  -component  of  the 
linearized  Fokker-Planck  equation  contains  all  positive  real  A^j  ,  For 
X  -  o  andw/-l  we  found  belongs  to  the  spectrum;  this  corres¬ 

ponds  to  a  shift  to  an  equilibrium  different  from  that  originally  postu¬ 
lated.  For  I  the  spectrum  is  (U+l  )-fold  degenerate,  corresponding 
to  the  (2X^1  )  different  spherical  harmonics  of  order  X  . 


/ 
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4.  Outline  of  the  Expansion  Theorem 

We  now  proceed  with  the  development  of  an  expansion  theorem 
based  on  (3,20).  Our  method  is  essentially  an  extension  of  the  theory 
due  to  Weyl  and  Levinson,^  to  include  singular  integrodifferential 
equations  with  Hilbert-Schmidt  kernels.  We  give  a  brief  outline  in 
this  section.  The  details  are  left  to  the  appendix. 

Since  the  eigenfunctions  of  the  kinetic  equation  are  bounded 
and  continuous  on  every  finite  interval,  it  is  natural  to  pursue  an 
expansion  formula  for  functions  u(c)  squat'e  integrable  on  the  interval 
/t;  os  c<oo  .  As  with  the  earlier  theory  we  first  establish  an  expan¬ 
sion  formula  on  a  finite  subinterval  S  of  i  .*  a  -  c  S  Q<<X} 
so  that  the  singularities  of  the  linear  operator  are  external  to 
£  .  The  expansion  formula  we  seek  is  then  obtained  by  taking  / ->  A 

in  a  suitable  manner.  In  the  following  we  will  mean  by  LA  the  integro¬ 
differential  operator  in  (3.20)  and  by  Lj  the  operator  obtained  when 
the  lower  and  upper  limits  of  the  integral  in  (3.2C)  are  replaced  by  cl 
and  t  ,  respectively.  In  the  followinq  the  index/  will  be  retained 
only  where  it  is  necessary  to  avoid  confusion. 

We  have  already  seen  that  P}  ?'  and  Cf  ai'e  continuous  on  <f 
and  that  Kf  Cc^)  is  bounded  and  integrable  on  the  square  k, 

i  8 

<k*c  ±b.  Tamarkin  has  shown  that  subject  to  these  conditions,  the 
solutions  of  i ;  f » -  A  (j>  which  satisfy  homogeneous  bondary  conditions 
at  c  -  cl  and  c  *  b  form  a  complete  orthogonal  and  normalizable  set  of 
eigenfunctions  on  £  ,  with  an  associated  denumerable  sequence  of 

reat  eigenvalues  Assuming  the  to  be  normalized,  the 

expansion  formula  on  ^  is  thus 
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*£)=  X  h*  (  ^  (3,26) 

W  =  0 

where  <a(c)  is  any  function  square  integrable  on  i  . 

We  now  use  the  Weyl-Levinson  theory  to  take  &  A  .  Since  the 
subsequent  development  of  the  expansion  theorem  is  in  every  respect  a 
duplication  of  the  earlier  theory,  we  will  display  only  the  salient 
features. 

Given  A  ,  the  most  general  solution  of  Lj  p  =  $  is  a  linear 

combination  of  the  two  linearly  independent  solutions,  My  h  •  h  • 
Thus  we  can  write 


(c> =  rs„i  i  (V-^)  f  <k(‘, -»;»>  (3-27) 

where  fjn|  and  are  complex  constants.  With  (3.27;,  (3.26) 


becomes 


t  i  .  /  OL. 


(3.28) 


W-©  J.jk  = 


Following  Levinson  we  define  an  Hermitian,  positive  semidefinite 
matrix  ^  ,  called  the  spectral  matrix,  with  elements  which  con¬ 

sist  of  step  functions  with  jumps  at  the  eigenvalues  A ^  given  by 

N*.  (**• +  0  )  '  fat  (x‘*  *°>  -  . 

Let  \  +o)  *  fy  (A)  »  and  let  be  the  zero  matrix.  We  use  the 

spectral  matrix  to  replace  the  infinite  series  in  (3.28)  by  a  Lebesgue- 


Stieltjes  integral 


oo  2 
—  0^  ^ 


(3.29) 
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where 


ck. 


As  £->A  (that  is,  o.-^-o  ,  ^  approaches  a  limit  matrix  ^  . 

To  find  let  X  ,  cj^o,  and  let  be  a 

solution  of  L*  t-AP  satisfying  the  homogeneous  boundary  condition 

Snsl*  £(cl)  f  Pi\N0(  ?M  KYO  -0) 

and  similarly  let  be  a  solution  of  the  same  equation 

satisfying 

OA  (i  tfl)  f  Ws/S  PA)  =o  . 

Clearly  wy  »  AO  and  silTlilarly  for  • 

As  <\-*o  and  i-*  o ©,  ^  and  approach  limiting  values  in  the  complex 
rv\  plane  denoted  respectively  bvy^A)  and  >*00$).  These  limiting 
values  are  clearly  determined  by  the  behavior  of  <£,  and  ^  for  small 
and  large  c  ,  for  A  complex. 

For  (}  and  ^  to  be  linearly  independent  it  is  necessary  and 
sufficient  that  their  Wronskian  equal  a  nonzero  constant,  say  one: 

p  [<H'-  <kv]  '  i-  (3-30) 

This  last  will  be  satisfied  if  and  ^  satisfy  the  conditions 

A,<*)  =  4x(s^)  -  o*<r 

P&)  4i'  f^A  )  *  -  c*»<r  P6)  <r  , 

where  2  is  an  interior  point  of  <£  and  o  i<r<7T.  These  conditions 
are  also  sufficient  to  ensure  that  ^  ^  are  entire  functions  of 


-  73  - 


for  each  fixed  c  on  <£  (this  follows  from  Tamarkin®). 

7 

With  these  properties  secured  we  can  find  the  limit  values  m# 
and  r*"  and  hence  the  limit  matrix  a*  ,  whose  elements  are  given  by 


/Wk  (X  ^  -fa'k  (}) 


where 


i 

-  W«o 


M,j  £ ) 


J.  Wa  t 
^  -  Moo 


M»o 
^0  -  V*oo 


To  find  the  M/k  we  need  asymptotic  forms  of  and  for 
large  and  small  c  .  These  are  given  by  (3.24)  and  (3.25),  as  may  be 
verified  by  direct  substitution.  Taking  to  be  asymptotic 

respectively  to  for  c  small,  we  apply  the  homogeneous  boun¬ 

dary  condition  to  and  then  cake  <^>0  to  find 


~  oo  (/ i tc) 


YA0  -  ~  Cot*  (A  -o)  m 


(3.32) 


Thus  for  /yf  O  only  Mu,  can  have  a  nonzero  imaginary  part  and  conse¬ 
quently  only  (p^  will  contribute  to  the  expansion  formula  (3.29).  When 
X-0  both  solutions  are  regular  atc^o  and  the  limit  matrix  is  not 
determined  until  we  specify  e<  .  The  boundary  condition  (3,22)  dictates 
the  choice  o(-0. 

For  c  large  we  tal  i 
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&  -  Vci/¥  ^  ’  ^/l/i7rl/i  \j  >  o 

and  find  by  integrating  (3.30): 


/4.J  vjj  44  6vf 


where  £#  is  a  constant  of  integration.  Applying  the  homogeneous  boun¬ 
dary  condition  to  2%  and  then  taking  is  ->*>  with 

VsA>0  we  find 


3/ exp  1^6^ 

=  - b - 


(3.33) 


Combining(3.31 )  through  (3.33)  we  have  finally 


c/f  A*)  *  dfaiQt)  =  ^>0-(3.34) 

Since  the  spectrum  is  empty  for^jco,  ^*?)is  constant  on  this  range 
The  expansion  formula  (3,29)  becomes 

OC 

*(‘)  =  (3.35) 

*  o 

with 


(A 


The  expansion  converges  in  the  mean  for  all  functions  a(c)  square 
integrable  on  (o,oo).  If  the  spectral  function  p  is  not  continuous 
at  ^  - o ,  this  point  will  contribute  to  the  integral  in  (3.35). 

We  return  to  the  description  of  perturbations  from  equilibrium. 
If  is  square  integrable  in  velocity  space,  the  crom 
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(3. 19)  and  (3.35)  we  have 

00 

e-£lf (i,r)  -  2  pp(X<)3U.^,  e' A"jr  (3.36) 

<"*  0 

wi  th 


and 


^ f *, A./.)  -*-'*-*  *6 6  A* )  tf7*.  *>  • 

The  functions  %jf  correspond  to  <f>z  in  (3.35)  and  are  the  solutions 
of  (3.20)  satisfying  (3.22). 

We  have  defined  the  density,  mean  velocity,  and  kinetic  tempera¬ 
ture  of  the  ion  gas  as  being  proportional  respectively  to  the  first 
three  moments  of  the  agui  librium  distribution  Fj^  .  Tf  is  continu¬ 
ous  at  « o  this  point  will  not  contribute  to  the  expansion  formul  i 
(3.36)  and  the  eigenfunctions  (3.9)  for  /\M<  -o  will  not  be  contained  in 
the  expansion.  By  virtue  of  the  conservation  laws  the  functions  (3.9) 
will  then  be  orthogonal  to  (3.36).  Thus  (3.36)  is  complete  only  if  f 
has  a  jump  at  0  for  X~°)  I  .  It  follows  that  the  exclusion  of 
(3.9)  from  (3.36)  yields  an  expansion  which  is  complete  with  respect  to 
all  square  integrable  perturbations  conserving  and  9  . 


5.  The  Electron  Kinetic  Equation 


We  have  developed  an  expansion  theor.  uased  on  the  uncoupled 
kinetic  equation  (3.3)  for  the  ions.  The  extension  to  the  electron 
kinetic  equation  is  straightforward  and  requires  only  a  little  algebra. 
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As  we  indicated  earlier  we  decouple  the  electron  kinetic 
equation  from  the  ion  equation  by  dropping  the  term  (1*4 
in  (3.2),  This  amounts  to  neglecting  the  effect  of  the  ion  perturbation 
or  the  electron  perturbation,  but  retains  the  effects  of  encounters  with 
ions  in  the  thermal  distribution. 

The  uncoupled  equation  conserves  electron  number  density  provided 
(3.4)  holds,  but  does  not  conserve  momentum  or  kinetic  energy  in  the 
electron  gai.  This  is  as  it  should  be,  since  a  substantial  portion  of 

the  electron  momentum,  and  a  small  amount  of  the  energy,  is  lost  to  the 

ions. 

Applying  the  methods  of  section  1  we  find  as  before  , 

.  and  for  A  -O  we  find  const.,  corresponding  to  (3.9). 

The  remainder  of  the  development  proceeds  as  before.  A  spherical  harmon¬ 
ic  expansion  yields  a  set  of  singular  integrodifferential  equations,  and 
the  transformation  (3.19)  brings  these  into  self-adjoint  form.  As 
before,  the  expansion  formula  has  the  form  of  a  generalized  Fourier 
integral . 

If,  for  example,  the  ions  are  protons,  then  we  can  take 
If  we  use  *(e  in  place  of  e( •  in  the  definition  of  "V  and  c.  t  then  the 

electron  equations  may  be  obtained  from  the  ion  equations  by  replacing 

7/0  in  (3.15)  by  T(c)t  ^ZT (<~Uxc)  where  *  «  *t/<{  *  **«/*£  . 

6.  Discussion 

We  ha.tf  used  boundary  conditions  obtained  by  requiring  the  solu¬ 
tions  of  the  kinetic  equation  (3.3)  to  be  consistent  with  the  conserva¬ 
tion  laws.  The  Hilbert  space  then  emerged  as  a  natural  function  space 
for  the  framework  of  the  mathematical  development.  The  question 
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g 

persists  (see  e.g.,  the  discussion  of  Uhlenbeck  and  Ford  )  as  to  whether 
square  integrabi lity  should  be  a  requirement  on  the  distribution 
functions  from  the  beginning.  In  the  light  of  the  present  work  this 
condition  does  not  appear  to  be  necessary,  and  for  our  purposes  it  would 
not  have  been  sufficient.  To  see  this  we  note  the  condition 
JolV  °°  leads  to 


'b**-  ^-3/2  =  o 

ir->o 


(3.37) 


which  is  weaker  than  the  corresponding  condition  (3.4).  Since  both 
solutions  of  (3.20)  satisfy  (3.37)  for  O  ,  it  would  be  possible  to 
have  an  expansion  theorem  for  solutions  of  the  kinetic  equation  which 
are  square  integrabi e  but  do  not  satisfy  the  conservations  laws. 


IV.  THE  SCATTERING  OF  PHOTONS  FROM  A  PARTIALLY  IONIZED  GAS 
1.  Some  General  Properties  of  the  Scattering  Function 

In  the  first  part  of  this  chapter  we  present  a  brief  review 
of  a  classical  derivation  of  the  scattering  function,  and  discuss 
certain  properties  of  an  equilibrium  gas  relevant  to  the  computation 
of  the  scattering  function.  The  rest  of  the  chapter  treats  photon 
scattering  from  a  partially  ionized  gas. 

The  photon  scattering  can  be  characterized1  by  a  cross-section 
describing  the  effective  area  that  a  particle  in  the  sample  presents 
to  an  incident  photon,  having  direction  JX  and  energy  tco  ,  for  the 
scattering  of  that  photon  into  a  small  solid  angle  about  the  direction 
and  into  a  small  energy  increment  about  to)'.  It  can  be  shown 
that,  neglecting  relativistic  and  dispersion  effects,  the  electron 
cross-section  is  given  by1 

cr  fajQ]  ~  SfejAto)  (41) 

where  (Tr  (&)  is  the  Thomson  cross-section  and  £  ,a  and  Aco  are 
given  by 

K  -»  k  -  i  1  Ci&Q  =  k'  {t*  k'l  &  to  -  co — oj/ 

with  co-=-  c  k  ,  ooLck',  and 

iM*  ~ -2: ,  =-&'- 
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‘  I'i  I  (If  . ' '‘’Ilf? 1 
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The  so-called  scattering  function  S  is  given  by 

S/V^)  -  Cit e-iU" 6* (s, i)  '4,2) 

.  aO 

where 

f  e  •  (4.3) 

The  normalization  is  such  that  N  is  the  total  number  of  electrons  in 
the  scattering  volume.  It  can  be  shown^  that  the  scattering  cross- 
section  for  ions  with  mass  ***  is  of  order  (*e smaller  than  the 
electron  cross-section  (4.1),  where  is  the  electron  mass.  We  will 
thus  neglect  the  photon  scattering  from  ions,  assuming  local  charge 
neutrality  in  the  scattering  system.  Scattering  from  neutrals  will  be 
considered  later. 

The  function  |9C  in  (4.3)  is  the  Fourier-transformed  electron 
density  operator.  Since  the  ierivation  of  (4.1)  was  necessarily  quantum 
mechanical,  it  follows  that  the  density  operators  should  be  described 
quantum-mechanically.  It  was  argued  in  reference  1  that  the  difference 
between  the  quantum  and  classical  descriptions  of  the  density  operators 
will  often  have  negligible  quantitative  significance.  On  this  basis 
the  somewhat  simpler  classical  description  was  employed.  We  will 
continue  to  assume  the  validity  of  this  approximation  here.  The  reader 
interested  in  a  quantum  mechanical  description  of  the  scattering  func- 
tion  would  do  well  to  consult  the  work  of  Rosenbaum,  Zweifel,  et  al.  J 
We  are  clearly  concerned  with  the  electron  density  operators 
yo  €^,t)  , where 


(4.4) 
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and 

=  2.  i(t- Ur-ir*/^) ,  (4.5) 

We  seek  ultimately  the  thermal  average  of  the  product  of  the  Fourier- 
transformed  electron  density  operators  as  displayed  in  (4.3);  i.e., 
with  (4,4)  we  wish  to  obtain 

To  compute  the  thermal  average  above  we  foil ov/  Osborn^  and  generate 
a  set  of  equations  for  the  phase-space  density  operators,  which  we  then 
solve  subject  to  certain  well-defined  approximations.  Since  the 
procedure  for  generating  these  equations  has  been  delineated  elsewhere,^ 
we  present  only  a  brief  summary  here. 

Assuming  the  dynamical  variables  of  the  system  obey  the  classi¬ 
cal  equations  of  motion,  we  have 


where 

H  =  T+v 


is  the  plasma  Hamiltonian.  The  symbol  {  }  means  Poisson  bracket,  and 
for  any  function  A  of  the  system  dynamical  variables, 


T  Ud  _  'M  Ml 

4  L>*?apf  >f> t  >*tJ  . 
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It  thus  follows  that 


and  hence 


h  *  - 

3  -  la 


Alt)*  e'tLA(°) 


In  particular,  the  phase-space  density  operator  for  particles  of 
kind  A  is  given  by 


n 

=  e ' H  Z  {(Wto) 


(4.6) 


If  the  plasma  Hamiltonian  is  taken  to  be 

«-ii  ,  ig£*iZf , 

A**  A/* 


f  £  2  vM (!».*-&) 


a,b  *,p 

a*  a 


(4.7) 


then  it  is  a  straightforward  matter  to  show  that  ,i)  satisfies 

the  equation 
3 a  *  3a  /I 

-ft 

-  i  ^  •  ,i  Z  (4.8) 


Equation  (4,8)  is  similar  to  eqn.  (11.19)  in  reference  1,  but 
is  now  generalized  to  include  any  number  of  species  in  the  scattering 
system.  Now  as  in  reference  1  we  let  the  average  of  4  be  F*  f 
i  .e. , 
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FVs.y.t}  =  <3*^0^,  :«.9> 

and  further  define  the  fluctuation  operator  Scj*  , 

gVf.Jr^-F *&£,*).  (4.10) 


Still  proceeding  as  in  reference  1  we  combine  (4.8),  (4,9)  and  (4.10) 
to  obtain  an  equation  for  the  fluctuation  operators  $*<j  *,  and  then 
approximate  this  equation  by  neglecting  terms  quadratic  in  the  fluctua¬ 
tion  operators.  We  obtain 

+* •  w '  k  k  •  h  2  ( 

'  £■  ft  ' 


fr) 


f  =0.  (4#11 


Equation  (4.11)  is  now  further  simp!., 'led  by  assuming  that  the  target 
plasma  is  in  the  thermodynamic  state,  and  furth  that  the  singlet 
densities  r  ,  F  are  independent  of  space  and  time,  and  are 
Maxwellian  functions  of  the  velocity.  Euqation  (4.11)  now  reduces  to 


Hi-V.  Ik* .  j.  \5? .  x  M 

}  t  ~  ^  ax 


Z  (dVc/v  v**fo  •&)  ^  =  o.  (4J2) 


It  follows  from  (4.2)  and  (4.3)  that  we  must  solve  the  system 
of  equations  (4.12)  for  the  Fourier-transformed  electron  fluctuation 
operators  as  functions  of  time  for  all  t  ,  -  <*>  <  t<Q<3.  To  this  end 
we  introduce  respective  Laplace  transformations^  for  t>o  and  t<o  ; 
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OO  L_ 


O'tloo 


-  lir  Wpept^?4jr,  p) 


ps  (F  1 1\> 
(T  >  OV 


cr-i«o 


^<o) 


(4.13) 


'  A>V  .  /  _ 


-*»  /Tfi'ao 


Sg  2  =  -iw  i e,p*prSg- 

ft<6) 

=■  o  f*>») 


p  =  CTf  l  V 
(T  <71. 


(4.14) 


where 


S$Afe,ir,t)~  Sg*  + 

and  further  introduce  the  Fourier  transformation 

r.  A 


(4.15) 


fe,jr#p)»  (4.i6) 


We  thereby  obtain  from  (4.12) 


w2  v*8fc)  Vi--'-  p)  =  *y/w).(4.17) 


Vie  now  divide  (4.17)  by  (  p  -  tK-«r  ),  integrate  over  jr  ,  and 


define 
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and 


Kj  (&,?)  3  (e,e,t) 


a  ,  .  (jvsRaVsjt 

*6  ^  "  J  p-f« -r  ; 


obtaining 


(4.18) 


(4.19) 


1  am  ^  -  ±  ( jv  •  <4-2o) 


Now  let  A  represent  the  determinant  implicit  in  the  system  of  equations 
(4.20)  and  let  be  the  cofactor  of  h*  .  Then  solving  for  we 
find 


± 


f,efejr °) 

P  -£*•£ 


(4.21) 


Finally,  we  multiply  this  last  expression  by  ,  thermal 

average  the  product*  and  integrate  over  \r*  .  These  operation-'  yield 


pWV,<Sj‘Ytljr,Jo)^v./jrijf)Pj>_  =  (4-22) 


where 


<r6<Y°)  s  <  f 


(4.23) 


We  now  perform  the  inverse  Laplace  transformations  on  (4.22), 
as  per  (4.13)  and  (4.14),  and  employ  the  Laplace  convolution  relations,^ 
obtaining 
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<T  +  t'«  +■ 

±J_f  Jill,  dufep) 

ini)  f-ia-i  &(*,?) 

<T  -i  oo 


(  dr 

(Vxr^  t>o) 


—  OO 


.  oo 


dr 

(r«r.ti:«=)  (4.24) 


-  oo 


where  we  have  defined 


T-‘.oo 


t  4  ><r+ 

-  (4.25) 


Writing 

p6  =  d+®  a® 

we  can  easily  show  (provided  the  interparticle  potentials  depend  only 
on  the  magnitude  of  the  separation),  since  and  A  are  functions  of 
the  only,  that  D®*  obeys  the  symmetry  relations 

DVs(!-)  =  P Vvtf  -  DV-£,r)*  (4.26) 

and  further 

pV?,!r)  =  P£Y;,r).  (4.27) 

Combining  (4.3),  (4.4),  (4.10),  and  (4,15)  we  have,  with 

(dVdvwVe‘-'^'-VM7rjrM7jf')  - 
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where  Ne  »  ne  x  scattering  volume* 

(ye(£tt)  =  (in)1*  he  Ne  U*)  + 


+  £  (4.28) 

e  ' 

Fourier  transforming  with  respect  to  the  time  variable  as  per  (4.2)  we 
find 

Ste,  =  fa)**'  Hi)  i  (&<*>)  +  -jt  X  GWVi/f-Jf-4.oV!) 


Ute'iK'£r  »*/&,*). 

-  oO 


(4.29) 


To  complete  the  description,  i.e.,  to  portray  a  given  experiment, 
we  must  specify  the  interparticle  potentials  and  the  quantity 
defined  in  (4.23).  The  complete  specification  of  this  latter  quantity 
requires  a  fairly  detailed  knowledge  of  the  scattering  system,  and  is 
generally  a  formidable  computational  ^ask.  Nevertheless  certain  general 
properties  of  (j^e(o )  are  readily  established  and  are  germane  to  a 
description  of  the  structure  of  S{£/  Aw) . 

Consider  the  thermal  average  of  the  product  of  the  time-indepen- 
dent  density  operators  for  the  species  (  ,  8  ).  With  (4.5)  or 

(4.10)  this  quantity  may  be  written  either  as 

N  A  Ne 

<a  -  <X  7  Z(q-q*)Uq'-q*)\  (4.30) 

or  as 

<  3  Vs)  3  6  (Q')\  *  F  Vo)  F  V®')  f 

+  <y  (ti)  f9  & 


(4.31) 
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FA  r> 

_  ^  ,  F'  are 

space- and- lime  independent  Maxwellian  functions  of  velocity  rM  , 

_  6 

.  For  later  convenience  we  denote 

GAt(Q^')  =.  .  <4-32> 

Since 

[  jJ®cJ3©'<§  Y®)g  =  A/Ya  (4.33) 

by  definition  of  the  density  operators,  it  follows  that 

j  JiQ</sQ'<J',YG’,£?/)  -  O.  <4-34) 

Our  interest  here  is  in  the  functions  the  Fourier  transform 

being  needed  for  (4.29). 

We  begin  our  analysis  by  separating  (4.30)  into  two  terms; 

ifSMS'K  fifa-B*i>r  t 

+  .  (4.35) 

%A  ^.-8 

We  take  the  system  Hamiltonian  H  *  T+V  to  be  as  given  in  (4.7),  and 
define 

2r  s  pV*e-T/s  >  ?„  =  JdV'e-''7® 

where  the  respective  integrations  run  over  the  coordinates  of  all  N 
particles  in  the  system.  We  now  write  the  second  term  in  (4.35)  as 
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N*,N 


<  l  $  {q-q*)  l(G'-sie)')>  » 

A-6  r 


=  (dV^Ve'^'A 
'  '^s0 


VA,bll 


=  mY*)mV)Z  (jV  v'e~v/sJ6'^)MW) 


=  MVf)wV!r')  ^(VO 


(4.36) 


where 


n  *  M  Y*)  =  fY  rtf) 


A  4-  U 

and  h  is  the  number  density  of  the  A  species. 

To  determine  the  functions  defined  in  (4.36)  we  take  the 

gradient  of  (4.36)  with  respect  to  X  . 


>X 


=  2  W  ^  w  ?»"'  e  '  V/®  J  -x p) 


(/  /UB 


N>fe 


(jV  Vi-'ik-tfi (&'-&>)  e  V/e^  •  (4.37) 
e  ^  '  / 


From  (4.7)  we  have 


<r  vi<: 


v  =  i2  Z  v^Oi'-rO'I  1  v^o^-n) 

c  hr  ' ~  ■  ~ 


c%V'" 
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N ' 


and  hence 

fft-S  II 

d  ot,  T) 


0V<*  \j  °~A 

,  x*  ‘  s  v*V!3r-.n) 


+  T 


z 


ax' 


(4.38) 


In  the  following  we  will  neglect  the  second  term  on  the  right  in  (4.38). 
This  term  is  proportional  to  the  force  exerted  on  a  particle  by  that 
same  particle,  and  we  expect  that  the  effects  of  neglecting  this  term 
will  not  be  manifest  in  any  observable  results. 

The  relation  (4.37)  now  becomes 


drv/V*,*7) 


3x 


na,  N*>.  id  D 


- 'i  I  L  J  j  j *x M  2V- 1 


D 

^  c<  i/  8  -  ^ 

(r^-  c<y  D  -  a 


-ill 

D 

/!>  ijj  8-8 

(T^*  |V  D 


?,-* If*-*-) i <v-*p)  jjiW6"->r)«J 

-V/0 


(*)e 


We  now  separate  this  last  into  two  terms;  one  for  <r  »  (  and  hence 

8  -  D  )  and  one  for  <r  ^  |3  : 


n/*  k/® 

”©  Z  W  V d3x " 2,'1 1 Whs -2“) 

(i  ^  ^  8  =  d 


_V/6>  W *%-£'>!) 
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-  x  2T  )  A^^)^*)I4V)«(4i39) 

D  «,A,r  ~  _,  _ v/(9 

^o(  ^6=A,  <t**p<  jj  O-A,  rf£  (*)  By'"  £ 

Recalling  tne  definition  (4.36)  of  (£.£')  and  introducing 


V/© 

nsAWft,2',«*)s2.  \  J^N  ?v_l  e*  ; 

Vir 

/b/rf  ^6-A|<ry*  j/  D^A/O*^ 
the  relation  (4.39)  h'romes 

3*. 


.  I  46  *-*'!) 

J_  *©  "»  - 55 —  + 


a 


£  (j  V  „  *«Y8)  j  /jr »)  .  0  .  (4.40) 

r%  /  *%*• 


AB 


To  compute  the  two-particle  correlation  functions  ,  we 


,  *  ABD  , 

wn  te  r>^  as 

ABD 


^  =  rAV  t  ri'n^Vi <',*')  +-  ne  m/Vs ,  £*)  + 


4«  40 


.r,  W,r)  +  C'ft.t'.B*) 

where  kj*^  is  a  three  particle  correlation.  Inserting  (4,41) 
(4.40)  and  neglecting  k,  then  yields  the  system  or  equations 

3th',6fe,S')  .  I  ,  46 

- r -  +■  ft  -  t 

3  2.  9  32 

*  ft'.*-).*. 


(4.41) 


in 


(4.42) 
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Thus  given  the  system  potentials,  the  correlations  may  in  principle  be 
computed  from  (4.42). 

A  fe 

Two  pertinent  properties  of  the  y\x  can  be  demonstrated 
without  specifying  the  potentials.  To  obtain  the  first  of  these  we 
introduce  a  cl  inge  of  variables  in  (4.42)  according  to 


Equation  (4.43)  is  clearly  invariant  under  the  transformation  r  r 
If  we  now  add  a  constant  vector  to  r  such  that 

Ii+al  =  in, 

A  6 

we  find  that  (4.43)  is  unchanged.  It  follows  that  fix  is  a  function 
of  |r|  alone;  i  ,e. , 

n i *  rH*Vf£i)  ~  *\Ab6 £-*'/).  (4?44) 

In  the  segue!  it  will  be  necessary  to  have  on  hand  information 
regarding  the  normalization  of  .  Since  the  normalization  is 

already  specified  by  (4.33),  we  merely  combine  (4.30),  (4.31),  (4.35), 
and  (4.36)  to  find  as  a  function  of  : 

+  .  (4.45) 


Integrating  now  over  (  X,  * ' ,  \r  ,  \r*  )  we  find,  with  (4.34), 


-  92  - 


or 

j  e/*V  ni*Vr)  *  6  ~  ^Ae)  . 


(4.46) 


We  return  now  to  our  discussion  of  the  scattering  function. 
Fourier  transforming  ,  we  have 

where  (r^/o)  ir  given  in  (4.23).  With  (4.44)  it  follows  that 


G-8efe,iri)f,/«0  =  0-^(-£.^,£\o) 


e>e 


(4.47) 


A& 

and  since  H  is  real  we  also  have 


* )-  fr*7wvo. 


Se*/ 


(4.48) 


Combining  (4.29  and  (4.47)  gives 


S  (hl,aco)=  S(-&,- AcS)  ■  (4-49) 

It  follows  easily  with  (4.27),  (4.29)  and  (4.48)  that  S(£, 4 u>)  is  real, 

5 

as  we  would  expect.  W.th  (4.28)  we  also  have 


Ge*(i  ,t)  =  G-ee(*-t)=  &"*(*, t).  (4.50) 

Combining  (4.2)  and  (4.50)  we  find 

S'  00  j. 

<lte-itCC*  (4.51) 

o 

Hence  when  (r$e satisfies  (4.50)  we  can  compute  the  scattering  function 
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~ee  . 

having  o.  ly  w'  for  V>  o  . 

At  this  point,  instead  of  inverting  the  Laplace  transform  as  in 
(4.24)  and  (4.25),  an  operation  straight- forward  in  principle  but 
generally  Herculean  in  practice,  we  adopt  a  procedure  due  to  Rostoker.6 
Thus  we  consider  first  the  identity 


_L  f  [Jt* 


*■  die 


t(  v'^)t 


5^ 


5  $+  +  S_ 


where 


.9  - 


3-%^  ti;p|  tfer  S(V0 


(4.52) 


(4.53) 


since 


riUK^)t  ft'p^ 

'o 

^Jt  ti(vLv)t  .  -i  P 


(4.54) 


with  P  indicating  principal  value. 

Now  let  S  +(&)p),  f-r+iv  be  the  Laplace  transform  of  some 
function,  ar j  consider  further  the  following  inverse  transformation 


t 


0~t  too 


r->o+ 


2W  S* 


(T-l®0 
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(T-H  *0 


J  c/pe pt  S+te,fO  « 


5“  -  t  « 


yoo  ,co 


=  ™  j  \jt'et(','~*)t Sfe,*') 

{.  J  04  *  Q 

=  £■  (!>'  Sfc«')  Cdt'e1'1'*1  Mi-t') 


*  eivtt(*,») ,  t  >  o  ;  -  «,  t«®.  (4‘55) 

_  OO 

In  arriving  at  (4.55)  we  have  used  the  definition  (4.53)  of  the  function 
*  Substituting  (4.2)  into  ( ^ , 55 )  we  have 

tZ*  w  (  dp«ftS,fe>>  U'  j*e"N)^« 

<r-i*  ®o  -  oo  _  oo 


'  .^je  cfW),  so,  (4-56) 

It  follows  that  and  the  Laplace- transformed  correlation  function 
for  t>  O  are  related  via 

tz*  %(*,*«*)  -  r^-vl,'£  (4.57) 

From  (4,53)  and  the  knowledge  that  the  scattering  function  is  real  we 
have,  with  (4.51^ 


AcS)  „  J—  &  si*  . 

° 1  <rs>ot  1  1 


We  now  combine  this  last  with  (4.3),  (4.4),  (4.10),  and  (4.22).  After 
Laolace  transforming  the  forward  scattering  term  S(n) 


and  noting  that 
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jQ. 

7TA/e  ^ 


^  (ir^  ^  N^(a) 


Q-  f  l  LCO 


(i*)*  **  <£(*)&  [u*<v)+  jf  ? 'to'] 


we  obtain 

S(x,Lco)  -  (tn )*  r\e  $(*)&(  + 


-F 


J3sJ  s' 

(7^t  i  Leo  -  tV  •  tr 


(4.58) 


2.  The  Electron  Scattering  Function  For  A  Partially  Ionized  Gas 

Ii  this  section  we  derive  the  electron  scattering  function  for 
a  system  competed  of  electrons  and  one  species,  respectively,  of  posi¬ 
tive  ions  and  neutral  atoms.  The  extension  to  a  more  general  multicom¬ 
ponent  system  is  straightforward  but  adds  considerably  to  the  algebraic 
complexity.  We  will  continue  to  assume  that  the  equilibrium  plasma  is 
characterized  by  a  single  temperature  conmon  to  each  species,  and  is 
free  from  spatial  gradients  or  external  fields.  Ue  present  two  differ¬ 
ent  treatments,  with  somewhat  different  results. 


A.  Reversible  Theory 


From  (4.58)  we  have,  ’.aglecting  forward  scattering, 

S(V^)=  -L  fit  -L(*> 


(*) 


Jl  A* 

-  A (l  +  -A,*-  A,,;  )(J ‘ej^  !  - 


A  -  (-^ca  / Aei )^/f~  Ag  )  ~l  Ga  1 
I  +  c.)  j  ltJ 


(4.59) 
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where 

A  —  I  -A.^g  ( I  "f  -A^  _A^t-Aj'^  ^  A 


i  A u(l*A„  -A*  A,„  -A^^f  -  A»j*) 


An 


■Au 


and 


5  )^<f3e^0)^O)>r- 

In  writing  (4.59)  we  have  made  use  of  the  identities 

•A  d  Aj*e  =  Atf 

-A+p  Ap  Afa  -a  -A/&,  Aty  A^. 


We  have  displayed  the  scattering  function  in  (4.59)  in  a  form  that  will 
facilitate  an  estimate  of  the  significance  of  the  terms  involving  the 
neutrals.  To  this  end  it  will  prove  useful  to  write  .A^in  the  form 


(T-»o  + 


A 


’  AB 


Aft 


9 


[■- 


oC 


Afo  - i/i  (due 


—  u 


71 


(4.60) 


u  +■ 

r\~  °° 


where 


*  ('ze/W)'^ 


and  the  sumbol  "/) "  next  to  the  integral  means  the  path  of  integration 
is  deformed  above  the  singularity. 

To  estimate  V  A  8fe )  when  either  or  both  of  the  pair  (  A  ,  6  ) 
is  a  neutral  atom,  we  assume  the  potential  may  be  approximated  by  a 
Yukawa  potential,  VA^(r)  - 
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where  <l  is  the  effective  range  of  the  potential, 

CL  *  I  0~B 

and  £ ^  is  a  constant  to  be  determined.  We  have  easily 

VA6/£)  =  Vucm  • 

With  -X  >>  lo~b  c~*.  ,  kx&7«  I  and  thus 

VA‘(S)*«*  <«•«!> 

To  estimate  the  constant  we  recall  that  the  center-of-mas^  differ¬ 
ential  elastic  scattering  cross-section  for  the  pair  (  A  ,  6  )  having 
relative  momentum  t  H  is  given  by7 

where  is  the  relative  mass.  From  (4.61)  it  Follows  that,  in  the 

energy  range  of  interest,  is  approximately  isotropic.  Writing 


(T*6  =  ^J(lO~At(e>)  «  lo~U  (T0  (c^1) 


8 

where  <r0  is  of  order  1-IOc*^  ,  we  have 


'A  6 


To  estimate  the  various  terms  in  (4.59)  we  assume  rv\***v  and 


take  to  be  the  mass  of  the  C'1,  atom,  ^  2.  5/  le~*  . 


.-25, 
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Assuming  for  simplicity  equal  electron  and  ion  densities,  we  have,  for 

X  ^  A  . 

Av 


A,  L_  ^ 


M 


(W*1)  *' 


A-.e  ^  _ 

Ac  ^  *>e  X.' 


Ac* 


*  2  *  to- "  £  rj/x6~)  ) 

v  •  (4-62) 


■ft 

Here  ©  is  the  laboratory  scattering  angle, 

UI--  li -i‘l  *  hx-1  | 

and  appears  only  in  the  argument  of  the  sine,  in  contrast  with  the 
kinetic  temperature,  also  denoted  by  9  . 

n 

Noting  that  the  quantity  in  square  brackets  in  (4.60)  is  of 
order  one  or  less  for  all  values  of  ^«o/k<£  ,  we  find  (for/\-  69*12  X 

/K. 


and 


A  n,  | 

L-AtJ 

1 

r»  .  n 

Mo, 

Ap*  | 

<T->o  +  j 

Mo,  >#  Aph 

'  -A* 

i  /  —  "If-  /  %  \  V~o 

*g  l.  6x10  (*~<  -*^r)  - r-  I  * 

Q(Jat) 


* 

V\C  I 


<ro]^[cA 
Q(jur)  L 


(4.63) 


@Y* 


; . 

V\  g  V*A  . 


It  is  clear  that  the  quantities  in  (4,62)  and  (4.63)  are  smaller  for 
larger  **h(or  )  and  are  largest  (about  a  factor  of  twelve  larger 
than  the  abofe  values)  for  a  hvdrogeneous  scattering  system. 
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To  determine  the  significance  of  electron-neutral  correlations 
in  the  scattering  function,  we  must  estimate  the  magnitude  of  the 
quantity  (see  (4, 45)  and  (4.59)  ) 


or 

6* 


[W O-  WVj 

jc/Ve4*  "L  [ n14f  (r)  -  v\‘ v>fJ 


(4.64) 


where  we  have  assumed  .  To  facilitate  computation  of  the 

pair  correlation  functions  Y\\  appearing  in  (4.64)  we  introduce  the 
assumption  that  the  average  distance  between  any  pair  of  particles  is 
large  compared  with  the  effective  range  of  the  charge -neutral  or 
neutral -neutral  potentials.  This  range  being  typically  of  the  order  of 
I  O'  the  assumption  implies  particle  densities  small  compared  with 
IOlH  .  Under  this  assumption  it  follows  that  the  contribution 

of  neutrals  to  the  correlation  between  charged  particles  may  be 
neglected.  We  thereby  obtain  from  (4.43)  a  pair  of  equations  for 

(r)  .  Taking  advantage  of  (4.44)  we  can  write  these 


i  >"*  4'  i  Kf?,  - 


<e(r ) 


where  is  the  ionic  charge  and  r\e-^,-wc.  The  equations  above  are 
similar  to  a  pair  of  equations  treated  by  Lamb^  for  a  singly  ionized 
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gas;  with  ^  -  I  Lamb's  equations  coincide  with  ours.  Following  La^b, 
we  first  perform  the  angular  integrations  and  then  differentiate  with 
respect  to  r  .  The  result  is  a  pair  of  differential  equations 


i'V.O  (4.65) 


t  "w 


The  system  (4.65)  permits  a  non-trivial  constant  solution. 

Taking 


A 


»•  > 


we  find 


A j  *■  , 


Lamb10  has  shown  that  for  distances  r  large  compared  with  e1* / &  , 

the  first  derivatives  in  (4.65)  contribute  negligibly  to  the  solutions. 
Neglecting  these  terms,  it  is  then  a  simple  matter  to  show  by  direct 
substitution  that  the  functions 

Pee  =  li.-r/z'#  f‘e  =  -'■/A* 

V  r  e  )  +  r 

satisfy  (4.65)  provided  8,  =  -  6-j. ,  and  where 

\  -1  '  t^t‘) 

D  9 

We  thus  take  the  solutions  of  (4.65)  to  be 

nf  (0  -  -  &Jbe-r/Ao  f  ^ 


r>>  Hi 
e 
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The  constants  Alt  6^.  are  determined  easily  with  the  normalization 
conditions  (4.46),  so  that  finally  we  have 

['  fr  +  I  - 


r\lr\c 


[»£  . 


(4.66) 


where  A/  =  N  e  f-  N L  . 

To  estimate  the  significance  of  electron-neutral  correlations 
as  per  (4.64),  we  write  r\*e  <^r)  in  the  form 

y\?fr)  =  [i i  *&)  H  '  jj] 

A 

where  7  is  a  function  that  we  expect  will  differ  appreciably  from 
zero  only  for  f  ^  a.  **  (o  0*.  &  of  the  electron-neutral 

potential.  The  normalization  condition  (4.46)  gives 

YO  *  jj  -  . 

Vv  6  Cg 

Substituting  the  abrve  form  of  r\v  together  with  from  (4.66) 

into  (4.64)  now  give:,  neglecting  the  contribution  from  terms  in  \/fsi 
(which  corresponds  to  the  neglect  of  forward  scattering) 

rh- 

With  icr5’^')"1  where  A  is  the  photon  wavelength,  we 

approximate  the  integral  above; 
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pVe'S-i  4Mfr)  =  ( p*(r) 

<K^  (o~a  e^, 

*  (jv  q*efr) 

ItK^ 

/V  ^**  +  *c)H  . 


Thus  finally  we  have,  with  A  =  J-X  (o~*  , 


tf;-0. 


We  a*^  now  in  a  position  to  estimate  the  magnitudes  of  the 
terms  involving  neutrals  in  (4.59).  First,  we  note  from  (4.62)  that 
A**  /Aft  is  negligible  unless  **/*e  *  lo,z  or  larger.  In  such 
cases  the  total  light  intensity  scattered  by  the  electrons  is  greatly 
exceeded  by  that  scattered  by  the  neutrals.  Since  our  interest  here  is 
in  the  influence  of  neutral  atoms  on  the  electron  scattering  function, 
and  since  we  do  not  anticipate  an  experiment  in  which  scattering  from 
electrons  could  be  observed  at  such  extreme  density  ratios,  we  will  not 
consider  these  extreme  cases  here. 

It  follows  from  (4.63)  that  (4.59)  may  be  further  simplified 


provided 


r.1^  ,  2? 

— - —  <<r  o 

9  (at) 


(4.67) 


Since  particle  densities  do  not  generally  exceed  /o*36~ra(e.g. , 

\o  <^~cVor  graphite;  \o,8v Gfar)  at  atmospheric  pressure), 
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we  would  expect  (4.66)  to  hold  for  almost  all  plasmas  of  interest  here. 
Thus  when  (4.62)  and  (4.63)  are  satisfied  (4.59)  reduces  to 

S  (£,  b'-o)  -- 


__L 

7 TSJe 


fia .  ^  . 


_L 

A 


(l  a)  Go 

<  - 


(4.68) 


where  now 

A  ^  I  4"  A  -f  ~  A  e*})  , 


The  dependence  of  the  scattering  function  displayed  in  (4.68) 
upon  electron-neutral  correlations  (the  term  containing  0^e  )  disap¬ 
pears  when 


<<  \o 


ll 


(4.69) 


If  in  addition  the  condition 


tc  (c~y)  /\w  .yj- 

— - - - - - - -  <<  |o 

Q%us) 


(4.70) 


holds,  then  Aev^yiKve«  /  and  (4.68)  reduces  to  the  result^  for 
a  fully  ionized  gas.  For  y\*  ^  iczl  ,  (4.70)  becomes 


G~o  v\e  (c — '•~5  3 
9x(*r) 


<.<  /<? 


(4.71) 


It  is  thus  apparent  that,  unless  very  high  neutral  and  electron  densi¬ 
ties  are  present,  together  with  relatively  low  temperatures,  the  present 
theory  does  not  predict  an  observable  effect  of  neutral  atoms  on  the 
electron  scattering  function. 
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Before  abandoning  our  quest  for  an  observable  effect  of  neutrals 
in  the  spectrum  of  electron-scattered  photons,  we  turn  to  a  somewhat 
different  and  more  realistic  formulation  of  the  scattering  function. 


B.  Irreversible  Theory 


It  was  noted  in  Chapter  II  above  that,  since  there  are  generally 
many  charged  particles  within  range  of  each  other  simultaneously  for 
most  plasmas,  the  effect  of  close  binary  encounters  may  often  be 
neglected.  In  this  case  one  may  represent  the  effects  of  particle  inter¬ 
actions  by  an  appropriate  electric  field  term  in  the  kinetic  equation. 

In  contrast,  we  would  not  expect  such  a  "field  representation"  to  be 
suitable  for  the  representation  of  encounters  between  particles  having 
ranges  of  interaction  that  are  small  compared  with  the  mean  interparticle 
distance.  In  the  following  discussion  we  adopt  a  scheme  outlined  in 
Chapter  II  designed  to  give  a  more  realistic  treatment  of  the  inter¬ 
actions  between  charges  and  neutrals  than  that  employed  above. 

We  begin  with  equation  (4,11)  for  the  fluctuation  operators  of 

XL 

the  A  species; 


* 


]_ 


.SL 

)\r 


(4.72) 


At  this  point  in  ref.  1  and  in  Part  1  above,  the  assumption  of  thermal 
equilibrium  was  introduced  for  the  target  plasma.  Before  doing  so  here 
it  is  convenient  to  exploit  the  difference  between  the  relatively 
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long- range  Coulomb  forces  between  charged  particles  and  the  relatively 
short-range  charged-neutral  and  neutral-neutral  forces.  To  this  end, 
and  with  an  eye  01  our  ultimate  goal,  we  multiply  (4.72)  by 
average  the  result,  '•nd  introduce 

(4-73> 

We  thus  obtain  from  (4.72) 


■jp* 

it 


+  vr-  _L  J  (jV'JV'  5V/>e('— 

L-.  I  A 


We  now  introduce  a  change  of  variables  according  to 


r  =  x  -  x ;  c  '  =  x/J-  x1 

W.  j  Wk  WV_  J 

so  (4.74)  takes  the  form 

.  J_  J  OwV' 

u  ~ 

■iJFVs-or^xit)  r  FYs;-;t)rys,s,t)7-c. 

In  writing  (4.75)  we  have  suppressed  the  dependence  on  the 
variable  x'  .  Our  identification  of  the  integrand  in 


(4  75) 


(4.76) 


as  the  equilibrium  phase-space  correlation  function  for  the  pair 
(A »  £  )  implies  that  P*e(°)  is  a  function  of  £  and  yf  only.  It 
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follows  from  (4,75)  that  P  (t)  is  a  function  of  r  ,  £  ,  and  t  alone. 
The  identification  of  (sAe(£l  i)  , 

CrAt  =  L’V  (4.77) 


as  a  time  dependent  correlation  function,  and  the  relationship  between 
the  scattering  function  and  the  space-time  Fourier  transform  of 
has  been  exploited  above.  The  description  of  r  by  a  linear  trans¬ 
port  equation,  as  in  (4.75),  has  been  suggested  on  the  basis  of  a  semi- 

intuitive  argument  by  Nelkin  and  Ghatax,^  and  has  been  employed  by 

■\  ■> 

them  and  by  Yip  and  Nelkin  L  in  a  study  of  slow  neutron  scattering 

13 

from  liquids  and  dense  gases.  Recently  Van  Leeuwen  and  Yip  have 

A  ft 

derived  a  similar  kinetic  equation  for  r  ,  for  short  range  potentials, 
from  the  cluster  expansion  of  a  one-particle  distribution  function. 

At  this  point  we  introduce  an  approximation  into  the  treatment 
of  the  charged-neutral  and  neutral -neutral  interactions.  This  leads  us 

ultimately  to  the  binary  collision  description  attained  by  Van  Leeuwen 

13  1112 

and  Yip,  and  employed  by  Nelkin,  et  al_.  *  We  thus  adopt  the 

treatment  outlined  in  Chanter  II  in  going  from  equation  (2.17)  to  the 

Boltzmann  collision  integral  (2.24).  Thus  the  terms  in  (4.75)  involving 

the  relatively  short-range  neutral  interaction  potentials  V*"/l 

are  approximated  as  in  Chapter  II  by  a  linear  Boltzmann  collision 


integral,  and  (4.75)  becomes 


an? ,  i£ie  -  i  y  Lw  . 

tJ  ^ 
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=  jc/V"  (j/i  |vr-jr*|irM  [FVf,*,*)r','Vr,sr.")i)  + 

+  F7r, /.",£)  n'7x,/.,t)-  FVr,^«r'7r,r;t)  - 

-  F7r,or*t)  ret/r^,t)7  (4-;8) 

) 

nl*ft 

and  similarly  for  1  .  At  this  point  we  introduce  the  assumption 

that  the  target  plasma  is  in  the  thermodynamic  state,  and  hence  that 

the  singlet  densities  F  *  ,  a^e  independent  of  space  and  time  and 

are  Maxwellian  functions  of  velocity. 

We  have  argued  in  the  preceding  section  that  electron-neutral 

correlations  should,  under  most  conditions,  contribute  negligibly  to 

14 

the  electron  scattering  function.  It  is  worth  noting  that  Salpeter 
has  argued  semi -intuitively  the.  pair  separations  which  are  small 
compared  with  both  the  photon  wavelength  and  the  Debye  lergth  contrib¬ 
ute  negligibly  to  the  spectrum  of  photons  scattered  by  electron 

density  fluctuations.  This  is  in  good  agreement  with  the  experimental 

1 5 

observations  of  Ramsden  anH  Davies.  This  suggests  that,  since  the 
electron-neutral  correlation  is  significant  only  for  separations  of 
order  )0 ~  or  less,  it  would  be  reasonable  to  ignore  electron-neutral 

correlations  in  our  computation  of  the  scattering  function.  A  somewhat 
different  argument  in  support  of  this  assumption  is  suggested  by  the 
observation  that  fluctuations  induced  in  the  neutral  distribution  by 
recoiling  electrons  should  be  insignificant  Referring  to  (4.73)  we 
see  th-_  neglecting  neutral  fluctuations  implies  neglecting  electron- 
neutral  correlations.  With  these  assumptions  (4.78)  becomes 
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n  +  -  '  ar 


J  jjvjvvft-ri)  rK(s',r,t)- 

h*etL  ' 


v\ 


^JV(lcLU-a:"l%r, 


and  similarly  for  P*"  .  It  is  apparent  that  the  neutral  species  colli¬ 
sion  operator  in  (4.79)  conserves  both  electron  and  ion  number  densities, 
but  does  not  conserve  total  momentum  and  kinetic  energy  in  the  system,  in 
contrast  with  the  collision  operator  in  (4.78)  which  satisfies  all  three 
conservation  laws.  This,  of  course,  results  from  our  neglect  of 
fluctuations  in  the  neutral  distribution. 

An  important  consequence  of  the  collision  approximation  manifest 
in  (4.78)  or  (4.79)  is  that  while  (4.75)  is  invariant  under  the  time 
reversal  transformation  t->~T  ,  jf ,  these  last  two  equations  are 
not  invariant  lie  will  return  to  consider  the  consequences  of  this 
irreversibility  shortly. 

We  now  introduce  a  further  approximation  into  the  collision 
description  in  order  to  avoid  the  complexities  of  the  collision  opera¬ 
tors  in  (4.79)  in  their  present  form.  We  thus  replace  these  operators 
by  the  linearized  version  of  the  simple  single  parameter  collision  model 
first  proposed  by  Chatnagar,  Gross,  and  Krook,^  and  discussed  In 
Chapter  II  above.  The  model  is  constructed  to  satisfy,  in  this  case, 
the  requirement  of  number  conservation  for  each  species,  preserves  the 
irreversible  nature  of  the  above  description,  and  provides  a  considerable 
simplification  for  the  subsequent  analysis.  Finally  then,  the  kinetic 
equations  that  we  use  are 
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Tne  parameters  y»'  ,  v'<  are  clearly  lectron-neutral  and  ion-neutral 


collision  frequencies.  Fourier-Laplace  transforminq  (4.80) 
as  before,  we  eventually  find,  with  (4.77) 

and  (4,81) 

G^(i,p)=  A-' (n-  A.y-v'-irr)  (  _ 

+  J  p- iT 

.-1. 

—  VV  ei  )  n  <  i 

/  p-  4  *  'if  r  v 

(4.8?) 

where 

)  p-iff'if+'i)**  } 

(4.83) 

A  »  (i  1 Ae  -  V  "ITT  )(/  +  A,v  -  vyr  )  -  AV  ^  «• , 

(4  84) 

and 


t/'Y?) 


& 


1-  {f>t» 


V\  A 


(4  85) 
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The  functions  P^^vC  o)  are  obtained  from  (4.76)  and  the  Fourier- 
transformed  equilibrium  phase-space  correlations.  Thus  with  (4.45)  and 
(4.66)  we  have,  neglecting  forward  scattering, 

nf7 >,£,<>)  =  H-lWj  7  (4.86) 

and 

P%,  f°)=  1  (4.8V) 

where  is  the  avarage  ionic  charge. 

As  a  result  of  our  having  employed  an  irreversible  theory  in 
the  description  of  the  correlation  function,  it  is  readily  shown  that 
this  quantity  is  not  symmetric  under  the  interchange  t in  con¬ 
trast  with  (4.28).  In  addition  G-ee(&,  t)  diverges  exponentially  for 
so  that  the  integral  in  (4.2)  does  not  exist.  To  overcome 

il  ip 

this  difficulty  we  follow  Nelkin,  et  al_. ,  *  and  prescribe  a  behavior 

for  negative  times  different  from  that  we  would  obtain  by  solving  the 
system  (4.78),  (  Pfe,  P'e)  for  t<  O  .  The  prescription  ensures  the 
convergence  of  (4.2),  yields  a  real  scattering  function,  as  it  must, 
and  the  result  is  symmetric  in  k  and  A  u>  as  required  for  classical 

5 

systems.’  The  prescription  is 

(srt)  =  (4.88) 

Fourier  transforming  (4,78)  for  both  specie?  we  can  show 

Gn(s,t)  «  &“(-£,+).  (4*89) 

With  the  symmetry  property  (4.83)  we  can  now  compute  the 
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scattering  function  as  before  using  only  the  Fourier-U;  lace  transformed 
electron  correlation  function  for  t  >  o  .  Thus  with  (4,59)  and 
(4.82)  we  eventually  obtain 

S  (stu)  =*-'  (k  to 


where 

h  ~ 


[‘  +  Sysktf  1'-^ 
1  * 


and 


^  '  0+f,-)[n-  (kaH'J 

It  will  prove  convenient  for  computational  purposes  to  write 
(4.90)  in  terms  of  dimensionless  variables.  We  thus  introduce  the  new 
variables 


A  - 


-I 


C  ^ 

i  - 


£  Cjo/ CO 


pe 


■*1*  2  V 


y« 


COc 


(4.91) 


uj 


pe 


v/h  ere 
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We  obtain  from  (4.83) 


cr  ->o+ 


-11*  *  -if.  (%)y  [”'*  (  fc  -  r*mft  vJ 


=  i. 

'  *W  ^  <2^ 


(4.92) 


where  the  complex  function 


z?  (x^'g) 

Q 

is  the  plasma  dispersion  function  tabulated  by  Fried  and  Conte  and 


The  scattering  function  now  takes  the  form 

S(f,A«)  -  fr«)J  (%)'a(V1h)-'  hi  cr'  M 


(4.93) 


-2i  'Gtfrz;%3 


(4.94) 


with 


£f AO  •  -z[t +t'2(t)])  out 
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The  scattering  function  displayed  in  (4.94)  will  differ  negli¬ 
gibly  from  the  result  for  a  fully  ionized  gas^  when  the  imaginary  parts 
of  the  arguments  of  the  functions  ft)  are  small  compared  to 

unity; 

«  i .  (4.9f) 

To  secure  an  estimate  of  when  this  condition  holds,  we  write  the 
collision  frequencies  V**  as 

*  (T/V**1)*"  C^Tx/o'1*  '4-96> 

A  2 

where  <ro  is  of  order  one  to  ten  cm  as  previously.  With  (4.91), 
(4.95)  becomes 


_ *  .  -  lb 

<T0  yT  x  lo 


«  I 


or 


f  *  ® 

"j. 


«  3 ,S*  lo 


i? 


(4.97) 


For  photons  emitted  by  a  ruby  laser  (  X  ~  69H3A  ),  this  can  be 


wri  tten  as 


CV*-1)  i'Y'"-*) 


£ 


(4.98) 


Comparing  (4.98)  with  (4.66)  and  (4.70)  we  find  that  the 
scattering  function  (4.94)  obtained  with  the  irreversible  theory  will 
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display  neutral  atom  effects  at  neutral  densities  at  least  five  or  six 
orders  of  magnitude  lower  than  those  necessary  to  see  these  effects  in 
the  reversible  result  (4.59).  In  the  following  section  we  describe  and 
discuss  the  photon  scattering  function  as  displayed  in  (4,94). 

3.  Photon  Scattering  From  a  Partially  Ionized  Gas 

In  this  section  we  discuss  the  results  of  the  previous  section 
for  various  plasma  configurations.  For  convenience  we  take  =  r**. 

To  obtain  quantitative  information  regarding  the  scattering 
function,  it  is  necessary  to  claculate  the  real  part  in  (4.9+) .  To  this 
end  we  write 

-  k  -n> ,  (4.99) 

and  define  )'fX  .  and 

=  ('-  (hi-i-X) 

bi  =  ( I  -  X;)  V  =  (^1  T') 

b3  s  (J'X  -  -fA'Hi )  ^  *  (i'h  4-  ke  ±;) 

(4.100) 

with  as  per  (4.93).  After  some  algebraic 
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manipulation  (4.94)  becomes 


S  (£,&*))-(**)  (l£) 


-i  /  wg_  ~  a  j- 


+■  A. 


(4.101) 


where 

Ne  *  +Ie<u)  ~ 

Nr  s  ftf{Jn+i)  ^F°'i  -I,*,)  +■ 

Ar  =  +  kfcj  <-  ^  .  (4.102) 

Using  tabulated  values^of  »  ^a  »  and  ^  we  have  com¬ 
puted  S  as  a  function  of  f  -  A«jA^e  for  different  values  of  V  and 
Je}^  •  Some  typical  results  are  plotted  in  Figures  4.1  through 
4.3.  To  enhance  our  understanding  of  these  results,  we  now  develop 
approximations  to  (4.103)  in  various  limits.  To  this  end  we  employ 
both  the  power  series  and  asymptotic  expansions  of  the  plasma  dispersion 
function.  These  are,  respectively,  with  , 

2Y t)  -  -  2t  JV  2*l/3  t  (4.103) 


and 


t7  /»  » 
a  /rr  1 

c_i  1-  / 


-  i  r  77  *'  ~t  -  i  j-i  -r  (  f/2t  )  t  104) 


where 


r  - 


0  1 1 
I  u  =  O 


(4,105) 
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Consider  first  S  as  a  function  of  /  .  For  I  collective 
effects  play  a  negligible  role  in  the  scattering  process.  With  (4.103), 
(4.101)  reduces  under  this  condition  to 

r— »  ■¥• 


*  (**)"*( lb)  hi 


i  -  e^Jt'  2f, 


* 

K 


;  O-YJfi 


(4.106) 


This  is  intical  in  form  with  the  result  of  Ghatak  and  Nelkin. 
(4,106)  reduces  to  the  well-known  ideal  gas  form 


11 


For 


^  —  o 


-  (0-),/lexp  [-”t(*“>MVze]  _  (4.10?) 


For*fe>0  ,  the  scattering  function  given  in  (4.106)  becomes 
narrower  and  increases  in  height  as  increases  in  magnitude.11 

For  I  ,  w\?«  /  ,  wi  and  ^ or 

less,  (4.101)  becomes 


%(*&«>)*  2  -  7  ^  ^  .  (4.ios) 


This  last  is  similar  to  (4,106)  but  now  the  scattering  electrons  are 
strongly  coupled  to  the  ions,  and  the  ion  mass  replaces  the  electron 
mass  in  the  scattering  function,  resulting  in  a  narrower  scattered 
photon  spectrum. 

Next  consider  3  with  the  condition 


«l)  ■Je  ~  Y/"f  - 


(4.109) 
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When  (4.109)  holds  we  find 

Sf*  A")  -  feW) 


W 


(4.110) 


In  this  case  the  center  of  the  scattered  photon  spectrum  is  dominated 
by  a  relatively  narrow  "ion  peak"  provided 

/  2  .  .  \l- 


or 


/**>  CiyA'  +>) 


(4.111) 


Under  the  appropriate  conditions  the  scattering  function  will 

have  a  resonance  at  the  electron  and/or  ion  plasma  frequencies.  Thus 

for/(  ^  3  ,^e« /  and  f'W  or  larger,  S  takes  the  approximate 

★ 

forms 

(&f 

(■f-o)  W  tl 


(*) 


}-(  ~ihT  ] 

Vpj 

(4,112) 


-  « 
a- 


(7e»»o) 


_/  (4.113) 


*  The  reader  will  note  that,  because  of  the  two  Afferent  asymptotic 
forms  (4.105),  the  form  (4.112)  cannot  be  obtained  from  (4.113)  with 

y+o . 
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The  resonance  at  is  apparent  in  (4.112)  and  (4.113).  It  is 

clear  from  (4,113)  that  the  effect  of  collisions  -5s  to  inhibit  the 
resonance  effect,  as  we  might  expect.  The  height  of  the  resonance  is 
inversely  proportional  to  while  the  width  increases  as 
increases . 

For  I  ,  ,  and  =  0  ,,  we  find 


(X) 


where 


exp  , 

[-hirZS] 

1 

( i 

[-W%T] 

(4.114) 


OJ 


£  *  jinslu 

rvv 


0 

and  is  the  electron  number  density  as  before.  For  we 

find  a  form  similar  to  (4.113). 

Finally,  consider  the  scattering  function  in  the  limit  of  very 
strong  collisions.  With  and  /  ,  S  takes 

the  approximate  form 

S  (&,**>)  * 

(**)"' (le)"1  (hi  - ^  (4.115) 

Hence  in  the  strong  collision  limit  collective  phenomena 
are  unimportant  (except  possibly  for  wa  f  ^  I  ;  see  (4.110)  ),  as 
we  would  expect.  The  effects  of  electron-neutral  collisions  dominate, 
and  the  approximate  form  (4.115)  of  the  scattering  function  is  similar 
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to  (4.106). 

In  Figure  4.1  we  have  plotted  the  normalized  scattering  function, 
for  a  single  ionized  (  •  =  /  )  gas  of  carbon  atoms,  vs.  ?  as  a 

function  of  the  parameter  X  for  Y  9  Y  (zero  effective  neutral 
density).  The  general  qualitative  dependence  of  $  on  the  value  of  4 
is  clearly  apparent.  forA<<!  ,  &  has  the  Gaussian  form  (4.107). 
For^t  /  the  resonance  at  f  y  I  is  present,  becoming  narrower  and 
higher  as  J  increases,  as  per  (4.112).  When  f  is  small  the  effect  of 
ions  dominates  in  S  as  suggested  by  (4.110)  while  an  incipient  ion- 
plasma  resonance  is  evident  in  the  vicinity  of  *  I  . 

For  /^>>  /  ,  S  approaches  the  "strong  coupling"  form  (4.108). 

In  Figure  4.2  we  display  the  normalized  scattering  function,  as 
in  Figure  4.1,  for  1,18  as  a  function  of  the  dimensionless  collision 
parameters  ,  Y  .  For  convenience  we  have  assumed  ; 

from  (4.96)  this  implies  equal  electron-neutral  and  ion-neutral  colli¬ 
sion  cross-sections,  i  .e.  a‘<,'vC  -  As  Y  increases  relative  to  z/'1 

collective  effects  are  seen  to  disappear;  the  resonance  near  —  euffL 
becomes  lower  and  broader  as  per  (4.113).  When^f»/  ,  $  approaches 
the  form  (4.115).  Figure  4.3  is  similar  to  Figure  4.2,  only  for 

To  understand  the  behavior  of  the  scattering  function  when 
»  we  have  computed  5  for  a  few  values  of  §  with  X~  /• 1 6  , 
when  (i)X^e^O,3tX^Y-0,t»  and  when  ( ii )/^e=0.  / ,  The 

results  are  displayed  below  in  tabular  form.  It  is  evident  that 
changes  in  the  ion  collision  parameter  have  only  a  very  small  effect  on 
5  in  the  vicinity  cf  f  —  /  .while  changes  in  the  electron  parameter 
similarly  do  not  significantly  disturb  S  for  f  small. 


h 7*»0 
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Figure  4.2  Electron  scattering  functions  for  a  singly  ionized  carbon 
plasma,  showing  the  effect  of  neutrals  for  l  =  (2*Aj))  ■  2.55.  The 

results  have  been  normalized  to  the  ideal  electron  gas  scattering  func¬ 
tion  at  5  =  ^u/oipe  -  0.  The  parameter  rje  is  defined  by  T)e  =  vne/oipe 
and  it  is  assumed  that  v06  -  (mi/mg)1/2^1.  The  broken  line  represents 
the  approximate  form  (eqn.  (4.115))  Ter  lr\e  -  10. 
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Figure  4.5  Electron  scattering  functions  for  a  singly  Ionized  carbon 
plasma,  showing  the  effect  of  neutrals  for  l  =  (2kAd)-1  *  1. l8.  The 
results  have  been  normalized  to  the  ideal  electron  gas  scattering  func¬ 
tion  at  5  »  ^>/<Upe  ■  0.  The  parameter  rje  is  defined  by  T)e  =  vne/ajpe* 
The  broken  line  represents  the  approximate  form  (eqn.  (4.115))  for 
JTie  ■  10. 
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Table  4.1 

Comparison  of  Scattering  Functions  for 
Different  Collision  Parameters 


10"3 

cn 

X 

o 

1 

ro 

1 

Sfo.  1,0,3) 

S  f°- l,o.l) 

1.21 

1.14 

0.986 

3  fo.l,o.S) 

Sfo.3,0.3) 

0.992 

0.85 

1.23 

S^>.  3,o.l) 

S(o. ) ,  o,  i) 

0.998 

1.09 

0.81 

Sfo.3,0.3) 

0.83 

1.12 

0.989 

This  is  as  we  would  expect;  when 


M  J'f 


^>1 


where  “^S^V^he  function  is  well  approximated  by  the  leading 
terms  in  its  asymptotic  form  (4.105)  and  electron  dynamics  dominate  ion 
dynamics  in  the  scattering  function.  Similarly  when 


Jt't  ~  «i 

ion  dynamics  play  the  dor-  ant  role.  These  observations  are  reflected 
in  the  approximate  forms  (4.110),  4,112),  and  (4.113)  of  the  scattering 
function. 

In  most  experimental  situations  the  scattered  spectrum  is 
observed  as  a  function  of  the  shift  Al  in  wavelength  from  the  wave¬ 
length  of  the  incident  photons.  Consider  for  example  the  differential 
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photon  scattering  cross-section  as  a  function  of  AX  for  ruby  laser 
photons  (/\B  6?*/3  A  ),,  for  a  singly  ionized  carbon  plasma  with 
h«»  lol6*sl  and 

/  =  (ift  fT(»*  ^  &)-'+  o(f)\(t.m) 

From  (4.1)  we  have,  converting  to  wavelengths, 

(li-4$-y*.,~,Sfc/H)rr&  (4.117) 

where  AX  -/A/— A  •  It  is  clear  from  (4.116)  and  (4.117)  that  ret  'r- 
ing  the  terms  in  AX/X  yieV'  a  cross-sec: ion  ch  is  asymmetric  about 
AX^O  .  for  constant  K  . 

Comparing  (4.116)  with  the  scattering  function  as  displayed  in 
Figure  4.1,  it  is  evident  that,  for  constant  K  ,  retaining  ^X/A  in 
(4,116)  has  the  effect  of  increasing  the  height  of  the  peak  at  £=  +  / 
and  lowering  the  peak  at  f «  - /  .  In  addition,  the  upper  peak  is 
shifted  nearer  AA  =  °while  the  lower  peak  is  shifted  further  away. 

Since  it  is  not  the  scattering  function  but  the  cross-section  that  is 

**> 

measured,  however,  we  must  also  account  for  the  factor  (I+-A X/X  )"J 
in  (4.117).  Clearly  this  factor  will  contribute  to  a  lowering  of  the 
upper  peak  and  an  increase  in  height  of  the  lower  peak.  The  two  effects 
are  thus  competitive. 

For  the  example  at  hand  we  have  found  the  net  effect  to  be  a 
5.9  per  cent  decrease  in  amplitude  of  the  upper  resonance,  and  an 
equal  increase  in  amplitude  of  the  lower  resonance.  In  addition,  the 
location  of  the  upper  resonance  is  shifted  about  1.3  percent  closer  to 
AX^O  while  the  lower  resonance  is  shifted  the  same  amount  further 
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away.  i_  \ 

We  have  so  far  neglected  contributions  to  the  scattered  photon 
intensity  due  to  elastic  (Rayleigh)  scattering  of  photons  from  the 
neutral  atoms.  For  simplicity  we  characterize  the  spectrum  of  neutral- 
scattered  photons  by  the  ideal  gas,  or  Doppler  scattering  function. 

The  ratio  of  light  intensity  scattered  by  neutrals  to  that  scattered  by 
electrons  is  then 

nn* r1  (Tne)l/letp  'f'Z&J 


•v 


(T^  ne  S  (£, 


I  Jt*  I 

r 1 


£L^ 
crj-  he 


(4.118) 


where  $  is  the  electron  scattering  function  discussed  above, 

and 


CT^  ^  )o~Z8 


17 


is  the  Rayleigh  cross-section.  It  is  evident  that  scattering  from 
neutrals  is  unimportant  when 


he  l  i/V 


«  [O 


(4.119) 


Moreover  when  scattering  from  neutrals  is  significant,  it  is  cl^ar  that 
the  effect  will  be  manifest  in  the  observed  spectrum  only  in  the 
vicinity  of  the  central  ion  peak;  i.e,,  when  lkco«co^t. 

From  the  foregoing  discussion  it  is  evident  that  the  scattered 
photon  spectrum  will  contain  electron  plasma  wings  whenever 
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A  5  ('ZfcAp)"1  *  I ) 


X 


*  5.  lx  io 


In  addition,  collisions  with  neutrals  eliminate  the  wings  when 


X^  »  A 


oj 


t* 


*  I  . 


It  thus  follows,  with  (4.96),  that  the  wings  should  be  observable 
whenever 


4.  Discussion 

In  the  preceding  analysis  we  have  employed  a  treatment  based 

on  a  temporally  irreversible  kinetic  theory  to  describe  the  effect  of 

neutral  atoms  on  the  spectrum  of  photons  scattered  from  electron 

density  fluctuations  in  a  partially  ionized  qas.  Our  analysis  was 

12  13 

quite  similar  to  that  suggested  by  Yip  et  al.  *  for  a  description  of 
neutral  particle  (i.e.  photons  or  neutrons)  scattering  from  moderately 
dense  neutral  gases.  We  remarked  in  Section  2  that  an  analysis  based 
on  a  reversible  kinetic  theory  predicted  neutral  atom  effects  only  at 
unusually  high  neutral  densities. 

18 

A  recent  experiment  reported  by  Greytak  and  Benedek  provides 
striking  quantitative  support  for  the  irreversible  treatment.  These 
authors  observed  the  spectrum  of  6328 ^  photons  scattered  from  thermal 
fluctuations  in  neutral  gases  near  standard  temperature  and  pressure. 
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Their  observation  of  a  symmetric  pair  of  spectral  lines  located  at 

A  is  in  excellent  agreement  with  the  theoretical  predic- 
12  19 

tion  of  Yip  and  Nelkin  *  for  scattering  from  neutral  gases.  These 
resuHs  clearly  contradict  the  predictions  of  the  reversible  theory, 
i.e.,  a  Gaussian-shaped  spectrum,  at  the  relatively  low  densities 
(  n  2 involved. 

We  have  assumed  in  cur  computation  of  the  classical  scattering 

Si 

function  that  each  of  the  particle  species  could  be  characterized  by 
a  Maxwell  velocity  distribution  with  a  common  temperature  for  all.  For 
many  experiments  this  assumption  is  invalid,  and  could  lead  to  errone¬ 
ous  conclusions.  The  extension  of  the  present  work  to  allow  for 
different  component  temperatures  is  straightforward  but  adds  consider¬ 
ably  to  the  algebraic  complexity.  While  such  considerations  are  beyond 

M  20  21 

*the  sccpc  of  our  purpose  here,  we  note  that  several  authors'  *  * 

have  investigated  the  effects  of  unequal  temperatures  for  a  fully 
ionized  two-component  plasma.  They  showed  that  the  scattering  function 
for  such  a  system  c*o  be  qualitatively  different  from  that  computed 

with  a  single  temperature  model,  In  a  recent  experiment  Kronast, 

22  14 

et  al  have  employed  Salpeter's  results  in  a  measurement  of  electron 

and  ion  temperatures  in  a  theta  pinch.  For  their  particular  experi¬ 
ment  they  found  Oe  ( . 

In  our  development  of  the  electron  correlation  function 
in  Part  2,  we  tacitly  assumed  that  the  inclusion  of  close  encounters 
(collisions)  between  charged  particles  could  be  neglected.  To  lend 
support  to  this  assumption  we  employ  a  simple  Krook  model  to  estimate 
the  significance  of  Coulomb  collisions.  From  the  analysis  of  Part  2 
it  is  evident  that  these  effects  should  be  negligible  provided 
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«  i 


(4.95) 


where  now 
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(4.120) 


9  o 

with  vc  the  Spitzer  collision  frequency^  for  charged  particles. 


v>c  * 


8ttv\  &  -A 


(4.121) 


Here /*a.  is  the  reduced  mass  for  the  pair  in  question  and  we  assumed 
for  simplicity  a  singly  ionized  gas.  Combining  (4.95),  (4.120),  and 
(4.121),  the  condition  for  neglecting  Coulomb  collisions  becomes 

aMA  <<  |  (4.122) 

A 

3  23 

where  A** A0  .  Since  A  is  generally  a  very  large  number,  we 

would  generally  not  expect  Coulomb  encounters  to  be  significant  here. 

24  25 

(Ron,  Dawson  and  Oberman  and  Fante  have  recently  estimated  the 
effects  of  C(  j  1  omb  encounters  on  the  electron  scattering  function  u:.ino 
somewhat  different  analyses  than  the  simple  model  employed  here,  li.ey 
found  the  inclusion  of  these  effects  produced  a  change  in  ^  of  the 
order  of  A”1  ). 

The  principal  result  of  this  chapter  is  the  electron  scattering 
function  for  a  partially  ionized  gas,  as  discussed  and  displayed  in 
Part  3.  It  is  apparent  that  for  given  vsl^es  cf  ,7*  ,  *  *  ,  J* 

end  A  t  the  scattering  function  as  a  function  of  is  uniquely 
determined,  Cven  sc.  we  would  not  \cect  that  a  single  experiment 
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could  serve  to  measure  all  of  these  quantities  for  a  qiven  plasma. 

Our  results  apply  for  instruments  with  infinitely  sharp 
spectral  resolving  power  and  uniform  average  plasma  density  over  the 
scattering  volume.  Average  density  nonuni t^rmi ties  and  the  finite 
resolving  power  of  instruments  together  with  the  natural  width  of  the 
incident  photon  beam  will  add  to  the  width  of  the  observed  spectral 
structure  and  it  may  be  necessary  to  take  these  into  account  in  a 
given  experiment. 

In  addition  to  the  extension  of  this  work  to  allow  for  different 
componenent  temperatures,  it  would  be  most  interesting  to  consider  the 
effects  of  magnetic  fields  and  small  spatial  gradients,  for  fully  or 
partially  ionized  gases. 


APPENDIX  A.  INTEGRALS 


Let©  be  the  angle  between  £  and  .  From  tl.e  generating 
function  relation  for  Legendre  polynomials  we  have 

l£-£T‘  - 


-  ‘-'I  (r )APj(«e)  (c'tc)  m 

/•o 

and  similarly  for  c  *  c'  .  Writing  and  employing 

the  addition  theorem  for  Legendre  H  lynomials,  we  obtain  (3.13). 

To  find  (3.14)  we  use  the  relation 

[j  +■  xl  -  2xyJ tX  -  j  Jx  x  jj*-xl-2xyj  A_ 

/  r  (A2) 

-  J  )  J*  [l  *  *'- I*]] 

Combining  (Al)  and  (A2)  we  have 


I  £  'S'l  -  C  X  t  C,X-  Zcc‘ 


Hi 


(c‘±c) 


and  similarly  for  c.* .  From  the  pure  recurrence  relation  for 
Legendre  polynomials  we  have 


(A3) 


<*&  =  jfc  P/.,  (‘o*  &)  +  PJtl  (<*®) .  <A4) 
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Combining  (A3)  and  (A4)  and  then  using  the  addition  theorem  as  before 
we  find  (3.14). 

Finally  (3.15)  is  obtained  fror  (3,14)  with  ^ 

with  the  Kronecker  delta. 


APPENDIX  D.  THE  SELF-ADJOINT  PROPERTY 


Let  l|  be  either  (i)  the  differential  operator  in  (3.20)  or 
(ii)  the  integrodifferential  operator  in  (3,20)  defined  as  in 
Chapter  III,  Part  4  on  the  closed  interval  s  b  #  o<  <*.  f 

b  <  .  Let  La  be  similarly  defined  on  (a,®0).  Then  in  either  case 

(i)  or  (ii)  there  exists  a  complete  orthonormal  set  of  functions  ^ 
on  ,  generated  by  with  homogeneous  boundary  conditions  at 

and  b  J  ^ 

Let  f  and  9  be  any  nonzero  functions  square  integrable  on 
(  o)  oo),  and  consider  the  inhomogeneous  problems 


t  A  n  *  is  v+  /vV  *  g  }  (B1 ) 

with  the  same  homogeneous  boundary  conditions  at  and  fe  as  in  the 
above  homogeneous  problem.  Let  'K.A'to  so  that  A  will  not  belong  to 
the  spectrum  of  the  set  \^}  •  Then  the  problems  (Bl)  have  nontrivial 
solutions 


Now  let 
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U  -  ft  *i  **  ^ 

^  4rfk.>  ‘ 

*  -  ft 
'  l  fk-,  *'-* 

with  u.y  and  Cj-  as  in  (B2).  Multiplying  tx*.  by  g*  and  by  f  and 
then  integrating  over  ^  we  have 


(>»*»*-  fl 

V  ->A  J.  jc*l  ~  A  \ 


V**  Jt 


f2 

-/*•  j,**i  ~~  X 


and  similarly 


=  f  J  &#>.  (84) 

V  J,*-1  A 

Since  (3.24)  and  (3,25)  are  asymptotic  solutions  for  both  the  differential 

and  integrodifferential  equations,  the  limit  matrix  f*  is  the  same  in 

either  case.  Taking  in  (B3)  and  (B4)  we  have  /°azz  and 


uy  f*  a*  /  ^ 

j  j  Sy-Wtot  =  y(*w/(c)jc. 

O  -/A  -A-  ~  ^  O 

After  taking  <a -*°o  and  employing  (Bl),  (B5)  becomes 
OO  .ac 

\(LL\/fiA^c  -  V  (La*)u*Jc 

Vi  •  A 


which  is  the  desired  result. 


APPENDIX  C.  PROOF  OF  THE  EXPANSION  THFOREM 


This  proof  is  based  largely  on  the  so-called  Weyl-Stone- 
Titchmarsh-Kodaira-Levinson  theorem,  as  outlined  by  YosidaJ  and  by 

9 

Coddington  and  Levinson,  '  Our  goal  is  an  expansion  theorem  for  real 

.00  x 

valued  continuous  functions  in  (0,00)  with  [Je  «  <*>  . 

The  expansion  functions  are  to  be  the  solutions  of  the  linear  integro- 
differential  equation 

(q+a')*+ 

Here  P*  P' ,  and  /<*  Kj  are  real  and  continuous  and  ?>o  ,  (0,00), 

,00  /OO 

and  5  )  deje' lKfet')l  <oo  ,  The  function  is  real  and  con- 

0  *Q 

tinuous,  is  regular  at  c^O  fo  r  X-O  ,  and  has  a  regular  singularity 
at  c*o  for  A  ^  o  . 

An  expansion  theorem  based  on  the  solutions  of 

* 

f  f  Kfr,*')  W)  =  °  (C1) 

'a. 

which  satisfy  homogeneous  boundary  conditions  at  o<a.  ,  fc<oo  has  been 
established  by  J.  D.  Tamarkin.  Our  task  here  is  to  extend  the  interval 
(  <k)  b  )  to  the  interval  (0,0c  ). 

Preliminaries,  Tamarkin's  Results 

Let  l £  represent  the  linear  integrodiffv-rential  operator  in 
(Cl).  The  following  properties  have  been  demonstrated  by  Tamarkin.'* 

(i)  For  a  fixed  complex  A  let  i  •  K  represent  a  pair  of 


I 
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linearly  independent  solutions  of  /.^rs-A^  »  real  for  ■  aal  A  »  satis¬ 


fying  the  conditions  ( 


4,  ■=■  i 

P/'s)  -o 


fa?sth)  »e 
tfs)fa '&,*)  =  I  . 


Then  4>,  ,  <^v  .  d>/  ,  fa?  are  entire  in  A  for  every  fixed  c  on  i  . 
(ii)  For  the  self-adjoint  boundary  value  problem 

L{  <j>  -  -  A<f> 

Co4*<I>(ck)  t  ?(<*)  <A7*0  *  O 

t  3 

there  exists  a  sequence  of  real  eigenvalues  jA^  and  a  complete  ortho¬ 
normal  set  of  eigenfunctions  In  terms  of  these  functions  the 

expansion  formula  for  any  «a  t  £VO  is 

*(*)  “  ta/O  (C4) 

*  /Ov 

To  extend  the  interval  to  Ufa0*)  we  proceed  as 

follows.  Since  t,  •  fa  form  a  basis  for  the  solutions  of  Lsf—ht 
we  can  write 

K  ft  =  rT-i  4  (C>*6»)  +■  rW  fa  <Vfw)  (C5) 

where  fj  t  are  complex  constants.  With  (C5),  (C4)  becomes 

rwC  m 

"  j,kr|  JCK 

Now  define  an  Hermitian,  positive  semidefinite  matrix  ,  called  the 


■Mmilttuifrt-iliHiMiiJKUUtiUfl'irtlHHtW  ftUtKHHVt.JiUi  iHri'  ifU  i  ■£:|i1|i:bi^-  •  n* » **! *1**!? 
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spectral  matrix,  with  elements  consisting  of  step  functions  with 

jumps  at  the  eigenvalues  A ^  given  by 

fth  ( 'S  to)  -  f>S)k  fah *»)  '  S'  C  . 

We  let  (o)  be  the  zero  matrix  and  define  away  from  the  eigen¬ 
values  by  *  fifA).  We  employ  the  spectral  matrix  to  replace 

the  infinite  series  in  (C6)  by  a  Lebesgue-Stieltjes  integral; 

(  Z  (C7: 

_  OC  ifM*l 


where 


Xik  "  ^ 


As  i  ♦  A  (that  Is,  <k.+o  ,  ^  approaches  a  limit 

matrix  .  Our  task  to  to  find  the  matrix  fA  and  to  prove  the  con¬ 
vergence  of  the  expansion  (C7),  (C8)  in  the  limit. 

We.yl's  Limit  Point  and  Limit  Circle  Theory 

For  any  number  >v»b  ,  the  expression  YtJ-4’+v*\, ^satisfies 
equation  (Cl).  We  now  choose  so  that  satisfies  the  boundary 


condition 


at  the  point  b  .  Then  must  satisfy 


-  137  - 


Note  that  the  one-point  boundary  condition  (C9)  has  not  restricted  A 
to  real  values.  Since  ^  ,  and  are  all  entire  functions 

of  A  » /A)  is  a  meromorphic  function  of  A  .  it  is  readily  shown 
that  every  zero  of  the  entire  function  is  real,  and  hence  all  the 
poles  of  w\/A)lie  on  the  real  axis  of  the  A -plane.  Consider  as  a 
function  of  A  ,  b  ,  and  .  If  we  let  «-*  ttf/b  and  maintain  A  and  b 
at  fixed  values,  we  can  write  as 

Ai  +  B 


C-i  t  D 


(CIO) 


Since 


|  AD  -  6 Cl  =  l  ^A,AWtY^A)-^A;A)^/^A)/PA) 
*  Vt'b  (h)^)  ¥  O 


where  is  the  Wronskian2  of  4,  •  ft  evaluated  at  b  ,  the 

transformation  (CIO)  is  a  one-to-one  conformal  mapping  which  transforms 
the  real  axis  of  the  £  -plane  into  a  circle  in  the  comp!  v*  /A) 
plane.  Therefore  if  V*  A  *  /  o  ,  then  varies  on  the  circle 

Ch(A)  with  a  finite  radius,  as  £  varies  over  the  real  axis  of  the  * 
plane. 

The  equation  of  the  image  of  the  real  axis,  .  is 

found  from  (CIO)*, 

or  ( Cl  1 ) 
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which  is  the  equation  for  (fO  .  It  follows  easily  that  the  center 
of  Ck  is  A*  , 


A 


9 


and  the  radius  is 

fb  (A)  - 


.  w.  g ,  a*) 

K  M,  if) 


j  wl  a, ,  <k)  I 
I  wk  di'i*) 


(C12) 


For  the  moment  let  satisfy  A^*  Lf 

with  .  Then  the  symmetry  of  L{  permits  demonstration  of  the 

Greens  formula, 

S 

Now  with  h\  -  A  and  by  virtue  of  (C2), 


Ws  (^i)  ^t)  ~  (4})$%)  ~  I-  -C13^ 

Further,  with  (pt  (c}A*)=  (p*fc}A)t  $*&/  A)  ,  and  making  use 

of  Greens  formula,  we  have 

b  /b 

I  p1  (c  A)l*dc  -  7co\  <h(‘,*)<h&)rf)<}c 

S  Sr_ 

-  i  - 

-  I  K  A  *),  <#,  **)] 


Combining  (C12),  (C13),  and  (C14)  we  obtain 
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[sm  isb  i  ’hfoA'ii'y*]  i  )«/\ 


(C15) 


Lemma  1.  If  #  then  the  interior  of  the  circle  M  is 
mapped  onto  the  lower  half  plane  of  the  £  -plane  by  the  transformation 
(CIO). 


Proof.  Since  the  real  axis  of  the  ?  -plane  is  the  image  of  the 
circle  Ch  ( A )  by  the  transformation  (CIO),  the  interio.  of  ^  (A)  is 
mapped  onto  either  the  upper  half  plane  or  the  lower  half  plane  of  the 
?  -plane  and  further,  the  point  at  infinity  of  the  wn  -plane  is  mapped 
onto  the  point  -  of  the  *  -plane. 

On  the  other  hand,  we  can  write 


^  PA) 4>v  Ay  A/ A )J  * 


,  _  t 

1  IMbj  A)  I' 

.  b  X 

__  <0  (  |  Px  (c)*)\  C k  (C16) 

'  - -i - - - >  O. 

This  means  that  -  f fa) 4%  fati/Pi/l/S)  belongs  to  the  upper  half  plane  of 
the  $  -plane.  Hence  the  point  at  infinity,  which  is  not  contained  in 
the  interior  of  c\,  ( A),  is  mapped  into  the  upper  half  plane.  This 
proves  the  lemma. 

Since  the  transformation  (CIO)  has  an  unique 

inverse  which  is  given  by 


p  A  \  Pi  (ti)  t- 

fkAA^M+^/V) 


(Cl  7) 


In  view  of  Lemma  1,  if  \wv  /l  =  £0  >0,  w\  belongs  to  the  interior  of 
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the  circle  /V\)  if  and  only  if  V*\^<o  ,  namely,  1(2 -g4*)  >0. 
From  (Cl 7)  it  follows  that 

Hi-t*)  ^  I  r_  p/'t)  liM”  i  HIM  f 
L  *■  4  <V> 

+  p/t)  ^/7fa/^)w<,f^/7k//')  -1 

+  **  faA)  J 

=  w,  4, *■*-<-*  4*) 

I  4,  *■  <*  <hll  ’ 

Therefore,  L,  2<0  if  and  only  if 


t  ^  i  m  4*  t  v*  >  0  .  (C18) 

By  Green's  formula  we  have 

2 so  (  I  d  £  — 

-  *  ft  ($,  +  >*  tij  4>, *+”<*£*)-  (p,  +  ^  ( ci 9 ) 

and  from  (r2), 

4-  (P^P,*)** 

1~  ($)&*)***  +  V^i  (4i)4i*)ImI% 


—  -  r  »w  * 

“  ”  ?  t  ^  w  , 


(C20) 


Combining  (CIS,  (C19),  and  (C20)  we  have 


Letnma__2.  If  do  *  A>Ot  then  wa  is  interior  to  the  circle  (a)  if 
and  only  if 
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I  1$,  fc,A)  ii(c  A)l  Ac  < 

3 

and  w  lies  on  Cw  if  and  only  if 

/*>  2 

)  1  <f,  f  IM  ^/<y  A)  I  Jc 

$ 

(Note:  it  is  easily  shown  that  Lemma  2  also  holds  when  co  *  ^  A<0 ,) 
It  follows  that,  if  w\  is  inside  ,  and  k'<b  f  then 

<W 

\  )  +  **  ^v|  Ac  < 

Hence  m  is  also  in  C^,  ,  even  though  the  centers  of  and  Cb/ 

may  not  coincide.  We  thus  have 

Lemma  3.  The  circle  C^,  contains  C\,  for  k  '<!>  and  A  ?o. 

It  follows  that,  as  b^oo  ,  the  circles  converge  either  to 
a  limit-circle  or  to  a  limit-print.  In  the  limit  circle  case  we  have 
from  (Cl 5)  that  is  lx(s.=o)  ;  the  same  property  is  readily  demonstra¬ 
ted  for  . 

Next  consider  the  boundary  point  ck.  ,  occl<  £>  .  For  an 
arbitrary  real  number  o(  ,  the  boundary  condition 

Co4  (A  Y'o^  +*  W\  ot  ?(°)  Y'n 
at  the  point  a.  with  fa  ,  determines 

M  /a)  -  -  ^  4/m)  +  PAM'fc.'O 
*  PfrWM) 

and  also  the  circle  CvVV)  ,  described  by  the  equation 

Wk(%^  K*)  »  O  . 


fit?  5 


(C21 ) 
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Similarly  as  in  Lemma  2  we  can  prove 


Lemma  4.  If  co  *  A  Jo  ,  then  v*  lies  on  the  circle  or  in  its 

interior  depending  on  whether 


CO 


or 


)<K 


Similarly  as  in  Lemma  3  we  can  prove 


Lemma  5.  The  circle  contains  Q  for  a><a‘  and  Va  hJo  . 

As  before,  asd.-*o  the  circles  converge  either  to  a  limit 
circle  or  to  a  limit- point. 


The  Limit  Matrix 

Recall  that 

hi  -  h(?t)  +Qi  + 

Let  </(  ,  <pr  be  solutions  of  satisfying  the  conditions 

cj>  Y s,  /i1)  =0  p/s)  4YY  a)  =  i . 

For  the  self-adjoint  boundary  value  problem  (  <><<*.  ,  y<o&) 

If  ^  »  -  A  ^ 

4*  6(  <^/a)  f  Pfa)  (f  'fa)  -  O 

t  (C/6  P/t>)  ^ Y i>)  =0 


(C3) 
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there  exists  a  sequence  of  real  eigenvalues  and  a  complete  ortho¬ 
normal  set  of  eigenfunctions  .  The  expansion  formula  for  any 

function  tA  £  fVi)  is 


(C4) 


Multiplying  (C4)  by  u? (<)  and  integrating  over  &  gives  the  Parseval 
relation 


CluMl'Jc  -  Z  I  ) 

'a  VA  1  A. 


(C22) 


Similarly  if  o\r fO;  (c)  e  then 

b  ,  b 


^  Jc  Uj  (c)  -  ^  (c \c  <A^  (c)  ^C23^ 

(K  V,  J ^  L- 


Employing  the  representation  (C5)  and  the  spectral  matrix  ^  ,  (C22) 
may  be  rewritten  as 


j  ju^fc/c  -  ^  ^  «( (A) 


(C24) 


__  oo 


where 


*  ^  Ac  t a(c)  <pk  (t}\)  . 


Applying  the  Parseval  relation  (C22)  to  any  continuous  function 
ia  on  (o^)  which  vanishes  outside  some  interval  S}  ,  contained  in 
£  ,  one  : btains 

\  ih&h'J*  -  1  Z  M  (c25) 

°  -OO  J  irss) 


where 
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=  (  6c  (ctX)  . 

Let  be  a  solution  of  L*4>  *- Aj> ,  ^  f\  <to 

satisfying  the  boundary  condition 


i"  <*  PA)  <V(°)  ~  0 

and  similarly  let  ^  be  a  solution  of  the  same  equation 

satisfying 


44yS  ^  A)  PA)  <^A)  -  °  . 

Then  and  lie  on  circles  c ^  and  in  the  complex  w  -plane 
whose  equations  are,  respectively. 


K  ( X,  %*)  =o,  wb  (\  r;)  -o  .  (C26> 

It  Is  easily  shown ^  that  Green's  function  for  the  boundary  value 
problem  (C3)  is  (  0  ) 


(  '1-Jt'A  c&c, 

)  >%//0  ~  V*w/V\} 

/  ^<k(£,)A)^(CJ^  c>  c, 

where  wi^  (a)  -  w\b  //\)  =  The  Parseval  relation  in  the  form 

(C23)  is  now  applied  to  the  functions 


6C  (<,c'}'') 


iti-y 


?Jb 


°>  1 


yielding 
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fliti  t£*  j  -  T 
).  y<-)‘ '  4- 


b  . 


v\ 


£f  t «  iff  Je  I  (C27) 

i'K  K{“  ^  L  L  J  • 


From  the  definition  of  Gr  it  follows  that 


(V;  A )  -  ^  "~7~v  c  -  S 

Hv/a)  -wib 


(C28) 


?bM 


<r  >  S 


/ 

cf*% 


b  / A ) 


Pfe)  £*v/W  -  Mb/A\7 


S 


A  fa  /\^  ^  >  5 

P/i)  l»*o,  (A)-  nh{A)J 


(C29) 


Using  (C28)  and  (C29)  and  Green's  formula,  the  integrals  in  (C27)  can  be 
evaluated.  For  example 

f  [ihh*)p‘] 

-  I»y4)_wk4)l  4  )  - 

-  W*fa(A*)j 

_  ' )  -  w,4)^ 

I  *V  /a)  -  /a)  I L 


(where  we  have  made  use  of  (C26)  to  arrive  at  the  second  step).  There¬ 
fore 
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«  V  [*^A>  -  *» 


(C30) 


Similarly 


(A)  -  fyjj  (A~X^S)  j  6f(c>s;  (<)dc  * 

=  w,  (K,  O  -  w$  (ytk)  C  )  t-  W-s  <%,  K*)  . 

-  w'JMft) 

=  Ws  [^-^),  J.fj  -  j/^(A,  kl) 

w  [»k  A).  A)}  r£* 


and  hence 


I  &c(eis,A)  k*„  (e)Je 


Tt». 

z^.-A 


(C31) 


In  arriving  at  (C31)  we  have  used  the  fact  that  ^4  ^4,  lf*)-4>which 
follows  since  both  X^  and  satisfy  the  same  boundary  condition  at 
^  ;  similarly  . 

Recalling  the  definition  of  the  spectral  function  we  combine 
(C27) ,  (C30) ,  and  (C31); 


M 


W  M(„  A) 


where 


(C32) 


We  can  similarly  show  that 
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f  dAik  (A)  _  U 

L  I  a- A  I1 


where  is  given  by  (C32)  and 

tofu  f/\)  =  M^2|  (A)  -  2 


/A )  j  /A) 
wiA  /h)  M  ; 


(C33) 


/a)  *k  /a) 

'  ill  ~  w^/a)  * 

From  Lemma  2  and  Lenina  4  we  have 

AJ  (  i  /C^  ^A]\l(}c  ■  /a)  , 

*S 

Thus  for  a  fixed  /I  ,  W /I  ny/A)  and  w>k /a)  are  in  oppostie  half 

planes.  Suppose  A *£  in  (C33).  Then  points  wy ({ )  lie  on  a  circle  ^ 
which  is  in  Ct/X  for**  j|  whereas  points  mk/V)  lie  on  Ck  which  is  in 
for  ^>3s/lt  Thus  there  is  a  constant  k(>  o  such  that 
*\  Ml*  if,  for  *|  ,  k>  ^  .  Since  /*')  and  wk  ({)  are 

uniformly  bounded  for  *  <  4.  *  t>  H  »  it  follows  from  (C33),  and  the 
definition  of  the  f'V^  that 


Thus  for  V  >  o  t 

_v 

This  last  together  with  Nk(o)  =0  gives 
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IfftM  I  <  }((l+AX)  -ao<A<°om  (C34 

We  are  now  able  to  prove  the  existence  of  the  limit  matrix  .  For 
this  we  need  the 

Hell^  Selection  Theorem.2  Let  *  *  / ,  2 , be  a  sequence  of 

real  nondecreasing  functions  -  ~<A<  ~  ,  and  iet  \i (A)  be  a  continu¬ 
ous  nonnegative  function  on  the  same  interval.  I  ' 

IkWl  *  Hfl) 

then  there  exists  a  subsequence  and  a  non-decreasing  function  k 

such  that 

|  h(A)\  £  tffo)  -oo<A<<^ 

and 

u. 

It  follows  from  (C34)  and  the  Helly  Selection  Theorem  that  there  exists 
a  sequence  of  intervals  ^  (Oj0*)  and  corresponding 

boundary  conditions  prescribed  by  ,  fa  ,  such  that  (A)  tends 
to  a  limit  It  is  easily  seen  that  the  limit  matrix 

/°A  »  Hke  fit  is  Hermitian,  is  positive  semi-definite,  and  is  of 
bounded  total  variation  on  every  finite  X  interval. 

It  remains  to  establish  an  explicit  formulation  for  f>^  .  From 
(C33)  we  have,  with  and  £  &  t 

b  M *  ( *  . 

co  )  (/*-Ar 

_  OC 
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Let  A,  ,  Ax  be  points  of  continuity  of  .  Then  integrating  the 
above  witiwj^O  held  fixed  and  finally  taking  co-^o  ,  we  have 

w  r  / 


A  |  -0* 


&  I  iV  (x-^)-^',(h^H 

_  Otf 


-  TT 


where  is  given  by 

Ma.i  1 


I 


Moo/*)  - 


h/]  M  -  1  ***>  /^)  +*  ^ 

Mi,Z  '  %u  *  0)^5) 


M 


*u/a)  a*,  A) 
avl  ~  v*«,  /ft)  -  vm0^a) 


If  both  points  <k  ♦o  ,  fe-^o©  are  in  the  limit;  point  case,  w# 
and  ^  are  unique  and  it  follows  easily  that  fe  is  unique.  If  either 
point  is  on  a  limit  circle,  the  spectral  matrix  is  not  unique  without 
the  specification  of  a  boundary  condition  at  the  point  in  question. 
Whether  a  particular  case  is  limit  point  or  limit  circle  is  readily 
determined  from  the  asymptotic  solution  and  the  expression  (Cl 5)  for  the 
radius  of  Ck  or  its  analog  for  CK  . 


The  Parseval  Relation  in  the  Limit  J-^A 


Consider  a  function  t*/0  having  a  continuous  second  derivative 
on  o<C<<x?  ,  and  which  vanishes  outside  some  interval  <T,  ,  contained 
in  $  .  Then  applying  (C25)  to  LfU.  we  have 
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[  lif«l  L  >  £  ^  f  Je  jTj  Jfa'  (C35) 


Applying  Green's  formula, 

00 


^  **  )  d  c  »  J  ( Lg  <j>j)  L\dc 

=  ~  ^  tyu  Jc 

>o 

“  “  ^  Mr/  ^  ^  ) 

and  hence  (C35)  becomes 


j  |l,«fje  -  \  ^  1.  %k*tj  A  > . 


(C36) 


J>»l 


Now  for  large, 

-“A  ,M,  % 

'  )z 

'  t  '  *  1. 


(f 


this  last  following  from  (C36) 

It  is  convenient  to  rewrite  (C25)  in  the  fonn 

P*1'"  ■  (£*C*0 1.“*  +v,  ■ 


(C25) 


J.Jcn 


Combining  (C25)  and  (C37) 
oo  /  *  A 


-  \  2  “ft  j  *  jl*.  ( |L f 

-A  Jl 


(C38) 
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To  take  the  limit  «£•*■  A  (that  is,  ,  b ■*►<»)  in  (C38)  we  need 

the  following 

Integration  Theorem.^  Suppose  J^KA)]1S  a  real,  uniformly  bounded, 
sequence  of  nondecreasing  functions  on  a  finite  interval  J^A-c  ,  and 
assume 

d-A-e  . 

W  ->-o0 

If  "P  is  any  continuous  function  on  (J*A-e.  ),  then 

We  established  earlier  the  properties  of  required  by  the 

integration  theorem.  Thus  letting  through  the  sequence  of  inter¬ 

vals  found  above,  it  follows,  using  (C38)  and  the  integration 
theorem,  that 


Now  allowing  /}->oa  ,  there  results  the  Parseval  equality 


00  » 

ll*M \%dc  *  \  7  (XkXj<lfak  (C39) 

»  I 

for  any  u\(c)  restricted  as  above.  We  new  show  that  the  Parseva 
equality  holds  for  any  t*(c)  in  First  suppose 

and  vanishes  for  c.  sufficiently  large  and  sufficiently  small.  Then 
there  exists  a  sequence  of  functions  14 h  t'  tX(°i°°)  possessing  continuous 
second  derivatives  and  vanishing  near  0-0  and  for  all  large  c  such 
that 
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'0s' 

i  Nk-  c/ =  o. 


Applying  (C39)  to  , 


i  |i*w-M*|VcU  *  (  2  . 

b  Joo  <k*l 


( C40) 


Since  the  left  aide  of  (C40)  tends  to  zero  as  *  ,  it  follows  that 


the  sequence  of  vectors  ,  ***  *  *7*)  *V-  where 

ao 

J  c  <v 

'  »  A 


(C41) 


converges  in  the  mean  in  .  and  since  the  latter  space  is 

2  ^ 

complete  there  exists  a  vector  ia  which  is  the  limit  in  the  mean  of 
this  sequence.  It  is  clear  from  'C41)  that  the  components  of  U  are  the 
continuous  functions 

,°o 


<Aj'  =  (  dc  u(c)  (j>j  ^A)  . 

Loo 


Returning  to  (C39) 


K  'o  L‘ohk'1 

=  (  Z  4*  4j*  J/Vk  , 

_oo  J/k-l 

which  proves  the  Parseval  relation  for  any  *(c)z£ V^00)  vanishing 
for  all  c  sufficiently  large  and  sufficiently  small.  Suppose  now 
that  t\(*)  is  any^  function  of  class  tY0/0*)  and  define 


C6")  - 


<*(c) 

O 


x-‘*3 

C<x  .  n  <c 


)  3 
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and 


a/ 

(A 


.00 


xcjj  ft  )  “  j  “xy  <pj  Jc 

/» 


'  )  frjA)  ole 

X 

and  similarly  for  Mrs(^:)  .  Since  {x<r)  ^  <  5  ) 

z  -  KsK^^k  -  ( [ + 1  j  h^jt) 


,*>  i 


-Oo  v)j  X  ~  I 


it  follows  that  the  set  of  vectors  t\x^  converges  (as  a  -+o  ■*<*>  ) 

in  the  mean  <f  $^A)  to  a  vector  function  ia  f  /  (fA)  .  By  letting 
X-=3>  o  ,  ^  -^00  in 

!  t  v  ^9/  =  SxK‘),lJ‘ 

- <0  J,k*l  * 


there  now  follows  the  Parseva!  equality 
.00  v 


i<«  Jy*  =/ 

for  any  tx,  ^<0  f  £  Y° ) 


v  °o 

Z  #k  Kj*  d/>/>4  M  ,  ( C42) 

I*  1/  —I  'a 


The  Expansion  Theorem  f -<r  the  Singula;  Interval 


With  the  Parseval  relation  established,  i;he  proof  of  the  expan¬ 
sion  theorem  may  now  be  given.  Let  D  1  hn^  and  define 

/7  1* 

«» r«)  ■ )  I  _  k  A)  ^).  <c«> 

_  ^  4jX~l 

If  and  nz(t)  are  in  ^Y0/00)  ,  then  the  relation  (C42)  implies 
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,oo 


,*>  x 


0<? 


Ux&toZ&dc  =  \  J  (  c/c  j  «jf  4  d< 

O  -.<*5  J)k~l  lg  'O 


(C44) 


which  follows  since  we  can  write 


-  j  *4j  +  |  + 

Now  consider  some  function  which  vanishes  for  cr  <  cr (  ^ 

AS 

<rt<e  ,  and  represent  the  transform  of  J*  by  the  vector  T  .  Multi- 

"rtf' 

plying  (C43)  by  /  and  integrating  we  have 

C  J<  *  C(  (7Z 

H  V  ”*7  I,k'l 


'  -1  f|k*l 
,14  ' 


Z  i*^') df 4Jk  A) 

-1  J'l*'-' 

-  r?fc 

J'Vj 

From  (C44)  for  m x  *  M  and  Mjj  -T  , 

*)  Yd/“Wk  .  * 

£  |  _oo  Jji*l 

Subtracting  (C45)  from  (C46),  and  using  the  Schwarz  inequality, 

i  o-«.>T*j*f  *  0^ 

1  ^  W;/(*l  '7  y,^=/ 

- 

/  /*?  \  ^  (  ~  'V-  I 

Nk 

=  (£- 


yAft 
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Applying  this  inequality  to  the  function  T(^)  given  by 


o<c‘l  6^ 


we  obtain 

(  /  -  (  i  ”  ^  )  Z^"  *** 

C(  1 

or  finally 

*  (C-  C)  i  <  «c  . 

v'  -t 

Since  the  right  side  does  not  depend  on  c,  ,  ,  the  above  holds 

with  ex~*> o  ,  Letting  ^  yields  the  expansion 

^  -  (  Z  A)*/w*  A) 

_loo  J;  *=» 

which  clearly  converges  in  the  mean  in 
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